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Abstract 

Supersymmetry (SUSY) remains one of the leading candidates for physics beyond the Standard 
Model, and the search for SUSY will be a central focus of future collider experiments. Comple- 
mentary information on the viability and character of SUSY can be obtained via the analysis of 
precision electroweak measurements. In this review, we discuss the prospective implications for 
SUSY of present and future precision studies at low energy. 
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I. INTRODUCTION 



With the advent of the Large Hadron Collider (LHC), the search for physics beyond the 
Standard Model (SM) will enter a new era. While the establishment of the Standard Model 
constitutes one of the major achievements in 20th century physics, particle physicists have 
long been aware of its limitations and have looked for clues as to the larger framework in 
which it is embedded. These limitations include the number of a priori unknown param- 
eters (nineteen), the absence of any requirement of electric charge quantization, the lack 
of coupling unification at high scales, and the instability of the electroweak scale under 
radiative corrections. From the standpoint of cosmology, the SM also falls short, as it 
provides no explanation for the abundance of either the visible, baryonic matter of the 
universe or the cold dark matter. The recent observation of neutrino oscillations and the 
corresponding implication that neutrinos have small, non-vanishing masses also points to 
physics beyond the SM as originally written down by Glashow, Weinberg, and Salam [lj]. 

One of the leading candidates for the larger framework in which the SM lies is super- 
symmetry (SUSY). For over two decades, the attractive features of SUSY have inspired 
particle physicists to explore its theoretical basis and phenomenological implications with 
great vigor. These attractive features include its elegant mechanism for producing elec- 
troweak symmetry-breaking and stabilizing the electroweak scale; its generation of coupling 
unification at the grand unification scale; its introduction of a candidate for the cold dark 
matter (the lightest supersymmetric particle, or LSP); and its possibilities for explaining 
the abundance of baryonic matter. In addition, SUSY is a rather generic feature of su- 
perstring theories, so one might expect a low-energy remnant of string theory to exhibit 
features of supersymmmetry. The presence of these elements that could resolve many (but 
not all) of the shortcomings of the SM have outweighed the costs of introducing SUSY, 
such as the additional large number of a priori unknown parameters, and have inspired a 
vast literature in particle physics during the past two decades. 

It is hoped among SUSY enthusiasts that experiments at the LHC will finally un- 
cover direct evidence for low-energy SUSY, and the variety of corresponding high-energy 
signatures have been discussed extensively elsewhere 0. In this review, we focus on an- 
other frontier in the search for SUSY: the high-precision, low-energy frontier. This low- 
energy frontier, which lies at the intersection of particle physics with nuclear and atomic 
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physics, has seen substantial recent advances in both experimental techniques and theo- 
retical analysis, making precision low-energy studies a powerful probe of SUSY. Roughly 
speaking, the sensitivity required for these studies to probe supersymmetric effects is given 
by <5susy ~ (a/ir)(M/rh) 2 , where M is an appropriate standard model mass scale and m 
is a typical superpartner mass 1 . The prefactor of a/ir appears because SUSY contribu- 
tions typically first arise at one loop order and are of electroweak strength 2 . A well-known 
illustration occurs in the case of the muon anomalous magnetic moment, (g^ — 2), where 
M ~ m M and where, for rh ~ M w , one has 5susy ~ 1CT 9 . Both the precision of the latest 
experimental measurement of (g^ — 2) [3] as well as the uncertainty in the theoretical, 
SM prediction, are at this level, making this observable an important means of accessing 
SUSY loop effects. Indeed, the ~ 2a deviation from the SM prediction reported by the 
E821 collaboration, gives the first tantalizing hints of low-energy SUSY. 

While there has been considerable recent attention paid to (g^ — 2) for these reasons, 
the high energy physics community may have less appreciation for the analogous power of 
other precision, low-energy observables to provide a window on supersymmetry. In the case 
of weak decays or electroweak scattering, for example, one has M ~ M\y, making 5susy ~ 
10~ 3 rather than the 10~ 9 in the case of (g^ — 2). As we discuss in detail throughout the 
remainder of this review, both the precision of low-energy weak decay and lepton scattering 
studies - as well as that of the corresponding theoretical, SM predictions - is now at the 
1CT 3 level or better, making an analysis of their implications for SUSY a timely endeavor. 
As with the study of precision electroweak observables at the end of the last decade, the 
study of low-energy precision electroweak processes can provide important information that 
complements what we may learn from the LHC or a future, e + e~ linear collider. Indeed, 
comparisons of the value of the top quark mass implied by precision data with the value 
determined by direct top quark observation at the Tevatron provided a significant test of 
the self-consistency of the SM at the level of one-loop radiative corrections and stands as a 
major success in the history of SM. Given the level of experimental and theoretical precision 
now available for the low-energy studies discussed here, one can anticipate a similarly useful 
comparison of indirect and direct search information in the LHC era. Moreover, there exist 
special cases - such as searches for lepton flavor violation or permanent electric dipole 
moments - where the SM predictions lie well below what will be achievable in the next 
generation of precision studies, where the expectations of SUSY effects lie well within the 
reach of these experiments, and where the physics reach of the low-energy measurements 
can significantly exceed what is accessible at the LHC or a linear collider. 

In the remainder of this article, we describe the recent experimental and theoretical ad- 
vances that have led to this situation and discuss the corresponding present and prospective 
implications of precision, low-energy measurements for SUSY. In doing so, we attempt to 
provide sufficient background material for readers from both the low-energy and collider 



1 As we discuss in later in this review, exceptions occur. 

2 The SUSY loop corrections scale as a/ir rather than a/in since the SU(2)x coupling that enters one 
loop amplitudes is related to the electric charge via the weak mixing angle as gi = e/ sa\Qw ■ The value 
of sin 2 6w ~ 1/4 compensates for the factor of four that would ordinarily appear in the denominator. 
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communities to be able to appreciate these developments. We begin with a brief review 
of low-energy SUSY, and refer readers to the excellent, more extensive "primer" by S. 
Martin Q] for additional details. Because many of the SUSY effects discussed here arise at 
loop level, we also provide a brief review of renormalization as it is applied to the observ- 
ables of interest here. The bulk of our subsequent discussion involves a review of low-energy 
charged current and neutral current experiments, searches for lepton flavor and number 
violation, tests of CP-violation and the corresponding implications for cosmology. Because 
members of the high energy community may not be so familiar with the phenomenology of 
these studies, we provide some background material while referring readers to recent, com- 
prehensive studies of low-energy precision tests j^]. In addition, the information on SUSY 
obtained from high energy studies is important in analyzing the low-energy sector, so we 
also give brief summaries of the present implications of high energy experiments (for recent 
reviews, see e.g., Ref. 0,0]). Finally, the reader will notice one significant, but intentional, 
omission from this review: a discussion of the present situation regarding (g M — 2). Because 
the recent literature on this topic is so vast and because there exist useful, recent reviews 
(see, e.g., Refs. 0, 0, S]), we believe that a truncated discussion in this article would be 
redundant and would not do justice to this important measurement. Thus, we refer the 
reader to the literature for a proper review of the muon anomalous magnetic moment, and 
concentrate on the other precision tests in the remainder of this article. 



II. MINIMAL SUPERSYMMETRIC EXTENSION OF STANDARD MODEL 



A. Introduction 



The Standard Model of elementary particle physics has been confirmed to high precision 
by a wide array of experiments. The strong, weak and electromagnetic interactions are 
described by SU(3) C x SU(2) L x U(l) y gauge interactions. At low energies, SU(2) L x U(l) y 
is broken to U(1) EM symmetry by the Higgs mechanism, which generates masses for the 
W ± and Z bosons, while keeping the photon massless. For this purpose, a complex scalar 
SU(2) L Higgs doublet H = (H + ,H°) is introduced. The electroweak symmetry is broken 
spontaneously when the neutral component H° gets a vacuum expectation value (VEV): 
(H°) = v/y/2. While three de grees of freedom for the Higgs doublet are eaten by W ± and 
Z (corresponding to their longitudinal degrees of freedom) , one physical Higgs boson, h, 
which is the real part of H°, remains in the spectrum. For the Higgs potential V(H) = 
m 2 H H^ H + \\W H\ 2 , the mass of the physical Higgs is related to the Higgs quartic coupling 
A and Higgs VEV v = 246 GeV: ml = 2Xv 2 . For A of order unity, the Higgs mass is around 
the electroweak scale of a few hundred GeV. 

Being a fundamental scalar particle, the Higgs boson can receive large corrections to 
its mass from quantum loop effects. Assuming the SM is an effective theory valid below 
a cut-off scale Ayv, the Higgs mass, m h , depends strongly on physics at the scale Auv- 
For example, any SM fermion / with Yukawa interaction (A/ /y/2)hff induces a one loop 
correction to the squared mass of the Higgs. The leading contribution to the mass of the 
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physical Higgs depends quadratically on A 



uv 



Am ^ = ufe [~ 2A?JV + 6m ' ln ( A uvM) + •••]. (!) 

where Auv is the cutoff scale, which could be as large as the Planck scale M p \ = 
(87rGNcwton) -1 ^ 2 = 2.4 x 10 18 GeV. A precise cancellation of 32 orders of magnitude between 
the tree level bare Higgs mass and the radiative corrections is needed to obtain a physical 
Higgs mass around electroweak scale. Such a high level of fine tuning is usually refered to 
as the "hierarchy problem" 3 . Finding a solution to the hierarchy problem points to new 
)hysics beyond the SM, such as supersymmetry 0, 10], extra dimensions little Higgs 
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[l2j |. composite Higgs ll3|, Higgsless models 14], etc. 



Supersymmetry - a symmetry under interchange of bosonic and fermionic degrees of 
freedom - is one of the most promising new physics scenarios among various proposals. For 
each particle in a supersymmetric theory, there exists a superpartner with spin differing 
by a half unit. When SUSY is exact, the masses and the gauge quantum numbers of 
superpartners are the same, and the couplings are related by the symmetry. These features 
protect the Higgs mass from receiving the problematic quadratic dependence on Auv as 
these contributions from fermionic and bosonic superpartners cancel. 

For example, the Ay V term in Eq. ([1]) from fermion / is cancelled precisely by the 
contribution from its scalar partners 4 S with mass mg and Higgs coupling As|iT| 2 |jS'| 2 [3|: 

Am 2 = [2A UV - 2m 2 ln(A uv /m|) + ...], (2) 

A cancellation of the quadratic divergence occurs when we employ the supersymmetric 
relation As = |A/| 2 and add Eqs. (0Q) and (J2J). The remaining Auv-dependence is only 
logarithmic level and fine tuning is no longer necessary. In Eq. (j2J) the logarithmic term 
proportional to m| arises from the tadpole graph containing the full quartic scalar inter- 
action. An additional logarithmic contribution to the Higgs mass arises from the diagram 
containing two insertions of the triscalar interaction \^2Xsv h\S\ 2 induced by the quartic 
interaction after electroweak symmetry breaking (EWSB). 

Am 2 = [-4m 2 m(A 2 v /m|) + ...], (3) 
The explicit dependence on to/ appears because 

X S \H\ 2 \S\ 2 = \X f \ 2 \H\ 2 \S\ 2 = \X f \ 2 + V2vh + ^ + ■ ■ ■ ) \S\ 2 (4) 
X f V2v = 2m f (5) 



3 The hierarchy problem has two elements: the large scale difference between M p \ and M wea k, and the 
need to cancel the radiative corrections to maintain a light Higgs scalar. The need for fine tuning to 
achieve this cancellation might be more accurately termed a "naturalness problem" . 

4 One Dirac fermion / has two complex scalar superpartners. Eq. ^ assumes the both scalar superpart- 
ners have the same mass mg. 
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with the "+ ■ • ■ " denoting the other scalar doublet degrees of freedom. Adding Eq. fl2j) and 
Eq. ([3]), the total contribution from a pair of complex scalars to the physical Higgs mass 
are 

Ami = T^S [ 2A uv - (2*4 + 4mJ) ln(A* v /m|) + ...]. (6) 



It is obvious from Eqs. ([I]) and ([6]) that in the supersymmetric limit, when rrif = ms, the 
logarithmic contributions from scalars and fermion also cancel each other. 

Besides providing an elegant solution to the hierarchy problem, SUSY has a variety of 
other attractive features. In the SM, the SU(3)c SU(2)x, and U(l)y gauge couplings come 
close to unifying at high scales, providing a tantalizing hint of grand unification. With TeV 
scale mass for superpartners, the SUSY (3 functions lead to coupling unification at a scale 



Mqut ~ 10 16 GeV - close to the Planck scale 15] . In addition, electroweak symmetry 
breaking can be generated radiatively with SUSY, due to the Oil) Yukawa coupling of 
top quarks and their scalar superpartners. Finally, SUSY provides viable particle physics 
solutions to problems in cosmology. In the minimal supersymmetric extension of the SM, 
for example, the lightest supersymmetric particle is a natural candidate for cold dark 
matter (CDM) if it is protected from decays into SM particles by a symmetry known as im- 
parity. Similarly, SUSY contributions to the Higgs potential and the introduction of new 
CP-violating interactions involving superpartners can make supersymmetric electroweak 
baryogenesis a viable mechanism for explaining the baryon asymmetry of the universe. 
In short, the theoretical and cosmological motivation for considering low energy SUSY is 
strong. 



B. The Minimal Supersymmetric Extension of Standard Model 

In the simplest supersymmetric extension of the SM - the Minimal Supersymmetric 
Standard Model (MSSM) - provides a useful framework for discussing the phenomenol- 
ogy of low energy SUSY. Although there is considerable interest in amplifications of the 
MSSM - particularly in the neutrino and Higgs sectors - we will concentrate on the MSSM 
throughout this article. In the MSSM, each SM particle is accompanied by a superpartner 
with the same gauge quantum numbers as given in Table [J for the matter fields and in 
Table [UJ for the gauge sector. The symbols for the SM superpartners are the same as the 
corresponding SM particles, but with a tilde on the top. For each quark and lepton, its 
spin superpartner is called a squark and slepton, respectively. The fermionic superpart- 
ner of each Higgs boson is called a Higgsino. Note that in MSSM, introduction of two 
Higgs doublets with opposite hypercharge - H u and Hd - is dictated by the requirement of 
anomaly cancellation among fermionic Higgsinos. In addition, the Yukawa interactions in 
supersymmetric models are derived from the superpotential, which must be holomorphic 
in order to be supersymmtric. In contrast to the SM, where the same Higgs gives mass 
both to the up and down type quark, the latter receive mass in the MSSM from the VEVs 
of the neutral H u and H d , respectively. Finally, the spin-1/2 superpartners of the SU(3)c, 
SU(2)i and U(l)y gauge bosons are called the gluino, Wino and Bino, respectively. 
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quark sector 


lepton sector 


Higgs sector 


SU(3) C , SU(2) L , U(l) y 


Q u d 
(3,2 i) (3,1 -§) (3,1 i) 


L e 
(1,2,-i) (1,1,1) 


(1,2,|) (1,2,-i) 


spin 
spin 1/2 


(ulJl) u* r d* R 
(u L ,d L ) u R d R 


ip,e L ) e|j 
(^e L ) 


(£+,tf°) (H° d ,H d ) 



TABLE I: Field content for the quark, lepton and Higgs sectors of the MSSM. 





gauge sector 


SU(3) C , SU(2) L , U(l) y 


(8,1,0) (1,3,0) (1,1,0) 


spin 1/2 


g W ± ,W° B° 


spin 1 


g W ± ,W° B° 



TABLE II: Field content for the gauge sector of the MSSM. 



The Lagrangian of MSSM can be written as 

C = £ gauge + ^chirai - v / 2(7[(0*T»A a + X^T a <j>)\ - i Yl OK^T^) 2 , (7) 

where i runs over the SU(3) C , SU(2) L and U(l) y gauge groups; denotes a spin-0 complex 
scalar field and tp is the corresponding fermionic superpartner; A a is the gaugino field for 
SU(3) C , SU(2) L and U(l) y , with ^ being the corresponding gauge coupling and T a is the 
hermitian matrix for the gauge group in the fundamental representation. The Lagrangian 
for the gauge fields £ gauge contains the kinetic term for gauge bosons and two-component 
gaugino spinors A a : 

C gmge = ~F^ a -iX^^D,X a , (8) 

where the metric is rj^ u = diag(— 1, 1, 1, 1), = (— l,cr), and is the gauge covari- 
ant derivative 5 . The Lagrangian for the matter fields £ cn irai contains kinetic term and 
interactions: 

= -D^*D^ - i^a^D^ + £ int , (9) 

where ip is a two component spinor for either left- or right-handed fermions and C m t can 
be obtained from the superpotential W 

Wmssm = uy u QH u - dy A QH d - ey e LH d + fiH u H d . (10) 

using 

Ant = (dw/dfafa)^ + (dw/d&ydw/d&y. (ii) 



5 Here, we have followed the conventions of Ref. 
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The first term in Eq. (fTTj) gives rise to the usual Yukawa coupling [from the first three 
terms in Eq. ffTUl) ]. and the Higgsino mass [from the last term in Eq. ffTUl) ]. The second 
term in Eq. ffTTf) gives rise to all the cubic and quartic scalar interactions. 

The general MSSM superpotential also includes baryon and lepton number violating 
interactions: 

W / al=i = -^\jkLiLj€ k + X' ijk LiQjd k + [ilLiHu, (12) 
W AB =i = ^Kjk^idjdk, (13) 

The simultaneous presence of non-vanishing A' and A" couplings allows for rapid proton 
decay that conflicts with present bounds on the proton lifetime. One way to eliminate 
such terms is to introduce a new symmetry called -R-parity, defined by conservation of the 
quantum number 

P R = (_l)3(B-i)+2 S) ( 14 ) 

where s is the spin of the particle. All SM particles have Pr = +1 while all the super- 
partners have Pr = — 1. If /2-parity is an exact symmetry, then all the terms appearing 
in Eq. ffl2]) and ffl3l) are forbidden and no dangerous proton decay can occur via these 
interactions. 

There are two important phenomenological consequence if -R-parity is exactly conserved: 

• The lightest supersymmetric particle is absolutely stable. 

• SM particles are coupled to even numbers of superpartners (usually two). 

If LSP is colorless and charge neutral, it can be a viable candidate for the cold dark matter, 
imparity conservation also implies that sparticles are produced in pairs in collider experi- 
ments and that each sparticle other than LSP eventually decays into final states containing 
odd numbers of LSPs. Moreover, for low-energy processes involving only SM particles in 
the initial and final states - such as those of interest in this article - supersymmetric contri- 
butions appear only at loop- level (e.g., virtual superpartners are pair produced). However, 
one may relax the constrains of i?-parity conservation while preserving proton stability via, 
e.g., forbidding baryon number violating terms in Eq. ffl3]) . In this case, the LSP is no 
longer stable and tree level SUSY contributions to low energy processes appear through 
imparity violating interactions. In what follows, we will consider the implications of both 
imparity conserving and i?-parity violating (RPV) supersymmetry. 

C. Soft SUSY Breaking 

If supersymmetry is exact, superpartners have the same mass as the corresponding SM 
particles. However, supersymmetry must be broken in nature because superpartners have 
not been experimentally observed at energies where they could be pair produced if they 
are degenerate with SM particles. In order to retain the exact cancellation of quadratic 
Auv dependence of the Higgs mass corrections, all the SUSY breaking couplings must be 
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"soft" (of positive mass dimension). After adding the fermion and scalar contributions of 
Eqs. f fTPj) . the remaining logarithamic correction to the Higgs mass is proportional to the 
soft SUSY breaking masses 6 : 

Am 2 h = -^[5m 2 ln(A 2 y /m 2 ) + ...] , (15) 

where we have taken m| = m 2 + 5m 2 s . Therefore, the soft SUSY breaking mass parame- 
ters (e.g. , 5m 2 s ) are should be below a few TeV to avoid reintroduction of the naturalness 
problem. Throughout this work, we will refer to this scale of SUSY-breaking mass param- 
eters as fa. A brief description of soft SUSY breaking parameters, SUSY particle mass 
spectra and interactions is given below. For a more detailed review of MSSM and related 



phenomenology, see Refs. |4|, [10 



In MSSM, the Lagrangian for soft SUSY breaking terms are 

Aoft = -~(M 3 gg + M 2 WW + M 1 BB) + c.c. 

-(u& u QH u - d& d QH d - e& e LH d ) + c.c. 

— Q^uiqQ — Dm$L — um^u* — Jmjj<? — era?e' — m 2 Hu H^H u — m 2 Hd H* L H d 
-(bH u H d + c.c.) (16) 

The first line gives the gaugino mass M i: i = 1, 2, 3 for U(l)y, SU(2) L and SU(3)c gauginos, 
respectively, and where the boldfaced quantities indicate matrices in flavor space. The 
second line gives the trilinear "A-term" that couples Higgs scalars with left- and right- 
squarks and sleptons. The third line gives the scalar mass m~ , m~ , and m\, for 

squarks, sleptons and Higgs scalars, respectively. Finally, the last line is the bilinear 
6-term, which couples up- and down-type Higgses. In principle, one may also include 
RPV soft interactions that correspond to the terms in the superpotentials H / al=i and 
Wab=i- However, pure scalar RPV interactions are generally not relevant to the low- 
energy observables discussed here, so we will not include them. 

The trilinear A-terms and the soft SUSY breaking squark and slepton masses are in 
general non-diagonal in the flavor basis, a feature that introduces flavor-changing-neutral- 
current (FCNC) effects beyond those that are GIM-suppressed in the SM. Moreover, after 
performing an appropriate set of field redefinitions, £ so ft - together with the /i-term in 
the superpotential - includes 40 CP-violating phases beyond those of the SM (for a use- 



ful discussion, see, e.g., Ref. [16j). In contrast to the effects of the CP-violating phase 
in the Cabibbo-Kobayashi-Maskawa (CKM) matrix, the effects of these new phases are 
not suppressed by the Jarlskog invariant (l7| and light quark Yukawa couplings. Thus, the 
interactions in £ so f t can lead to unsuppressed FCNC and CP-violating effects at low en- 
ergy. On the other hand, both FCNC and CP violation have been tightly constrained by 



6 There are additional logarithmic contributions proportional to the square of the triscalar coupling, a / , 
defined belowl' 
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experiment. Attemps to reconcile these two phenomenological implications of £ so ft with 
experimental bounds on FCNCs and CP-violation are known as the "SUSY flavor" and 
"SUSY CP" problem, respectively. 

A detailed discussion of the SUSY flavor and CP problems appears in Section IVI1 
However, for purposes of illustration, we consider one approach to the flavor problem in 
which it is assumed that and rriq, m?, m|, and m| are diagonal in flavor basis and 
that a u , ad, and a e are proportional to the corresponding Yukawa matrices, y u ,d,e- The 
conventional parameterization thus gives - after diagonalization of the Yukawa matrices - 



A u\ Vc , a d = A d \ y s , a e = A e \ y \ . (17) 






Specific SUSY breaking scenarios have been studied in the literature, which could solve 
the FCNC and CP-violation problems introduced by the soft SUSY breaking terms. If 
SUSY breaking is mediated from an unseen, high energy "hidden sector" to the visible 
weak scale sector (superparticles in MSSM) via gravity, it is called Gravity Mediated SUSY 
Breaking (SUGRA) [l9(. If SUSY breaking is mediated from the hidden sector to the 
visible sector via gauge interactions, it is called Gauge Mediated SUSY Breaking (GMSB) 
20|- Recently, Anomaly Mediated SUSY breaking (AMSB) |2l| and gaugino Mediated 



SUSY Breaking scenarios [22J have also been considered. These models relate the large 
number of parameters in £ so ft to a few parameters associated with SUSY-breaking physics 
at high scales. We will occasionally refer to these model-dependent relations throughout 
this review. 



D. Superparticle Spectrum 

The superpartner mass spectrum emerges after diagonalization of the relevant mass 
matrices. For squarks and sleptons, the mass matrices contain three components. In the 
(Jl, /r) basis, there are mass matrices for the LH fermion superpartners, M RR 
for the fn] and matrices Mj R that mix the two. The LR mixing matrices arise only 
after electroweak symmetry breaking, and to the extent that the triscalar couplings af are 
proportional to the Yukawa couplings as in Eq. (1171) . one expects the effects of this mixing 
to be relatively small except for the third generation sfermions. In Section IIVI we discuss 
low-energy tests of this expectation. In flavor space, each of these matrices is 6 x 6 (3 x 3 
in the case of the sneutrino). Since an extensive discussion of the flavor problem appears in 
Section IVIl we will assume momentarily that the M^ B (A, B = L,R) are flavor diagonal 
for purposes of illustration. In general, one has 

M£ L = + mj + A f (18) 

M^ R = m? + mJ + A f (19) 

with 

A f = ( l( - Q f sin 2 6 W ) cos 2/5M| (20) 
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and 



A f = Q f sin 2 9 W cos 2(3M 2 Z (21) 

M 2 = M 2 = { V ^ f Sin ^ ~ /iYf COS ^ ' " ~ type sfermion (22) 
f [a f cos /3 — yuyf sin /3] , d — type sfermion 



Here m 2 is the mass matrix for the corresponding fermion specie, // and Q / are the third 
component of isospin, and the charge of the fermion, respectively, and tan f3 is the ratio of 
the neutral Higgs vevs tan/5 = (H®)/(H®). 

The diagonal elements depend on the unknown soft SUSY breaking parameters rriq, 
m 2 , etc. while the off-diagonal elements depend on the supersymmetric parameter /x, the 
soft-triscalar coupling af , v and tan j3. Assuming no flavor mixing among different sfermion 
generations, the sfermion mass matrix reduces to a set of 2 x 2 matrices for each flavor. 
The corresponding mass eigenstates Fi t 2 are mixtures of the Jl,r, with the mixing angle 
0; 
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F, 





(23) 

In the more general situation where one allows for flavor mixing, the 6x6 diagonal 
mass matrix is given by 

(M 2 ) diag = Z f tM f 2 Z f (24) 

where 

for each species of sfermion. Hence, a given sfermion mass eigenstate Fj is given in terms 
of the flavor eigenstates fj as 7 

fj = Z 1 / Fj (26) 

where I = 1, 2, 3 indicate the left-handed (LH) flavor states 8 fij and / = 4, 5, 6 refer to the 
right-handed (RH) flavor states /r x _ 3 - Hence, the simultaneous presence of non- vanishing 
Zy and Z 2 ^ would indicate flavor mixing among first and second generation LH sfermions, 
while having both ZJ ^ and Z^ ^ would imply mixing among the LH and RH 
sfermions of the first generation. Note that for sneutrinos, the indices I,j run over only 
1,2,3 and we have no left-right mixing in this case. 



7 Here, we follow the notation and conventions of Ref . 18( . 

8 For scalars, the handness simply indicates that they are the superpartners of the left-handed or the 
right-handed fermions. 
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The gauginos and higgsinos mix with each other since both are charged under the 
electroweak gauge group. The mass matrix for the neutral states ip° = (B, W°, H%, H®) is 



M 



N 



( Mi -c p s w M z s p s w M z \ 

M 2 cpc w M z -S(3C W M Z 

-cpSwMz cpCwMz -fx 

\ sps w M z -spewMz -fi J 

Here we have introduced the abbreviation sp = sin/5, c@ = cos/3, sw 
cw = cos 8w- The mass matrix can be diagonalized by a 4 x 4 unitary matrix N: 

rVLo diag 



(27) 



sin^vi/, and 



N'MjjN -1 , 



(2f 



The mass eigenstates are called neutralinos Xi = Nijip®, i = 1 . . .4, with m x o < m x o < 
m x o < m x o. In the limit that Mz <C M 1; M 2 , \fi\, each neutralino is a pure gaugino or 
Higgsino state, while in general, the neutralino is a mixture of gauginos and Higgsinos. 
Similarly, the mass matrix for charged gaugino and Higgsino ij}^ = (W + , H£, W~ , Hj) 

is 

M 2 V2s M w 



(x 


X T \ 






o ) 





Mi 



The mass eigenstates are called charginos xt 



(29) 



\Z2cpMw yU 

= 1,2 (m± < m ±), which are related to 

Al A2 



the gauge eigenstates by two unitary 2x2 matrices U and V that diagonalize the chargino 
mass matrix: 



Al 

xt 



V 



w~ 



Xi 

X2 



u 



U*XV~ 











(30) 



The SU(3)c* gluino is a color octet fermion, and does not mix with other particles in MSSM. 
Its mass is parametrized by M 3 as defined in Eq. f|T6|) . 

In gravity-mediated and gauge-mediated SUSY breaking models, there is a unification 
relation for SU(3)c, SU(2) L and U(l)y gaugino masses M 3i2) i: 



M, 



M 1 



^2 

j 

a s a 2 Q-i 

where a s , a 2 and a x are related to the couplings of SU(3)c, SU(2) L and U(l)y via 

a (?9 2 5 a 5 gy 2 

a i = - 



47T 



«2 



sin 2 fl 



47T 



3 cos 2 6*vk 3 An 



(31) 



(32) 



This relation holds at any energy scale to one-loop order 9 . In particular, at electroweak 

Lgino 

(33) 



scale, if we take a s = 0.118, a = 1/128 and sin 2 Q w = 0.23, the ratio betweem gaugino 



masses are 



M 3 : M 2 : M x « 7 : 2 : 1. 



9 We do not discuss possible threshold effects at the GUT or Plank scale. 
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However, such a unification relation need not hold in the most general MSSM, and, M 3 , 
M 2 and Mi could be completely independent of each other. 

The MSSM has two complex Higgs doublets, H u and H d , which give mass to up and 
down type fermions, respectively. The potential for the neutral Higgs fields is 

V = (M 2 +mk)I^T + (H 2 +<)l^ (34) 

As in the Standard Model, the minimum of the potential corresponds to a non-zero VEV 
for the neutral Higgs fields, thereby breaking electroweak symmetry. Let us write (H®) = 
v u j y/2, (H d ) = v d j \J2. The sum of v 2 and v d is related to the Z boson mass and the gauge 
couplings as v 2 u + vj = v 2 = AM^/(g 2 + g' 2 ) « (246 GeV) 2 . It is convenient to write the 
ratio of v u and v d as tan (3: tan (3 = v u /v d , where v u = v sin (3 and v d = v cos (3. 

The two complex Higgs doublets contain eight real scalar degrees of freedom. After 
EWSB, two charged and one neutral degree of freedom are the would-be Nambu-Goldstone 
bosons, G ± and G°, that are eaten by W ± and Z to become their longitudinal modes. We 
are, thus, left with five physical Higgs bosons: two neutral CP-even Higgses, h° and H°; 
one neutral CP-odd Higgs, A ; and a pair of charged Higgses, H ± . When m^o 3> Mw, 
h° is the SM-like Higgs. The tree level Higgs masses can be obtained via expanding the 
potential around the Higgs VEVs and diagonalizing the 2x2 mass matrices. One finds 

m 2 A0 = 26/ sin 2/5, (35) 
m 2 H ± = m 2 A o + M^, (36) 

m 2 h o <HO = \ (m 2 A0 +M 2 zT y/ (m\ + M 2 f - AM 2 m 2 A0 cos 2 2/3J . (37) 

The mass eigenstates can be written in terms of the gauge eigenstates as 

A = v^cos/? Im[#°] + sin/3 lm[H$\), (38) 
H + = cos $ H+ + sin/5 Hj*, H = cos/5 H+* + sin j3 H~[, (39) 

H°)~ V Vsina cosa J \ Re[H° d ' - ' ' 



where 



v d 



sin 2a m 2 A0 + M| cos 2a m 2 A0 - Ml 



sin 2(3 m 2 H0 — m? ' cos 2(3 m 2 H0 — mj 



(41) 



The tree level Higgs masses are determined by only two parameters: b and tan (3 (or m^o 
and tan/5). For large m^o, A , H° and H ± are heavy and decouple from low-energy 
observables. On the other hand, the tree level mass of the light CP-even Higgs h° is 
bounded from above: m^o < Mz, which has already been excluded by the current LEP 
Higgs searches [23[. However, the mass of the h° receives large radiative corrections from 
third generation quarks and their superpartners due to the large top Yukawa couplings, 
allowing for a mass large enough to be consistent with present direct search limits. The 
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dominant contribution is from the stop loop: 



AK^1,V-Vln(^) • (42) 

For stop masses around 1 TeV, this correction can push the m h o above the current exper- 
imental bound. Detailed two loop calculations for the light CP-even Higgs mass indicate 
that the upper bound is about 135 GeV. Uncertainties of a few GeV arise from neglected 
higher order effects and the experimental error in the top quark mass j24] ]. 

The couplings of the Higgses to the (s) quarks and (s)leptons are proportional to the 
Yukawa couplings and are, therefore, non-negligible only for the third generation. For the 
low energy precision measurements where only the light quark and leptons are involved, the 
contribution from the Higgs sector can almost always be neglected. However, the details 
of the Higgs sector do affect the phenomenology of SUSY CP-violation, so we will provide 
a brief summary of the status of the experimental searches of MSSM neutral Higgs bosons 
in Sec lVHTl 



E. SUSY Interactions 

The gauge interactions in MSSM can be obtained from the usual SM gauge interactions 
by replacing two of the SM particles with their superpartners. For example, the coupling 
of the SU(2)l gauge boson to quarks is (before EWSB and diagonalization of the quark 
mass matrices) 

- gQ^l ■ W,Q (43) 
while the corresponding gauge- squark-squark interaction is 

igd^-.W^Q + c.c. (44) 
Similarly, supersymmetry leads to a squark-quark-SU(2)i gaugino interaction 

- V2g(Q^Q) ■ W + c.c. (45) 

The corresponding Feynman rules are illustrated in the diagrams of Fig. [TJ _ Additional 
gauge boson-gauge boson-squark-squark interactions appear via the D fl Q*D fl Q term. 

The other fermion-fermion-gauge boson, sfermion-sfermion-gauge boson, and fermion- 
sfermion-gaugino interactions follow a similar pattern. Similarly, the supersymmetrized 
gauge boson self-interation leads to gauge boson-gaugino-gaugino coupling, as shown in 
Fig.d 

The Higgs Yukawa interactions are obtained from the superpotential Wmssm in Eq. ffTUl) . 
It also give rise to Higgsino-quark-squark interactions via the first term in Eq. ffTTj) and 
Higgs-squark-squark interactions via the second term in Eq. ffTTl) . The soft SUSY break- 
ing trilinear A-term give rise to additional Higgs-squark-squark couplings. The coupling 
between the Higgs and the lepton sector can be obtained similarly. 
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FIG. 1: Feynman diagrams for supersymmetric (a) gauge-quark-quark, (b) gauge-squark-squark 
and (c) gaugino-quark-squark vertices. 




FIG. 2: Feynman diagrams for supersymmetric (a) tri-gauge boson coupling, (b) gauge boson- 
gaugino-gaugino coupling. 



Including the effects of EWSB leads to modifications of these expressions, due to (a) 
diagonalization of the quark mass matrices, leading to the presence of the CKM matrix 
in Eqs. (I43 f 45 p ; left-right mixing (as well as possible flavor mixing) among sfermions, 
leading to the presence of mixing matrices Z 1 ^ in Eqs. (I44f45p : and mixing of gauginos and 
Higgsinos into the charginos and neutralinos, leading to factors of the matrices iVy, V!y, 
and Uij etc. in Eq. (1451) . The Feynman rules for these interactions appear several places in 
the literature and we do not reproduce a complete list here. Throughout this article, we 
generally follow the conventions given in Refs. 0,[25|. 



F. i?-parity Violating Interactions 

Additional B- and L- violating interactions may appear in the MSSM if i?-parity viola- 
tion is allowed. Rapid proton decay can still be avoided if we only turn on B or L violating 
terms, but not both simultaneously. The RPV terms in the Lagrangian can be obtained 
from the superpotentials Eqs. (IT2|) and (I13p via Eq. fill I) . For low energy process where 
light quarks are present in the initial and final states, the RPV terms that are of interests 
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are Yukawa-type interactions: 

Crpv, AL=i=\jk{^LiLje\ + LiLjc\) + X' ijk (LiQjdl + LiQjd\ + LiQjd\) + /i-Ljif u ;(46) 

Crpv, ^B=i=K 3 Md)d\ + u\d)d{). (47) 

Note that here we follow the notation of Ref. j4| in which fermion fields are denoted by 
two-compoments Weyl spinors. 

fJ'L Ve 



FIG. 3: Tree-level PR-violating contributions to muon decay. 

Such terms will contribute to the low energy SM process via the exchange of heavy 
scalar quarks or scalar leptons. For example, Fig. [3] shows the RPV contribution from the 
purely leptonic term proportional to \i 2 k to the muon decay amplitude that determines 
the Fermi constant, G^. After a Fierz reordering, the resulting four-fermion amplitude has 
the same structure as the tree-level (V — A) x (V — A) amplitude of the SM, but with 
a normalization determined by the ratio of |Ai2fc| 2 to the square of the exchanged slepton 
mass. More generally, for momentum transfer q 2 <C m 2 , the correction to low-energy SM 
amplitudes from RPV interactions can be parametrized in terms of the following quantities: 



I A 



A ijk (f) = r ' ijkl - > (48) 

with a similar definition for the primed quantities and double primed quantities. 

The quantities Ay*, etc. are constrained by other precision measurements and rare decay 
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searches. A summary of the current experimental bounds can be found in [26(. Here, we 
up-date our earlier global analyses of low-energy constraints on RPV obtained from the 
low-energy observables in Table HTT1 For each observable, we indicate the sensitivity to the 
various A^ k (f) along with a reference to the chapter in this article where a more detailed 
discussion appears. 

The results of our fit are particularly sensitive to tests of the unitarity of the first row 
of the CKM matrix, discussed in Section IIVI As we discuss in there, the status of first 
row CKM unitarity is presently unsettled, so we provide a fit for two scenarios: (a) using 
a value of the kaon decay form factor /+(0) obtained from large Nc QCD studies, leading 
to a deviation from CKM unitarity by roughly two standard deviations; (b) a value for 
/+(0) obtained from recent lattice QCD simulations that implies agreement with unitarity. 
The resulting 95% C.L. ranges for the A^ and A'- k under these two scenarios are given 
in Table 13 
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Quantity 


A' llfe (4) 


A' lfel (<z£) 


Al2fc(4) 


A' 21fc (4) A 


2M(4) 


Value 


Discussion 


s\v ud \y\v ud \ 2 


2 





-2 








-0.0032 ± 0.0014(a) 
-0.0002 ± 0.0015(b) 


Sec. ED 
Sec. ED 


5Q$/Q$ 


-4.82 


5.41 


0.05 








-0.0040 ± 0.0066 


Sec. El 


$R e /fl 


2 








-2 





-0.0042 ± 0.0033 


Sec. ED 


SG^/G^ 








1 








0.00025 ± 0.001875 


Sec. [V] 


SQ e w/Q e w 








-29.8 








0.14 ±0.11 


Sec. |V] 


5R U 








-0.21 


0.22 


0.08 


-0.0033 ± 0.0007 


Sec. |V] 


5Rp 








-0.077 


0.132 


0.32 


-0.0019 ±0.0016 


Sec. El 



TABLE III: P^-violating contributions to 5\V ud \ 2 /\V ud \ 2 , 5Q$/Q$, 5R e/fl , 5G^/G^, 5Q P W /Q P W , 
5Qw/Q\¥i Rv an d R&- Here <5|V^d| 2 /|V^d| 2 means the possible correction to the value of \V ud \ 2 ex- 
tracted from beta-decay that is allowed by first row CKM unitarity tests. See text for description 
of scenario (a) and (b) for <5| I4d| 2 /I^ud| 2 - Columns give the coefficients of the various corrections 
from A£. fc and Ai2k to the different quantities. Next to last column gives the value of the cor- 
responding quantity extracted from experiment assuming only Standard Model contributions to 
the relevant process. Final column gives section of this review containing relevant discussion.. 





(a) large iV c QCD 


(b)lattice QCD 


AW4) 

A' m ($) 
A m (e^) 

A 2 i fc (4) 
A' 2fel (4) 


- 0.0020 
- 0.0017 
0.0013 - 0.0039 
0-0.0015 
- 0.0013 


- 0.0024 
- 0.0019 
0.0006 - 0.0031 
- 0.0014 
- 0.0009 



TABLE IV: 95% C.L. ranges for the A^. obtained from fitting to the low energy observables 
listed in Table [TTTl 

G. SUSY Searches 

Both direct and indirect searches for superparticles have been carried out in various 
experiments [23]. Sparticles can be pair produced at e + e~, pp and pp colliders via the 
intermediate 7*, Z* , gluon or sparticles. Each sparticle subsequently decays into energetic 
lepton, jets plus an LSP for _R-parity conserving scenarios. In most cases, the LSP is a 
neutral weakly interacting particle (e.g., a neutralino) which travels through the detec- 
tor without depositing significant energy. Therefore, typical signatures consist of some 
combination of jets, leptons, possibly photons, and large missing energy. 

The large electron-position (LEP) collider at CERN has completed its running at center 
of mass energy more than twice the mass of Z boson with a few hundred pb _1 integrated 
luminosity. No events inconsistent with the SM have been reported and all visible sparticles 
with mass up to half the Z mass have been excluded. Data taken at center of mass energy 
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above Mz can set stronger limits on the neutralino, chargino, squark and slepton masses, 
although limits would depend on the interplay between the sparticle masses, cross sections 
and decay branching ratios. Charginos are excluded up to 103 GeV except in the case of 
low acceptance or low cross section. The limits on the slepton masses are based on the 
pair production of 11 with / — > lx\, which gives a lower mass bound of about 80—100 GeV. 
Limits on the stop and sbottom masses are about 90 GeV, which varies with the left-right 
squark mixing angle. 

The Tevatron Run I has accumulated about 110 pb _1 of data forpp collisions at ^fs = 1.8 
TeV. Pairs of squarks and gluinos can be produced because of the large strong interaction 
cross sections. Signals of energetic multijets plus missing transverse energy have been 
searched for, with null results being reported. The exclusion region in (rriq, rrig) has been 
derived, and mass larger than 300 GeV has been excluded if m g = rrig. The gluino mass 
has been excluded up to 195 GeV for any squark mass. Charginos and neutralinos can be 
produced via qq' — > xfxt- Leptonic decay of both X\ an d X2 leads to trilepton signals 
which reduces the background significantly. The same sign dilepton signal is possible for 
charginos produced in the squark and gluino cascade decay. Bounds on charginos and 
neutralinos, however, are not as strong as the LEP results. The upgraded Tevatron Run II 
with a/s = 2 TeV and designed integrated luminosity L = 2 fb _1 could cover large regions 
of SUSY parameter space. The Large Hadron Collider (LHC) at CERN of pp collision with 
Ja = 14 TeV and L = 100 fb~ x would cover the region of squarks and gluinos with mass 
up to a few TeV. 

Precision measurements of Z-pole observables, My/, m t and several low-energy observ- 
ables could constrain the MSSM mass spectrum. The bounds from the current precision 
Z-pole measurements on the SUSY parameters are discussed in Sec. IVHI In addition, 
FCNC and CP violation experiments impose strong limits on the possible flavor structure 
and CP-violating phases of MSSM, which heavily constrain the structure of soft-SUSY 
breaking parameters (see Section IVI]) . In the i?-parity conserving MSSM, sparticles could 
affect the low energy measurements via radiative corrections. In general, loop-induced 
SUSY effects in low-energy observables are proportional to a/ir(M/m) 2 , where M is the 
relevant SM particle mass. For m ~ Mw and M ~ Mw, probing these effects, there- 
fore, requires precision of better than one percent. The analsyis and implications of these 
precision measurements constitutes the bulk of the remainder of this review. 

III. RENORMALIZATION 

In this review, we focus on possible manifestations of SUSY in low energy processes 
that may complement the information obtained from high energy collider searches. The 
low energy processes of interest here further break down into two broad classes: (a) those 
allowed by the Standard Model and for which there exist precise SM predictions, and (b) 
processes that are either forbidden by SM symmetries or that are highly suppressed. In 
order to discern the effects of SUSY in the first category of observables, one generally 
requires knowledge of SM predictions at the level of one-loop (or higher) radiative correc- 
tions. In the case of SUSY models that conserve R-parity, SUSY contributions to precision 
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electroweak observables will also appear solely via loop effects. 

In order to make meaningful comparisons between experimental results and predictions 
based on supersymmetric extensions of the SM , one must compute and renormalize super- 
partner loop effects using the same framework as adopted for the SM predictions. Doing 
so requires care, however, as one must employ a regulator that preserves supersymmetry. 
Here, we use dimensional reduction (DR), wherein one works in d = 4 — 2e spacetime 
dimensions while retaining the Clifford algebra appropriate to fermion field operators in 
d = 4 dimensions. Renormalized quantities are obtained by introducing counterterms that 
remove the factors of 1/e — 7 + ln47r that arise in divergent one-loop graphs - a subtrac- 
tion scheme known as modified dimensional reduction, or DR. This scheme represents 
a variation on the more familiar modified minimal subtraction (MS) commonly used for 
the SM. In MS renormalization, one regularizes divergent graphs using dimensional reg- 
ularization, which differs form dimensional reduction by continuing both the number of 
spacetime dimensions as well as the Clifford algebra into d = 4 — 2e dimensions. Doing 
so entails an explicit breaking of supersymmetry, however, as dimensional regularization 
effectively changes the number of fermionic degrees of freedom relative to those of their 
bosonic superpartners in d 7^ 4. In contrast, DR retains the correspondence between the 
number of bosonic and fermionic degrees of freedom. 

Most low energy precision electroweak observables are mediated at lowest order by 
the exchange of a virtual gauge boson (GB), so we consider first the renormalization of 
GB propagators and GB-fermion interactions. Low energy processes generally involve the 
lightest quarks and leptons, and since the corresponding Higgs-fermion interactions are 
suppressed by small Yukawa couplings, we will not discuss renormalization of the Higgs 
sector in detail here 10 . To set our notation, we first discuss renormalization relevant to 
charged current (CC) processes and subsequently discuss the neutral current (NC) sector. 



A. Charged Current Processes 

Radiative corrections to CC amplitudes naturally divide into four topologies: (a) W- 
boson propagator corrections; (b) corrections to the 1^-fermion vertices; (c) fermion prop- 
agator corrections; and (d) box graphs. These corrections have been well-studied in the 
SM and we refer the reader to the extensive literature on the subject for further details 



(see, e.g., Refs. [28|, |29|, |30|, l3l|). As our emphasis in this review lies with the effects of 
supersymmetric particles, we show in Fig. [H illustrative superpartner contributions to each 
type of radiative correction. 

One loop corrections to the IV-boson propagator, fermion propator, and P^-fermion 
vertices are divergent, so one must carry out the appropriate renormalizaton. After such 
renormalization, the W-boson propagator, iD^ u {k) takes the general form in the Feynman 



10 A brief discussion of MSSM neutral Higgs searches is given in Section [VIIH and more extensive review 
can be found in Ref. [Ql. 
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w w 

(b) 




(c) (d) 

FIG. 4: Representative supersymmetric corrections to charged current observables: (a) VK-boson 
propagator corrections; (b) vertex corrections; (c) external leg corrections; and (d) box graph 
contributions. 



gauge 



{k ) + Ln U D ww 



where the transverse and longitudinal projection operators are given by 



T, 



[IV 



■'/IV 



9/j.u "I" k^k u / k 
k^kyjk 



(49) 



(50) 
(51) 



and Dy^(k 2 ) are finite scalar functions and the hat indicates quantities renormalized in 
the DR scheme. In low-energy processes, effects associated with the longitudinal term are 
suppressed by light fermion masses, so we will not discuss the component further. The 
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renormalized transverse component is given by 



k 2 -M 2 w + ti T ww (k 2 ) 



(52) 



where Mw is the finite part of the bare W-boson mass parameter appearing in the renor- 
malized Lagrangian after electroweak symmetry breaking and Yl^ w (k 2 ) gives the finite 
loop contribution after DR subtraction is performed. Both Mw and Il^r W (k 2 ) depend on 
the t'Hooft renormalization scale fi. However, the physical W-boson mass - defined by the 
value of k 2 giving [D^yyijt 2 = M^)]" 1 = - is /i-independent. The finite residue Z w of 
the pole in £>ww * s gi ven by 



Zw — 



i + n 



{M 2 V 



WW\ 1V1 W 



The corresponding expression for the renormalized, inverse fermion propagator is 



Sj\k) =ft-m f + \A L {k 2 )^ + B L {k 2 )] P L + \A R {k 2 )^ + B R (k 2 ) 



R 



(53) 



(54) 



where Pl,r are the left- and right-handed projectors and the Ai tR and B^ R contain the 
finite loop contributions. The physical fermion mass is given by 



m f = 



m f - -B L (m)) - -B R (m) 



1 + \A L {m 2 ) + l -A R {m 2 ) 



(55) 



while the residue of the pole is 



1 + A L (m})P L + A R (m 2 f )P R 



(56) 



Note that for CC interactions in the SM, the left-handed (LH) components given the domi- 
nant contribution to physical amplitudes, as the presence of right-handed (RH) components 
will be suppressed by factors of the fermion masses 11 . 

The renormalized vertex functions for CC amplitudes are relatively straightforward. We 
illustrate using the muon decay process [T — > u^W^ , for which the tree-level amplitude is 



9 - 



V2 



(57) 



After one-loop renormalization, one has 



iM^ c + iM 



cc 

vertex 



1 + F v (k 2 ) v^PlV 



(5c 



11 For example, the weak magnetic moment operator in the SM is chirality-odd and is generated by one- loop 
vertex corrections that contain single insertions of the Yukawa interaction. 
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where is the running SU(2)x gauge coupling and Fy(k 2 ) is the finite part of the 
one-loop vertex correction. 

The vertex and propagator corrections outlined above will contribute to the four-fermion 
amplitudes that describe low energy CC processes of interest to us, such as /i- and /3-decay. 
Additional, but finite, SM one-loop contributions are generated by box graphs involving 
the exchange of two vector bosons. To the extent that the external masses and momenta 
are small compared to the weak scale, the box contributions will have the form of a product 
of two left-handed currents, (V — A) <g> (V — A). In the case of \i~ — > z/ M e~z/ e 



iM 



cc 

box 



2M 2 V 



(59) 



where the + • • • indicate terms whose structure differs from the {V — A) <g> (V—A) structure 
of the tree-level CC amplitude 12 . In the SM, such terms will be suppressed by factors of 
m^/M^. Superpartner loops can lead to relatively unsuppressed non-(V^ — A) <g> (V — 
A) contributions in the presence of mixing between left- and right-handed sfermions (see 
Section HVj). 

Including the box contribution along with the other renormalized one-loop contributions, 

~ 1 /2 

taking into account the factors of 1/Z, that arise in the standard reduction formulae, 
and working in the k 2 << limit, one has 



lMl C ee + lACL + ^propagator + iM 



CC 
box 



2M 2 W 



1 + 



LL ww 



(0) 



M 2 W 



(60) 



+ i£(0) + i^(0) + F£(0) + F£(0) + 8 } 



box 



x z/ M 7 P L \i ejxP L u s + 



or 



iM cc 

J 1 one— loop 



2M 2 W 



1 | n w(0) | l 



iy^ x P L fi e"i\P L Ve + 



(61) 



where 5vb denotes the fermion propagator, vertex, and box graph contributions. 

The resulting rate for muon decay - including the bremstraahlung contribution - is then 
given by 



1 



m, 



fr 



96vr 3 \ 8M^ 



Ivl w 



brem 



(62) 



m, 



192^ [1 + 6QED] 



Here, vg is the v-spinor for the electron antineutrino. 
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where r M is the muon lifetime and the second equality defines the /i-decay Fermi constant, 



G^, and where 



x a f 25 



7T~ 



+ 



(63) 



denotes the contributions from real and virtual photons computed in the Fermi theory of 
the decay. Thus, one has 



G. 



\/2 8M, 



M w 



(1 + Af A 



(64) 



where Syg is given by 5ys but with the Fermi theory QED contributions subtracted out. 

Along with the fine structure constant and the Z-boson mass, Mz, the value of is 
one of the three most precisely determined parameters in the gauge sector of the SM. Thus, 
for purposes of computing other electroweak observables, it is conventional to express g 2 
in terms of G^, Mw, and the correction Ar M : 



9 



> _ 8M W G^ 



V2 1 + Af A 



(65) 



As a particular application, we consider the corresponding amplitude for the /3-decay d 
ut~v e : 



iM 



f3— decay 
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2M W 



V ud (1 + Af p ) *y*P L d e^ x P L v- e 



(66) 



■ G» 
V2 



(1 + Afjs - Affj) m 7 a (1 - 75) d e 7A (l - j B )u B 



where we have omitted terms that are suppressed by quark masses and where Arp includes 
virtual photon contributions, in contrast to Ar M . 



B. Neutral Current Processes 



The renormalization of neutral current (NC) amplitudes follows similar lines, though 
additional features arise due to mixing between the SU(2) j, and U(l)y sectors (for early 



computations in the SM, see, e.g., Refs. [32|, |33j, [34|, [35|). The general structure of the 



renormalized amplitude for the neutral current process £ + f^£ + f is 

*M^ loop = -i^fa c {J?)l^(fc + &y 6 ) I f lx (gf v + gf A j 5 ) f + b0X , (67) 

where £ and / denote the lepton and fermion spinors, respectively, and "+box" denotes 
the box diagram contributions. The quantity /?nc is a normalization factor common to 
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all four-fermion NC processes that can be expressed in terms of gauge boson masses, the 
fll v (k 2 ), and Ary 



P(k 2 ] 



NC 



1 + 



Re Il T zz (M 2 z ) U^ w (0) 



k 2 -M 2 + iM z T z 
K T zz(k 2 ) - n T zz(M z 



k 2 -M 2 



Mf 



M 2 W 



where 



(6* 



(69) 



Ml = Ml-nl^(M|) 

and MzTz = ImIF| z (A; 2 ). For k 2 < Mf, Y Z = 0. Representative superpartner contribu- 
tions to Tlzz are shown in Fig. [51 
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FIG. 5: Representative supersymmetric corrections to neutral current observables: (a) Z-boson 
propagator and Z-j mixing conributions; (b) vertex corrections; (c) external leg corrections; and 
(d) box graph contributions 
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The renormalized vector and axial vector couplings of the Z Q to fermion f - §y and 
g^ A - can be expressed in terms of the weak mixing angle, 9w, and an associated uni- 
versal renormalization factor, along with process-specific vector and axial vector radiative 
corrections: 

g f v = 2l( - 4£(fc 2 ,/i) sin 2 § w (ji)Q f + X f v (70) 
gf A = _ 2 // + . (71) 

where l( and Qf are the fermion isospin and charge, respectively; sin 2 9w{n) = S 2 (/i) 
defines the weak mixing angle in the MS scheme: 

sm 2 § w ^) = ; (72) 

9W + 9 W 

and Ay A are process-dependent corrections that vanish at tree- level. Here g and g' are the 
SU(2)i and U(l)y coupling, respectively. 

As with the running QED and QCD couplings, d(/i) and a s (/x), respectively, the running 
of the weak mixing angle is a prediction of the SM and provides a useful benchmark for 
precision studies in the NC sector [36[ . A renormalization group-improved SM prediction 
for s 2 (/x) in the MS scheme has recently been carried out in Ref. [371 ]. where logarith- 
mic contributions of the form a n ln n (n / /j,q) (with fiQ being a reference scale) have been 
summed to all orders. Additional subleading contributions of the form a n+1 ln n (/x//io) and 
aa" +fe ln n (/x//xo) with k = 0,1,2 were also included in that analysis, and a refined esti- 
mate of the hadronic physics uncertainty associated with light-quark loops at low scales 
performed (see below). The results are shown in Fig. [6J where the scale /i has been chosen 
to be Q = ^/jFf for a process occurring at squared momentum transfer k 2 . The reference 
scale has been chose to be /io = Mz and the running of s 2 (Q) normalized to reproduce 
its value at the Z°-pole : sin 2 ^(M^) = 0.23122(15). The discontinuities in the curve of 
Fig. [6] correspond to particle thresholds, below which a particle of the corresponding mass 
decouples from the running. The change in sign of the slope at Q = M\y arises from the 
difference in sign of the gauge boson and fermion contributions to the f3 function for s 2 (fi). 
Note that threshold matching conditions in the DR-scheme will differ from those in the 
MS framework due to the differences in continuation of the Clifford algebra into d = 4 — Is 
dimensions 38. 39l . 



For purposes of this review, it is particularly interesting to quote the value of the the 
running weak mixing angle at Q = 0[37j]: 

sin 2 § w (0) = 0.23867 ± 0.00016 (73) 

where the error is dominated by the experimental error in siio 2 §w(Mz) and where the 
value of s 2 (/x) at the two scales differs by roughly three percent. In Section [V] we describe 
a variety of low-energy NC experiments designed to test this running. 

By itself, s 2 (^) is not an observable since it depends on the renormalization scale. One 
may, however, define an effective weak mixing angle that is /^-independent and may in 
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FIG. 6: Calculated running of the weak mixing angle in the SM, defined in the MS renormalization 
scheme. Also shown are the experimental results from APV, neutrino DIS (i/-DIS), parity vio- 
lating asymmetry measurement at E158 (Apy), the expected precision of Qweak and the lepton 
forward-backward asymmetry measurement at CDF (App)- This plot is taken from Ref. 401 ] . 



principle be isolated experimentally by comparing experiments with different species of 
fermions: 

sin 2 6 w (k 2 ) cS = k(k 2 , n) sin 2 9 w (n) (74) 

where the quantity k(k 2 ,[i) describes a class of electroweak radiative corrections that is 
independent of the species of fermion involved in the NC interaction. Contributions to 
k(k 2 ,n) arise primarily from the Z-'y mixing tensor: 

nf 7 (fc 2 ) = nl 7 (A; 2 )T^ + n| 7 (A; 2 )L^ . (75) 

Note that in general, the functions n^ 7 (/c 2 ) depend on the choice of electroweak gauge 
parameter, so to arrive at a gauge-independent £(fc 2 ,/i), a prescription for removing the 
gauge-dependent components of IT| 7 (/c 2 ) must be employed [ij]]. 

For processes involving \k 2 \ <C M§, contributions from light fermions to k(k 2 ,fi) can 
lead to the presence of large logarithms when one chooses fi = Mz- The presence of 
these logarithms can spoil the expected behavior of the perturbation series unless they are 
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summed to all orders. To illustrate, consider the amplitude for low-energy, parity-violating 
M0ller scattering: 

M e P e v = ^=p NC m c v9 e A ej^e (76) 

with 

Qw = Pnc(O) g e v 9 e A = Pnc(O) [-1 + 4/c(0, /i)s 2 (/i) +X f v + + As 2 )] +■■■ (77) 

being the "weak charge" of the electron and with the + • ■ • indicating box diagram con- 
tribution and terms of order (a/47r) 2 . At tree-level (k — > 1, Xy A — > 0), the weak charge is 
suppressed, since s 2 is numerically close to 1/4: Q^ tree ~ —0.1. Inclusion of one-loop SM 
radiative corrections reduce the magnitude of Qyy by nearly 40 %, owing largely to the near 
cancellation between the first two terms in Eq. (1771) and the presence of large logarithms in 
k(0, fi) when \x is chosen to be M z as is conventional j36|]. Given these two considerations, 
one would expect the relative size of two-loop corrections to to be considerably larger 
than the nominal a /Arc scale. 

In order to improve the convergence of the SM prediction for Qyy, one would like to 
sum the large logarithms to all orders. The use of the running sin 2 9\y(n) provides a means 
for doing so. By choosing fi ~ Q in both k(k 2 ,fi) and sin 2 9 w (fj,), using the requirement 
that their product is /i-independent as per Eq. ( !74|) . and solving the RG equations for 
sin 2 #jy(/i) as in Ref. [13], one effectively moves all the large logarithms from k{k 2 ,jj) into 
sin 2 6\y(n) and sums them to all orders. The result is a SM prediction for sin 2 9w(k 2 ) eS 
with substantially smaller truncation error than would be obtained by the naive application 
of perturbation theory to one-loop order. 

In the case of superpartner loop contributions to low-energy observables, it is sufficient 
to include their effects solely in the form factor k(k 2 , fi) while choosing fi = Mz (illustrative 
contributions to il^ 7 are shown in Fig. [5]). In addition, one should include their effects in 
the value of sin 2 §w(M z )- One may adopt two different strategies for doing so: 

(1) Include their effects implicitly in the value of sin 2 9\y(Mz) that is obtained from fits 
to precision Z-pole observables. To be consistent, such fits must include the effects 
of superpartner contributions to 0(a) electroweak radiative corrections, and to our 
knowledge, such an extraction has not been carried out using LEP and SLD data 
in a way that does not rely on a model for SUSY-breaking mediation (see, e.g., 



Ref. [42j, |43| and references therein) . 



(2) Use the requirements of electroweak symmetry to compute superpartner contribu- 
tions to siia 2 9 w(M z ) explicitly. Specifically, using 

e 2 ^) = g 2 ^)s 2 ^) (78) 

M 2 W = M 2 z c 2 , (79) 

writing 

d(/i) = a + Aa(fj) (80) 
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where a is the fine structure constant, employing Eqs. (I65|I69I) . and choosing /x = M z 
we obtain 



Af (/ x) = Af M+ -- . (82) 



s 2 (M z )c 2 (M z ) = -= — (81) 

1 Z) V Z) [1 — Af(M z )] 1 ' 

where 

Aa IL T zz (M l /x) 

Thus, by computing the superpartner loop corrections to the various terms in Eq. (182j) 
and employing Eq. (ED), one may determine the predicted shift in s 2 (M z ) explicitly 
for a given set of SUSY parameters. In the remainder of this article, we follow this 
second strategy. 

The corrections contained in the X VA are fermion-specific, containing the Zff vertex 
and external leg contributions. Representative superpartner contributions are shown in 
Fig. [51 In the case of low-energy interactions involving one or more charged fermions, an 
additional contribution to the X v is generated by 7 exchange and involves the so-called 



anapole coupling of the fermion [4J, |45 



^anapole = d v F^ , (83) 

where Fa is the dimensionless anapole moment and M is an appropriate mass scale. The 
interaction £ a napoie generates a contribution to the fermion martrix element of the electro- 
magnetic current: 

(P'l ^? M (0) \p) = U(P) [FiK + (84) 



(fc 2 7 M - 75 + — — a^k u ^ 5 U(p) 



where k = p' — p and where Fx, F2, Fa, and Fe give the Dirac, Pauli, anapole, and electric 
dipole form factors, respectively 13 . Both the anapole and electric dipole couplings to the 
photon are parity-odd, while the electric dipole coupling is also odd under time-reversal. 
Since only weak interactions can give rise to a parity-odd photon-fermion coupling, we 
choose M = M z in Eq. (1531 . 

From Eq. (1531 one sees that the anapole coupling gives rise to a contact interaction 
in co-ordinate space, since dyF^ = with j M being the current of the other fermion 
involved in the low-energy interaction. This feature is illustrated in momentum space in 
Fig. 0, where the contribution of the anapole coupling of fermion /' to the scattering from 
fermion / is shown. The I/tk^ term in Eq. (|84|) gives a vanishing contribution due to current 
conservation, while the k 2 term cancels the 1/k 2 from the photon propagator to yield a 



13 



Note that the overall sign of the anapole term in Eq. ([51)1 differs from the convention used in Ref. 46] 
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four-fermion contact interaction proportional to (k 2 )Qf. Since this interaction involves 
a coupling to the vector current of fermion /, it corresponds to a contribution to g(r. 

Note that at low energies for which \k 2 \ <C M|, the pref actor in /S(/c 2 )nc becomes k 2 - 
independent constant, signaling that the Z-exchange contribution is also a contact inter- 
action. In this regime, one has no experimental, kinematic handle with which to separate 
the anapole and Z-exchange contributions 14 . Indeed, the coupling Fa itself depends on 
the choice of electroweak gauge, and only the complete one-loop scattering amplitude that 
includes all 0{a) electroweak radiative corrections (including Fa) is gauge-independent 
(see Ref. jiBJ and references therein). Nonetheless, when classifying various topologies of 
the one-loop corrections, it is useful to separate out the anapole contributions, and in doing 

so, we note that superpartner loop contributions to Fa are gauge-independent. In what 

f' 

follows, we focus on the low-energy // , in which the F A -contribution to the product of 
vector coupling §y and the axial vector coupling g A is given by 

(gWv) =-lGc 2 s 2 Q f F% . (85) 
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FIG. 7: Anapole contributions to the NC interaction between two fermions. 



In studies of precision Z-pole observables, it has been useful to characterize possible 
corrections to the gauge boson pro pag ators from new heavy particles in terms of the so- 
called oblique parameters, S, T, U 47, 48, 0, 50, 51, [H, 53, 54]: 



In contrast, at k 2 ~ M|, the anapole contribution becomes negligible in contrast to the resonating 
Z-exchange amplitude. 
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where the superscript "New" indicates that only the new physics contributions to the self- 
energies are included. Contributions to gauge-boson self energies can be expressed entirely 
in terms of the oblique parameters S, T, and U in the limit that M NEW S> M z . 

However, since present collider limits allow for fairly light superpartners, we do not work 
in this limit 15 . Consequently, the corrections arising from the photon self-energy (Il 77 ) and 
7-Z mixing tensor (n^ 7 ) contain a residual /c 2 -dependence not embodied by the oblique 
parameters. Expressing the contributions to p and sin 2 9 w (k 2 ) eS = k(k 2 , p) sin 2 9w(p) in 
terms of S,T, and U we obtain: 
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17) 



where k 2 is the typical momentum transfer for a given process. For low energy interactions, 
k 2 — > 0. Note that we have included in 8 sin 2 9$ both the the contribution from Il Zl (k 2 ) / k 2 
that enters k(k 2 , p) as well as the shift in s 2 (Mf ) obtained from Eq. (1ST)) as discussed above. 

In analyzing SUSY radiative corrections to low-energy observables, Eqs. (1ST)) provide 
a useful means of incorporating constraints on new physics from precision Z°-pole ob- 
servables. For example, 8p SVSY is highly constrained by bounds on T obtained from such 
observables. In contrast, [8 sin 2 9w(k 2 ) eS ] SVSY is less stringently constrained. As we discuss 
in Section |V] below, the unconstrained contributions to the effective weak mixing angle can 
lead to relatively large effects in some low-energy NC processes. 



15 It is possible to extend the oblique parameterization in this case with three additional parameters [5^. 
For the low-energy observables of interest here, this extended oblique approximation is not especially 
useful. 
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C. Theoretical Uncertainties in Electroweak Radiative Corrections 

An important consideration in exploiting low-energy, precision electroweak observables 
as a probe of SUSY is to ensure that the theoretical uncertainties associated with SM con- 
tributions are well-below the level of possible SUSY effects. The SM uncertainties generally 
involve one of two considerations: (i) neglect of higher order electroweak contributions, and 
(ii) contributions from strong interactions. While an extensive discussion of these consid- 
erations goes beyond the scope of the present article, we give here a brief overview of the 
strategies employed to address them. 

Nominally, one expects the one-loop contributions to quantities such as Af^, k, etc. to 
be of order a/n ~ 10 -3 , so that neglect of two- and higher-loop effects is well justified 
for the present level of experimental sensitivity. Moreover, since SUSY loop contributions 
must generally decouple in the rh — > oo limit, one expects the relative magnitude of their 
contributions to be 



where M is the relevant SM mass and rh is a generic superpartner mass. For weak processes, 
such as \i- and /3-decay, one has M — > My/-, and to the extent that rh is not too different 
from the weak scale, one would expect <5susy loop to be comparable in magnitude, or slighter 
smaller than, the scale of one-loop, SM electroweak corrections. Thus, one would expect 
neglect of two-loop SM contributions to be a justifiable approximation. As discussed 
above, however, exceptions may occur when the one-loop SM contributions contain large 
logarithms, when the tree-level SM amplitudes are suppressed, or both. In such situations, 
the summing terms of the form a n ln n (////xo) is essential, and the RG equations can be 
employed for this purpose. 

Reduction of theoretical uncertainties associated with QCD corrections is generally more 
challenging. Short-distance QCD contributions can be treated using the operator product 
expansion (OPE), and the resulting correction to a given order in a s achieved with suffi- 
cient effort. In the case of PV electron-proton scattering, for example, the one-loop WW 
box contribution is anomalously - but not logarithmically - enhanced, and its effect on 
the proton weak charge, Q^r, nearly cancels that of the large logarithms appearing in k. 
Since the semileptonic, WW box graphs involve hadronic intermediate states one might 
expect relatively important QCD corrections to the one-loop amplitude. In this case, the 
loop integral is dominated by dominated by high momentum scales (k 2 ~ M^), so the 



for a total QCD correction of ~ —0.7%. 

A more problematic situation arises for one-loop corrections that sample momenta of 
order the hadronic scale. To illustrate, we again consider PV electron scattering. For both 
PV M0ller and elastic ep scattering, light quark loop contributions to 11^ lead to hadronic 
uncertainties in k(0, //). Traditionally, light quark contributions have been computed by 
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relating fl£ to the a(e + e~ — > hadrons) via dispersion relations 571] . much as one does in 
computing hadronic vacuum polarization contributions to the muon anomalous magnetic 
moment. In the case of II^ 7 , however, additional assumptions regarding flavor symmetry in 
the current-current correlator are needed in order to make use of e + e~ data. Recently, these 
assumptions have been examined and more stringent bounds on the hadronic uncertainty 
in £(0,/i) obtained [37 . 

For semileptonic processes, additional hadronic uncertainties appear in box graphs that 
contain one 7 and one weak gauge boson. In contrast to the situation for the WW^-box 
graphs, the 7Z loop integral samples all mometum from the hadronic scale to the weak 
scale. Neglecting the short-distance, perturbative QCD corrections, one finds 
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(90) 



where A is a scale characterizing the transition between the perturbative and non- 
perturbative domains and C 7 ^(A) parameterizes contributions to the loop integral from 
momenta ^|/c 2 | < A. The coefficient of logarithm in Eq. (|90|) is determined by short 
distance dynamics and can be calculated reliably in perturbation theory. However, the 
values of both A and C 7 z(A) are sensitive to long-distance scales and have not, as yet, 
been computed from first principles in QCD. A similar contribution arises in neutron, nu- 
clear, and pion /3-decay. An estimate of the theoretical uncertainty associated with these 
contributions had been made by varying A over the range 400 < A < 1600 MeV. Re- 
cently, Marciano and Sirlin observed that for the box, both the pQCD corrections to 
the logarithmic term as well as the value of A could be obtained by comparison with the 
theoretical expression for the Bjorken Sum Rule using isospin symmetry [3 lfl. As a result, 
these authors have reduced the previously-quoted theoretical error by a factor of two. The 
analogous treatment of the 7Z box is more complex, since one cannot obtain the isoscalar 
contribution from isospin arguments. In both cases, the more refined estimates of the un- 
certainty associated with the low-energy constants C 1 z and C 7 w remain to be performed. 
Fortunately, the impact of the uncertainty in SQ^dZ box) due to C^z is suppressed by 
the overall factor of 1 — 4s 2 ~ 0.1. 

When discussing the implications of various low-energy observables for SUSY, we will 
also summarize the current situation regarding hadronic uncertainties in the SM predic- 
tions. 



IV. CHARGED CURRENT PROCESSES 



Historically, the study of low energy charged current (CC) processes have played an 
important role in developing the SM, in determining its parameters, and in testing its 
self-consistency at the level of one-loop radiative corrections. Indeed, the observation of a 
parity- violating asymmetry in the /3-decay of polarized 60 Co [58[ and /i + -decay[59[ confirmed 
Lee and Yang's hypothesis of parity violation in the weak interaction [60] and pointed the 
way toward the V — A structure of the weak charged currents. Measurements of the 
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muon lifetime yield the parameter that is one of the three independent, experimental 
inputs needed for the gauge sector of the theory. Studies of nuclear /3-decay give the 
most precisely-known element of the CKM matrix - V u d, while measurements of branching 
ratios for kaon leptonic decays yield a precise value for V us that allows for stringent tests of 
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the unitarity property of the CKM matrix 16 (for recent discussions, see Refs. [6ll . 
Comparisons of the widths r[7r — > fiv^)] and T(7r — > ez/ e (7)] have provided tests of the 



universality of CC leptonic interactions at the few parts per thousand level [64j. M- The 
theoretical interpretation of these precise measurements in terms of the SM has required 
comprehensive calculations of one-loop radiative corrections to the tree-level amplitudes. 
The vast majority of this work has been carried out by Sirlin and Marciano, dating back to 
the classic treatment within the current algebra framework by Sirlin 28| . The implications 
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for various SM extensions have been analyzed extensively by Herczeg and others [6( 

Interest in precise studies of low energy processes remains high, as reviewed recently in 
Ref. 0]. While an extensive survey of the field can be found in that work, we highlight 
recent developments that motivate the study of CC processes from the standpoint of SUSY. 
Experimentally one has seen: 

i) New measurements of the Michel parameters that characterize the spectral shape, 
angular distribution, and polarization properties in polarized //-decay (68|, 6^, 7(j (for 
a recent global analysis, see Ref. fn\ ) 



ii) New efforts to measure r M with an order-of-magnitude improvement in precision [7 

i3 



iii) Recent Penning trap measurements of "superallowed" nuclear /5-decay Q-values |74j . 



751 ] with significant implications for tests of first-row CKM unitarity 



iv) New measurements of the neutron lifetime, r- [76| . and decay correlation coefficients 
that determine V u d in a manner free from possible nuclear structure ambiguities 771 . 



v) Extensive new measurements of kaon leptonic decay branching ratios that could imply 
significant changes in the decades-old value of V us (recently reviewed in Ref. 63]; also, 
see below) 



vi) Improved precision in the pion /3-decay branching ratio 81 



vii) New efforts to measure the ratio R e /n = T[ir — > ei/ e (7)]/r[7r — > ^1^(7)] [82], [83 



The theoretical interpretation of the semileptonic decays has also been sharpened 
through 



16 The value of V u b is also required, but its magnitude is too small to be relevant. 
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A new analysis of strong interaction uncertainties in AfY that associated with Wj 
box graphs that reduces the theoretical uncertainty in the extraction of V U( i from 
/3-decay rates by a factor of two[3~D] 

Computations of the 0(p 6 ) loop corrections to the kaon decay form factor f+(t) 
whose value at t = is needed in order to extract V us from kaon decay branching 



ratios 1 84, 85 



iii) New analyses of the 0{p ) counterterm contributions to /+ (0) using large Nc 
QCD[86| and lattice QCD computations (see below) 



Given the level of activity in this area, consideration of the implications for SUSY is a 
timely activity. In reviewing these implications, we also provide additional background on 
the theoretical and experimental issues for the benefit of readers who may not be familiar 
with the field. We begin with the purely leptonic CC interaction that gives rise to muon 
decay and follow with an extensive discussion of low-energy semileptonic CC processes. We 
divide the latter discussion into several parts: (1) general considerations for semileptonic 
CC processes; (2) pion leptonic decays and the related implications for SUSY; (3) neutron 
and nuclear /3-decay; (4) pion and (5) kaon /3-decay; (5) implications of first row CKM 
unitarity tests for SUSY. 



A. Muon Decay 



As discussed in Section lITTj the measurements of the muon lifetime generally do not, by 
themselves, provide information on non-SM physics. Rather, the value of G„ as extracted 



from using Eq. (I<32"j) provides a key input into SM predictions for other observables, and 
the presence or absence of deviations from these predictions provides information about 
various SM extensions. For example, one may use of Eq. (l8"Tj) for this purpose, treating 
a, G/j,, Mz, and sin 2 9 w(Mz) as independent, experimentally determined quantities and 
computing the correction Af(Mz) in the SM[87|]. The degree of self-consistency among 
these quantities in Eq. flHTj) will lead to constraints on any deviation of Af(M^) from its 
SM value associated with possible new physics, such as SUSY radiative corrections to the 
/i-decay amplitude or the Z-boson self-energy [see Eq. f[8"2"j) ]. 

Studies of the muon spectral shape and angular distribution, along with the electron 
polarization, can provide additional handles on non-SM contributions. The spectrum and 
polarization are typically described by the eleven Michel parameters 88|, [89| . Four of them 
(p, r], 5, £) characterize the spectral shape and angular distribution: 
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where x = |p e |/|Pe|max! = cos _1 (]3 e • s M ), P^ is the /x ± polarization, and h(x) and 5 (a;) are 
momentum dependent radiative corrections. Five additional parameters (£', 77", a/A, 
ft/ A) are needed to describe the electron transverse and longitudinal polarization and two 
more (at' /A, ft' /A) parameterize the T-odd correlation of the final state lepton spin and 
momenta with the muon polarization. The parameter 77 also characterizes deviations of 
the muon lifetime from its value in the pure V — A Fermi theory, as can be seen from the 
corresponding modification of the total decay rate: 
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In the Standard Model, one has p = 5 = 3/4, P M £ = 1, and 77 = 0, so that from 
measurements of r M one obtains a fractional uncertainty in the Fermi constant of AG M / = 
9 x 10~ 6 . Allowing for 77 7^ due to possible non-SM contributions to the decay amplitude 
and constraining such effects with experimental determinations of the Michel parameters 
can result in a larger uncertainty in G^. A recent measurement of transverse positron 
polarization from p + decay yields rj = (71 ± 37 ± 5) x 10 -3 , leading to an increase in the 
relative error on G M by a factor of 40 [70]. The results of this experiment also lead to new 
values for the parameters rj", a' and ft'. 

New measurements of p and 5 have also been completed, and a new global analysis of 



the Michel parameters has been carried out in Ref. [7l|. It is conventional to analyze the 



results in terms of the effective, four-fermion Lagrangian 



£J i ~ decay 



4GV 
V2 



(93) 



where the sum runs over Dirac matrices T 7 = 1 (S), 7" (V), and a al3 (T) and the subscripts 
e and p denote the chirality (R,L) of the final state lepton and muon, respectively 17 . The 
SM has g\ L = 1 with all other couplings vanishing, leading to the SM values for p, 5, P M £, 
and 77 noted above. The sensitivity of the Michel parameters to non-SM interactions is 
obtained by expressing them in terms of the general set of couplings g]^. For example, one 
has 



\9 V LR \ 2 + \9 V rl\ 2 + 2 



\9l R \ 2 + \g T RL \ 2 ) + ^ (g s LR g T L * R + g s RL g T R l) 



(94) 



In order for SUSY to affect the Michel spectrum or lepton polarization in a discernible 
way, it must generate contributions to the effective Lagrangian (193!) other than those asso- 
ciated with the g\ L term. If one assumes conservation of R-parity, then such contributions 
can be generated at the one-loop level via the box graphs of Fig. [8] (a) . The corresponding 



17 The use of the subscript "//" to denote both the chirality of the muon and the flavor of the corresponding 
neutrino is an unfortunate historical convention. 
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amplitudes contribute to both g\ L and g RR . One has [90 



9RR, loop 



Z l J*Z 2 JZfZf* (N* 2 - tanOwN^) N n (N f2 - tan6 w N fl ) NJ n 



(95) 



_ z iJ* Z 2j Z u Z 4i* ( N * 2 _ tan6wN * i ) Njl ( Nf2 _ tau^vx) ^1 

where the Z^ , Zfj , Uij, and iVy are the sneutrino, slepton, chargino, and neutralino mixing 
matrices, respectively, defined in Section HT1 and 



T n (a, 6, c, d) 



i 

JQ 



1-x 



dy 



l-x-y 



dz [ax + by + cz + d(l — x — y — z)\ n . (96) 



The term in Eq. f[9"3"l) involves scalar couplings between the right-handed (RH) charged 
leptons and left-handed (LH) neutrinos. Although CC interactions are associated with 
the SU(2)i sector of the MSSM, couplings to the RH charged leptons can arise via either 
lepton Yukawa interactions or mixing of the LH and RH slepton weak eigenstates into mass 
eigenstates. In obtaining Eq. (1951) we have retained only contributions associated with the 
latter, as reflected in the presence of the matrices Z^ l , I = 1, 2 in each of the terms. 




V 

I 
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FIG. 8: Box graphs contributing to (a) the muon decay parameter g RR and (b) (3 decay parameters 
a RRi a RL' an d al RL- 



As discussed in Ref. [90|, the products Z^*Z^ and Zf^Z^* lead to lepton flavor changing 
amplitudes at one-loop order and are, thus, highly constrained by searches for lepton flavor 
violating processes such as \i — > ej. Thus, for practical purposes, one may neglect the last 
two terms in Eq. fl95|) . In contrast, the first term in g RRloop is flavor diagonal but requires 
the presence of left-right mixing among first and second generation sleptons. The product 
r?2i* Z 5i a | go en t ers the SUSY contribution to a^ = (g fl —2)/2, the muon anomalous magnetic 
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moment. Since the magnetic moment operator is chirality odd, a non- vanishing SUSY 
contribution requires the presence of left-right mixing either through the muon Yukawa 
coupling or smuon left-right mixing. Consequently, the degree of left-right mixing in the 
first term in Eq. (1951) is constrained to some degree by the experimental results for a^. 
Taking these considerations into account, the authors of Ref. 9(| find that contributions 



to g^ioop as large as a few xlCT 4 are possible, with the largest effects occuring when 
one of the two smuon mass eigenstates becomes light. Contributions of this magnitude 
would imply that either lies well above the electroweak scale or that all but the SM- 
like, lightest Higgs boson would decouple in order to avoid the presence of charge and 
color-breaking minima in the scalar potential. 

From the standpoint of the Michel parameters, g RR contributes quadratically to the 
combination of parameters that governs the spectral shape and spatial asymmetry 

1 - £ = + \\9 S RR ? + \\9 S lr - 2<?L| 2 (97) 

as well as the parameter £' that enters the energy- dependence of the outgoing lepton 



longitudinal polarization. Experimentally, one has 9l|, |92|, |93 



P£- = 0.99787 ±0.00082 (98) 
f = 1.00 ±0.04 (99) 
leading to \g RR \ < 0.067 at 90% confidence from the recent global fit of Ref. [71[ under the 



assumption that P M = 1. Given the maximum magnitude of g RR loop it appears unlikely 
that one will probe SUSY contributions using these quantities. In contrast, the parameters 
T], if', and P'/A carries a linear dependence on g RR i oop - The impact on the parameter 
7] that characterizes the energy-dependence of the isotropic outgoing lepton spectrum is 
particularly interesting: 

V 2\9 V LL \ 2 ' 1 j 

where the "± ■ • • " indicate contributions from the other gj that are not generated in the 
MSSM. The present limit on rj obtained in Ref. [701 is about two orders of magnitude above 



the SUSY expectations, and a substantial improvement in precision would be needed to 
probe this parameter at an interesting level from the standpoint of SUSY. 

While a direct probe of g RR loop via measurements of the Michel parameters may be 
challenging, its contribution to rj could be large enough to affect the extraction of G M from 
Tp. As indicated by Eqs. (I92|ll00p . the fractional shift in due to this parameter is 

AG^ _ m e ^ e 9 ll9rr, loop . , 

G n "V \9ll\ 

so that contributions to g R R to of order a few xlO -4 would lead to ppm effects in the 



value of the Fermi constant. Given the objectives of the new PSI experiments |72l. l73l|. this 



36 



correction may have to be considered if future collider experiments discover superpartners. 
In principle, one might also probe this SUSY-induced non-(V — A) <g> (V — A) operator 
using Eq. (IHTj) if the uncertainty in Mz and the weak mixing angle can be improved by one 
and three orders of magnitude, respectively. While the latter appears to be an especially 
daunting task, obtaining a commensurate reduction in the theoretical uncertainty in Af 
would likely be even more difficult. 



B. Semileptonic Decays of Light Quark Systems: General Considerations 

As with the analysis of /z-decay, the theoretical interpretation of semileptonic decays 
requires careful attention to electroweak radiative corrections. Moreover, the presence of 
low-energy strong interactions introduces additional complications not present for purely 
leptonic decays. Nonetheless, considerable progress has been made in reducing the theo- 
retical uncertainties associated with non-perturbative QCD effects, allowing one to derive 
information on SUSY and other possible new physics scenarios from precise studies of these 
decays. 

In analyzing both the 0{a) electroweak corrections and non-perturbative QCD effects, 
it is useful to return to Eq. (1661) . replacing the quark spinors by the corresponding field 
operators: 



V2 



Semileptonic = ~ £ (1 + ^ - Af M ) e 7 A (l - J B )u e « TA (l - J 5 )d + h.C. . (102) 



The computation of decay observables requires taking matrix elements of the effective oper- 
ator between initial and final hadronic states. Note, however, that the 0(a) correction Af ^ 
can only be computed reliably in perturbation theory down to a scale A na d at which strong 
interactions between quarks become non-perturbative. Non-perturbative effects arising 
from QCD dynamics below this scale lead to a dependence of the radiative corrections on 
the structure of the initial and final hadronic states as well as on the spacetime properties 
of the hadronic charged current. It is conventional to include these pro cess- dependent 
QCD contributions to the 0(a) radiative corrections in pro cess- dependent corrections. To 
this end, we write the decay amplitudes for various semileptonic processes of interest here 
as 



MJ 2 = 


G ii 

~T2 Vud 


(1 + A# 


- AfJ £ 7 A (1 


-7s) 


w t (0| u^xlbd 7T~) 


M% = 




(1 + Af J 


- Af J e 7 A (l 


-15. 


We (f\ ui\d \i) 


M" GT = 


~T2 ud 


(1 + Ar£ 


- AfJ e 7 A (l 


-lb] 


)f e (f\u f yyysd\i) 



(103) 

where Adp {Mq T ) denote the Fermi (Gamow- Teller) amplitudes for nuclear, neutron, or 
pion /5-decay that involve hadronic matrix elements of the charged vector (axial vector) 
current and M.\ 2 denotes the amplitude for pion leptonic decay. The Af J, Arg, and 
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Ar^ denote the corresponding process-dependent 0(a) radiative corrections 18 . Since we 
have factored out the 0(a) contributions to the amplitudes explicitly, the hadronic matrix 
elements appearing in Eqs. fl 103[) involve purely strong interaction dynamics. 

For economy of notation, it is also useful to define a set of process-dependent Fermi con- 
stants that encode the 0{a) corrections and information on the hadronic matrix elements: 



G7T 

A 



G/jVud (1 + Arjf - ArJ 
G,V ud (l+ArJ-Af M )^(0) 
GfjYud (l+Arf-AfJ)g A (0) 



(104) 



where gv{q 2 ) and g A {<l 2 ) are the nucleon vector and axial vector form factors, respectively, 
defined in Eqs. (11401) below. The pion decay matrix element appearing in M.J 2 a l so contains 
a dimensional form factor, F v {q 2 ), but it is conventional to keep the dependence of the decay 
rate on F n = F n (m%) = 92.4 MeV explicit rather than absorbing it in G\. 

IIV BL 1 Pion Decay Leptonic Decay 



The purely leptonic channel tt + — > fi + u('~f) yields a value for the pion decay constant 
F n that provides a key input for the analysis of chiral dynamics. Theoretically the SM 
prediction for r[7r — ► ^(7)] has been computed to one-loop order (30|. In contrast to the 
situation for /i-decay, where QED corrections have been factored out before extracting 
from r M , the correction Ar A contains QED corrections and the associated infrared 



divergences, 
decay rate 



Consequently, one considers the SM prediction for the total, infrared finite 



r[7r+ -> e + p,(y)} = t\w + -> tv,\ + r[7r+ -> tv a \ 



ml, 



mi 



(105) 



+ brem 



where the "+ brem" indicates 0{a) bremstraahlung contributions. 

The most recent analysis of r[7r~ — ► has been carried out by Marciano and 

Sirlin, yielding the result (30f 
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z inn 
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7T J L 
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7T /X 



(106) 



x 



1 + -F(x) 

7T 



where the quantities proportional to a arise from both the electroweak and QED radiative 
corrections entering G\ and contributions from real photon radiation. Note that the cor- 
rections are manifestly process- and hadron structure-dependent. The quantity A x = 47rF n 



18 The value of the correction ArY for the neutron differs from the corresponding correction in pion /3-decay. 
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is the scale of chiral symmetry breaking associated with the onset of nonperturbative dy- 
namics; the Ci(ji) denote low energy constants that parameterize presently incalculable 
non-perturbative QCD effects and that depend on the renormalization scale /i; F is a cal- 
culable function of x = mj/m^; and the + • • • indicate additional terms that are suppressed 
by m|/A^. Both the function F(x) and the terms containing the Ci(fi) arise from QED 
corrections to the decay rate for a point-like pion 19 . 

The large logarithms ln(M^//i) in Eq. f)105|) have been resummed to all orders in powers 
of [(a/n) ln(Mz//i)] 2 using the renormalization group, yielding the electroweak correction 
factor Sew{^, Mz) that replaces the factor 1 + 2( a/7r ) ln(Mz//f) (30j| ■ For fi = m p one 
has SE\v{ m pi Mz) = 1.0232. The authors of Ref. |30J estimate that the uncertainty in 
T[ti — ► y,v{'y)\ associated with the Ci(fi) is ±0.56%. This uncertainty dominates the error 
in F-K since the pion lifetime and leptonic branching fraction are known to ±0.02% and 
±0.00004% uncertainty, respectively. One thus obtains 

F n = 92.4 ± 0.025 ± 0.25 MeV (107) 

where the first error is associated with the value of V u d and the second with the effects of 
the Ci{m p ). 

Since one cannot presently compute F n from first principles with a precision comparable 
to the uncertainties quoted in Eq. (I107p . a measurement of r[7r — > ^(7)] by itself does not 
provide a useful low-energy probe of SUSY effects. On the other hand, inclusion of SUSY 
corrections to the rate could alter the extracted value of F n . To illustrate this possibility, 
it is convenient to write Eq. (11061) as 



r[7T+ - £ + p e ( 7 )} 
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x 



{1 + (2 [Mi - ArJ + brem) SM + 2 (At? - Af,) ncw } 



where (2[Af^ — Af p ] + brem)sM denotes the 0(a) SM corrections to the rate appearing 
on the RHS of Eq. (11061) and 2(Af^ — Af M ) new indicate the corrections from new physics. 
Since the latter generally involves the exchange of heavy particles, one does not encounter 
infrared singularities in the new physics corrections and need not consider the corresponding 
bremstraahlung contributions. 

To illustrate the impact of these corrections in SUSY, we first consider the case of RPV 
interactions, neglecting left-right mixing among sfermions. In terms of the Aijk(f) and 
A' ijk (f) defined in Eq. flUD, we have 0, [95 



(Af£ - AfJ 



RPV 



Ai2fc(e 



R) 



(109) 



where the subscript "f denotes the generation of the final state leptons. From a fit to 
other low energy observables discussed in Section III Fl we obtain 

(a) -0.004 < (Ar^ 



(b) -0.003 < (Af~ 



Af ) < 
Af M ) ncw < -0.0004 



-0.001 95% C.L. 



95% C.L. 



;no) 
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We have adopted a slightly different normalization convention from the one used in Ref. 30 1 
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for the partial rate involving final state muons, where we have required the Ay fc , A[- k to be 



non- negative according to the definition of Eq. (|48[) . Here (a) or (b) refers to the different 
value of S\ V^d| 2 /| Vud\ 2 in Table. IHII that was used in the fit. As a result, we obtain a small 
shift in the central value for the rate (11061) and an additional uncertainty of roughly ±0.5%, 
with a corresponding ~ quarter percent increase in the uncertainty of F n . It is interesting 
that the estimated hadronic structure uncertainty associated with SM radiative corrections 
is comparable to this RPV effect. 

In order to obtain a probe of new physics using tt£ 2 decays, one must attempt to circum- 
vent the uncertainties associated with the hadronic matrix element that is parameterized 
by F n . To that end, it is useful to consider the ratio of partial rates for final state electrons 
and muons 



R 



r 



7T 



e/fi 



Ml)] 



whose measurement provides a useful test of lepton universality of the CC weak interaction 
and its breakdown due to 0(a) corrections or new physics. From Eqs. (11061) and (11081) it 
follows that 



R e/ , = ^ [4r^rJ {l+(2[A^(e)-Af^)]+A brem ) SM (112) 

+2 [A^(e)-A^(/i)] ne 

= R™{l + 2[Art(e)-Ar^)] n J 

where the Ar^(£) indicate the corrections for final state lepton £, the 'A brem " indicate 
the difference in the bremstraahlung contributions to the rate for final state electrons and 
muons, and R^j^ is the SM value for the ratio. The SM contributions have been computed 
in Ref . M 



pSM _ ^\ 

~ ml 



m; 



mi 



m 



7T 
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(113) 

The terms proportional to am\ are numerically insignificant and have been omitted from 
this expression. To obtain a precise numerical prediction, the authors of Ref. [301 included 
structure-dependent bremsstrahlung corrections and performed a renormalization group 
resummation of all [(a/V) ln(m e /m / J] n corrections that enter F(m e /m n ) — Fim^jm^^ 
leading to 

R s e Jl = (1.2352 ± 0.0005) x 10~ 4 . (114) 

where the error is dominated by theoretical uncertainties in the structure dependent 
bremsstrahlung contributions. 

The ratio R e /, has been measured by groups at TRIUMF 64] and PSI (65|, yielding 
the world average 



E>cxp 

R s f 

e/n 



0.9966 ±0.0030 ±0.0004, 
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where the first error is experimental and the second is the estimated theoretical uncertainty. 
A new measurement that has been approved to run at TRIUMF aims to reduce the exper- 
imental uncertainty to a level comparable with the theoretical error [8~2| . A measurement 
at PSI with a similar improvement in precision has also recently been approved[83]. 

The new measurements could provide significant tests of possible SUSY contributions 
to the breaking of CC lepton universality. In the case of RPV contributions, for example, 
one has 

[A^(e)-A^(/i)] Rpv = A' llfc (4)-A' 2U ,(4) . (116) 

Thus, the precise measurement of R e /^ provides an important input into the global analysis 
of the RPV corrections A^ fc , etc. (see Table UTTj) . The global analysis leading to Eq. (11 10p 
uses the present average (1115)) . and the new measurements should reduce the range on the 
RPV-related uncertainty in the value of F n . 

When Pr is conserved, lepton universality can be broken by superpartner loop cor- 
rections when the first and second generation sleptons have unequal masses. SUSY loop 
corrections to the Wud vertices, light quark propagators, and VT-boson propagator are 
identical for the tc + — > e + z/ e (7) and n + — > // + 1/^(7) amplitudes, thus canceling out of the 
ratio R e /n at O(o). In contrast, corrections to the Wtv t vertices, lepton propagators, and 
box graphs can differ in the two cases, so that R e /^ provides a probe of the non-universality 
of the first and second generation sleptons. An analysis of these corrections has recently 



been performed in Ref. 96]. In the limit that the charginos and neutralinos are approx- 
imately pure gaugino and higgsino states (a limit achieved in the absence of electroweak 
symmetry-breaking), the SUSY correction to R e /^ is dominated by box graphs. The vertex 
and external leg corrections for the W + £ + i>£ vertex sum to zero as required by the Ward 
Identity 20 . The resulting correction is 



A r>SUSY 

XJdSUSY uv ud 1 iv±w \ I F I "L\ -ft \ _ F l ""L 2 "ft I I / 11? N 

orL e/n j^SM , ■- > \ 1 r I 1 I i rl ' i r2 I I i rl ' \ rl \ ( ' 




where F(x,y) is a loop function associated with the box graph with F(l,l) = 1 (corre- 



sponding to m^. = tUq = M 2 ). When the two sleptons are equal in mass, Ai?fF = 0, 



whereas when m^^ >> mi., one has 

\5R™ SY \ = 1.7 x 1(T 3 [ ^ ] F[-^:^k\ , (118) 

so that the correction can be as large as a few x 1CT 3 in this case. 

Allowing for significant gaugino-higgsino mixing can lead to can lead to logarithmic 
enhancements of SRfj 1 ^ as the vertex and external leg corrections no longer sum to zero 




20 When D R is used, the SM and superpartner corrections individually contribute //-independent constants 
of equal magnitude and opposite sign. 
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in this case. The resulting contribution is 



e/ti \ V+L 47rs 2^ y m 2 



2 

, 2-2V* 1 U* 1 N i2 N l2 



-V* 1 U; 2 N l2 N i3 /V2 + U; i V* 2 N* 2 NyV2\ + (119) 

where the pre-factor in the square brackets that contains the chargino and neutralino 
mixing matrices vanishes in the limit that M\ i2 and A* are niuch larger than My/ and and 
can be as large as 0.5 for gaugino and higgsino mass parameters of order 100 GeV. The 
resulting correction can, thus, be as large as 

\5R% SY \ < 1.3 x 10- 3 ln ( ^ ) . (120) 




A numerical scan over MSSM parameters performed by the authors of Ref. 1961 is consistent 
with this bound. Given that the new TRIUMF and PSI experiments hope to reduce the 
experimental uncertainty to the level of the present SM theory uncertainty (~ 5 x 10 -4 ), one 
could expect to see a departure from the SM expectations of several standard deviations if 
the masses of the selectron and smuon differ by more than a factor of ~ 2 and the mass of 
the lightest chargino is less than 300 GeV. Moreover, in the case of such non-degeneracy 
among sleptons, the sign of 8R^ SY provides an indication of whether the lightest first or 
second generation slepton is heavier. 

IIV B1 2 Neutron and Nuclear (3 Decay: General Features 

Studies of nuclear and neutron /3-decay have yielded both important information about 
parameters of the SM as well as constraints on physics beyond it (for comprehensive re- 
views, see Refs. (62l . l66l In particular, the most precisely-known element of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix that characterizes the misalignment of quark 
weak interaction and mass eigenstates in CC interactions has been determined from "su- 
perallowed" nuclear /3-decays. The latter involve transitions between spin-parity J 71 " = + 
nuclear states that are mediated solely by the vector charge component of the CC weak 
current. To the extent that the initial and final + nuclear states are states of pure isospin 
and the energy transfer to the outgoing lepton pair is negligible compared to typical nuclear 
scales, the transition matrix element is independent of nuclear structure. In contrast, the 
rates for nuclear decays involving initial and/or final states with non-zero spin - including 
the decay of the neutron - depend more strongly on the details of hadronic and nuclear 
structure. In these cases, the extraction of precise information on the electroweak sector 
of the SM or on new physics generally requires measurement of both the decay rate and 
one or more decay correlation coefficients. 

In all cases, the overall decay rate in the SM is characterized by the the reduced half 
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life or u ff value, given by 
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(G V ) 2 M 2 + (G^) 2 M 2 T 
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h(2n 3 \n2){hcf / {m e c 
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where t is the half-life; / is a factor that takes into account the outgoing /3-particle wave- 
function in the presence of the nuclear Coulomb field; the Fermi matrix element Mp 
is the nuclear transition matrix element of the vector current charge operator Jq(x) = 
u'(x)d(x) + h.c, while the Gamow- Teller matrix element Mqt involves the spatial compo- 
nent of the axial vector current operator 21 . 

It is often useful to consider the differential decay rate for the decay of a polarized 



nucleus, for which one has|93l (see also Refs. 



62|,|6f 



621) 
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[122) 



where N(E e ) = p e E e (E — E e ) 2 , E e (E v ) and p e (p u ) are the (3 (neutrino) energy and mo- 
mentum, respectively; E is the endpoint energy; and J is the polarization of the decaying 
nucleus, a is the (3 polarization, and V = a/1 — (Za) 2 . The coefficients of the various 
correlations involving lepton momenta and nuclear spin depend on the structure of the 
underlying lepton-quark weak interaction. The correlations parameterized by A, B, and 
G are odd under parity (P) and even under time- reversal (T); the D-term is T-odd but 
P-even; the R correlation is both T- and P-odd; and all others are P and T-even. 

The various correlation coefficients in Eq. f)122p carry complementary dependences on 
the ratio of Fermi constants 



P 



y _G" A _ 9a(0) ( 
G^ v 9v(0) [ 



(123) 



as well as on operators that depart from the (V — A) <g> (V — A) form of the SM, CC low- 
energy current-current interaction. In the SM, the correlation coefficients can be expressed 
in terms of A alone. In the case of neutron decay, for example, one has: 



1-A 2 
1 + 3A 2 ' 



A = -2 



A(l + A) 
1 + 3A 2 ' 
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' 1 + 3A 2 



(124) 



21 For neutron decay, Mp and Mqt simply indicate the matrix elements of the vector and axial vector CC 
operators without respect to their spacetime components. 
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with analogous expressions for the coefficients N and G. The quantity b appearing in the 
so-called Fierz interference term is zero for purely vector and axial vector interactions, 
while N vanishes for pure vector transitions. The T-odd correlations are zero in the SM 22 . 

As with /i-decay, it is useful to describe non-(V — A) ®{V — A) possible departures using 
an effective, low-energy Lagrangian. In analogy with Eq. (193]) we write 

AC 

£/3-dccay = ^ ^ Q 7 gjv^ -p^ _ (125) 

The SM gives 

&LL — Vud (1 + Af/3 — Af M ) , (126) 

with all other 5 = 0. Note that in hadronic matrix elements of the SM, (V — A) ® ( V — A) 
operator, non-perturbative QCD effects will lead to a renormalization of the axial vector 
current and to differences in hadronic contributions to the radiative corrections ArY and 
Arj as indicated in Eqs. 01030104]) . 

Tree-level, supersymmetric contributions to a\ L can be generated by RPV interactions. 



One has 94, 951. 



Ar^' A - Af n 



RPV 



A' u& (4) - A m (4) . (127) 



Note that because the RPV interactions preserve the (V — A) <g> (V — A) structure of the 
low-energy CC interaction, the RPV contributions to ArY and Ar A are identical and, thus, 
do not affect the value of A. As we discuss below, probes of these contributions can be 
obtained using a combination of /?-decay, kaon semileptonic (K^) decays, 5-meson decays, 
and the unitarity property of the CKM matrix. The results of these unitarity tests have 
been used as input for the global analysis of RPV corrections A^, A^. fc , etc. summarized 
in Table EU 

Supersymmetric loop corrections can also contribute to [Af^' A — Af M ]. As with the 
RPV corrections, loop SUSY loop corrections to ArY and Atq are identical in the MSSM, 
since the superpartner mass scale lies well above the hadronic scale. Hence, the presence of 
these corrections leaves the value of A unchanged from its SM value. An analysis of these 



corrections has been carried out in Ref. [98| using the MSSM with R-parity conservation. 
Doing so requires computation of the corrections illustrated in Fig. H] (see Section [ITT]) . In 
general, carrying out an analysis of these corrections on the vast number of MSSM param- 
eters is a formidable task. Typical analyses resort to one of two strategies to contend with 
the large number of parameters: either one carries out a model-independent analysis by 
randomly generating a large number of MSSM parameter sets and computing the correc- 
tions for each set, or one reduces the number of independent parameters by resorting to 
a model for SUSY-breaking mediation in which the soft SUSY parameters are determined 
from a small number of parameters at a high scale and their RG running to the weak scale. 



Final state QED interactions can induce non- vanishing contributions to D that mimic the effects of bona 
fide T-violation. 
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In the case of the corrections to Gy, however, cancellations between corrections to AfY'" 4 
and Ar^ allow for an alternate approach. Specifically, contributions to the P^-boson prop- 
agators - flyyw ~~ cancel entirely between the two terms 23 , as do corrections to the first 
generation lepton propagators and Wev e vertices. As a result, [Af^' A — Af M ] S usY-ioo P de- 
pends only on the differences in the W u v v , and W ud and external leg corrections as well 
as the box graphs. The authors of Ref. [98( found that these simplifications allowed comple- 
tion of a model-independent analysis without resorting to randomly generated parameter 
sets. We review the results of this analysis below. 

One-loop radiative corrections in the MSSM may also induce non-zero scalar and tensor 
interactions via the box graphs of Fig. [8] (b). These loop-induced non-(V — A) ® {V — A) 



interactions have recently been analyzed by the authors of Ref. [90(| , who obtained 

$° = \U kl \ 2 Z^Z^Z]-Z'r \N n ?Ti (M x o,M xp M^M Jm ) (128) 

~8f = - a ^ V * U jl V£Zh i *Z«Zl m Z?» \N kl \ 2 M xt M x oF 2 (m x +, M x o, M lm ) 
where the notation is similar to that of Eq. fl95l) . The corresponding operator coefficients 



are 

4r = (129) 
n s - -?n T - P ] 

Note that, in contrast to the situation with /x-decay, the presence of SUSY loop-induced 
scalar and tensor interactions relevant to /3-decay requires only flavor- diagonal L-R mix- 
ing among scalar fermions. These loop-induced scalar and tensor interactions generate 
contributions to the Fierz interference parameter, b, of Eq. (11221) ; the energy-dependent 
components of the neutrino asymmetry parameter B and spin-polarization coefficient Q'; 
and the energy-independent component of the spin-polarization coefficient N. We discuss 
the prospects for future /3-decay experimental probes of these loop-induced interactions 
below. 



C. Superallowed Nuclear Decays 

For the superallowed Fermi transitions one has Mqt = 0, leaving only the dependence on 
Mp in the /t-value. To high degree of accuracy, Mp is independent of the details of nuclear 
structure, making these transitions an ideal venue for determining Gy and, thus, V u d- As 
Eqs. (11041) and (I12ip indicate, the determination of V u d from these transitions requires 
both experimental and theoretical input. Experimentally, the ft values for twelve different 
transitions have been measured with uncertainties ranging from ~ (3-25) x 10" 4 |6l|, 



Up to negligible corrections arising for the different kinematics of each process. 
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while the value of CL is known to ten ppm accuracy (assuming rj = as in the SM). As 
discussed in Refs. @8, a determination of the experimental half lives requires three 



distinct experimental measurements: the total decay half life, the branching ratio for the 
decay to the + ground state of the daughter nucleus, and the energy release in the decay, or 
Q-value. New measurements of the Q-values for several superallowed decays ha ve re cently 



been completed, leading to shifts in the ft values in some, but not all, cases [75|, [100J. The 
impact of these new measurements on the overall fit to the twelve decays awaits completed 
analyses of Penning trap measurements of the Q- values of 26m Al and 42 Sc. 

Theoretically, one requires computations of the matrix element Mp as well as the SM 
radiative correction factors Af^ and Af^. In the limit of zero energy transfer and exact 
isospin, the Fermi matrix element is 

M F = (I,I z ±l\J \I,I z ) = [(I T Iz)(I±Iz + l)} 1/2 ■ (130) 

For the most precisely-known cases, one has I = 1,I Z = so that M F = y/2. For 
realistic nuclei and finite energy transfer, one must apply small nuclear structure-dependent 
corrections 

M|(l + S R )(1 - 5 C ) = [(I =F I Z )(I ±I Z + 1)] (131) 

where 5c is a correction that accounts for isospin-breaking and 6r is a nucleus-dependent 
contribution to the 0(a) electroweak radiative corrections. After applying these one ob- 
tains a corrected ft value 

Tt = ft(l + 5 R )(1-5 C ) ■ (132) 

It follows from Eqs. ( 1121ft and ( 11321) that Tt should be the same for each nucleus, since this 
quantity depends only on the nucleus-independent Fermi constant G v and the universal 
matrix element in Eq. (11301) determined solely by the isospin quantum numbers. For 
historical reasons, this prediction is described as a consequence of the conserved vector 
current (CVC) property of the semileptonic, CC weak interaction. The corrections 5r,c 
have been computed by the authors of Refs. (6ll . 99| using the nuclear shell model (NSM) 



and applied to the twelve best-known superallowed transitions, leading to an average Tt 
value 

Tt = 3072.7 ± 0.8 s (133) 

where the error includes the estimated theoretical uncertainty associated with the correc- 
tions 5 R and 5c 24 - This result for Tt implies consistency with CVC at the 0.026% level. 

The result in Eq. ( 11331) provides a test of the inter-nuclear consistency implied by CVC 
but does not account for the possibility of a nucleus-independent systematic theoretical 
error associated with the calculated nuclear corrections. To allow for this possibility, 



the authors of Refs. [6l|, |99| compared the NSM Tt valu es w ith those obtained by us- 



ing the 5c Hartree-Fock calculations of Ormand and Brown [10 11 ] . Since the latter exist for 



nine out of the 12 most accurately measured superallowed transitions, Towner and Hardy 



The x 2 P er degree of freedom for the average is 0.42. 
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(TH)averaged their values with those of Ormand and Brown (OB) for these cases, yielding 
the average 

Tt = 3073.5 ±1.2 s (134) 

where an additional theoretical uncertainty given by half the difference of the TH and OB 
values for Tt. In what follows, we use the average (11341) in the discussion of Vud- 



A value for G v can be extracted from Tt by employing Eqs. (1121 |l 1321) and subsequently 
used to determine V u d by applying the SM radiative correction difference Af # — Af „ ap- 



pearing in G v . The dominant uncertainty entering the latter is the theoretical uncertainty 
associated with hadronic contributions to the correction Af ^ arising from the W'y box 



graph. To O(a) one has [28 



AfJ(iy 7 box) 



a 
8^ 



M4|h+c>(A) 



(135) 



Here, the leading logarithmic term is generated by short-distance contributions to the loop 
integral and can be computed reliably in the SM. The the constant C 7 vi/(A) parameterizes 
contributions to the loop integral below a momentum scale A. An estimate of C 7 vk(A) was 
given by Marciano and Sirlin [29( using nucleon intermediate states in the box diagram, 
and this estimate had been retained by subsequent authors for many years. An estimate 
of the uncertainty in this contribution had been obtained by varying A over a reasonable 
range, yielding an uncertainty of ~ ±0.00038 in Ar^. 

Recently, Marciano and Sirlin have reduced this uncertainty by a factor of two over 
previous estimates by relating the asymptotic part of the W'y box integral to the Bjorken 
sum rule to include perturbative QCD corrections through C(a^) and by employing large 
A^-based current-current correlators to treat the resonance region [3l|. The result is an 
uncertainty of ~ ±0.00019 in Ar . Using the results of this new analysis and an up-date 



i — I i—i /3 ■ , 

of the fit in Refs. [61|, [99j by Savard et al [100] one has 



1^ = 0.97377(11) (15) (19) (136) 

where the first error arises from combining the experimental error in Tt and nuclear struc- 
ture theory uncertainty; the second error is associated with nuclear Coulomb distortion ef- 
fects; and the final error is the theoretical hadronic structure error associated with C 7 w(A). 

In light of this new result, the uncertainty associated with Af^ and with the nuclear cor- 
rections 8r,c are comparable. In order to test theoretical, nuclear structure computations 
of these corrections, new measurements of additional superallowed decays are underway 
with nuclei in which the magnitude of the calculated nuclear corrections are larger than the 
corrections for the nine most precisely measured transitions used to obtain the result (11361) . 
These studies involve even-Z, Iz = —1/2 parent nuclei in the mass range 18 < A < 42 
and four odd-Z, Iz > 62 + parent nuclei. Obtaining refined measurements of half lives, 
branching ratios, and Q-values for both series of decays is challenging, though use of Pen- 
ning trap techniques have allowed for new, highly precise mass measurements for the light 
cases. A more exte nsive review of these challenges and prospects for meeting them can be 
found in Ref. [l02 . 
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In addition to providing a determination of Gy, studies of superallowed transitions also 
provide a probe of the Fierz interference coefficient bp (the subscript "F" denotes those 
contributions allowed for Fermi decays). In terms of the effective operator coefficients aj s 
one has 

bF = ±lllRe ( a nL + a RR ) (137) 

9V V a LL J 

independent of the details of the nuclear matrix elements (the upper and lower signs cor- 
respond to P~ and (3 + decay, respectively). Here, gy and g$ are the nucleon vector and 
scalar current form factors, respectively [see Eq. (|140p below]. An global analysis of super- 
allowed decays leads to the &f = 0.0026(26) [6l|, [99|], implying stringent bounds on scalar 
interactions. In order to probe the non-(V — A) ® (V — A) interactions generated by 
superpartner loops, roughly an order-of-magnitude improvement in sensitivity would be 
required. Specifically, from Eqs. (I128|I129I) we have 

bp = ±^ ( 9S\ ReZ )m Z Am* 



3tt \g v 



U kl \' Z?Z$ \N jt \ A M Z T i ( M x o, M xt , .1 /,,..!/, 



-UjiVZZpZg \N kl \ 2 M 2 z M xt M x ^ 2 [M xp M xl ,M^ M^)] . (138) 

The prefactor 2a/37r is (9(10~ 2 ) while the product of M§ and the loop functions can be as 
large as 10 _1 for fh ~ Mz- For nearly maximal L-R mixing among sfermions, the product 
of rotation matrices ZjfZf?* etc. is 0(1). Thus, contributions to bp as large as ~ 10~ 3 can 
occur in the regime of maximal mixing and superpartner masses of order the weak scale. 
Measurements designed to observe bp at the few xlO -4 level would provide interesting 
probes of L-R mixing among first generation scalar fermions. Note that non-vanishing 
results at this scale would disfavor the alignment hypothesis of Eq. (fTTj) . which, for of 
order the electroweak scale, implies that L-R mixing for the first two generations is Yukawa 
suppressed. As in the case of the muon decay parameter g^ R , large L-R mixing for the 
first generation also implies that the masses of the Higgs bosons H°, A , and are super 
heavy, leaving only the light SM-like Higgs h° as an experimentally accessible degree of 
freedom. On the other hand, null results would provide added experimental plausibility to 
the alignment idea. 



D. /3-Decay Correlations 

Experimental studies of the (3 spectral shape, angular distribution, and polarization can 
provide information on both non-(V — A) <g> (V — A) interactions as well as an alternate 
means of obtaining V u d- As a specific illustration, we consider polarized neutron decay, for 
which both the vector and axial vector components of the CC contribute to the ft value 
in Eq. ffT2ll with 

M F = 1 (139) 
M 2 GT = 3g 2 A 
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and where gA is the hadronic axial vector coupling that characterizes the strong interaction 
renormalization of the axial vector quark current that enters the neutron decay matrix 
elements. The latter are given by 



(p|ti(0)7„d(0)|7i) = U P (P') 
(p\u(0)wd(0)\n) = U P (P') 



igM{q 2 



2,171 7\r 



gpjq 2 

m N 



<V9 



q»i5 



u n {P) 

Un(P) 



(140) 



At q 2 = one has gy(0) = 1 and <?m(0) = Kp - k„ according to the CVC property of 
the vector CC, whereas gA = <?a(0) ~ 1.26 and gp ~ 8.5 (Here, we quote a value for gp 
taken from chiral perturbation theory that is in agreement with the results of ordinary 
muon capture experiments but that differs fro m the result obtained from radiative muon 

Neglecting the small g 2 -dependent 
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capture. For a recent discussion, see, e.g., Ref. 
corrections associated with the gu and gp terms, Eq. (I140p leads to the matrix elements 
in Eq. (11391) . Nucleon matrix elements of scalar and tensor operators, ud and uo^d, 
associated with non-SM interactions are parameterized by analogous form factors, g$ and 
gx, respectively. 

At present, it is not possible to compute with the 0.1% precision needed for a 0.1% 
determination of V u d from the neutron lifetime, r n . Consequently, the axial vector contri- 
bution to ft in Eq. (11211) must be separated experimentally from the vector contribution. 
Doing so requires measurement of a neutron decay correlation coefficient appearing the par- 
tial rate Eq. (I122p . The coefficients A, a, and B carry a dependence on A [see Eq. (1124p ], 
so that their measurement can yield the requisite separation. The most precise value of A 
has been obtained with a 0.6% measure ment of the A parameter by the PERKEO collab- 



oration, leading to A = —1.2739 ± 0.0019 104]. Since the publication of that result, a new 



value for r n has been obtained at ILL that differs from the previous world average by more 



than six standard deviations 76] . Including that result and performing a one-parameter fit 
to neutron decay measurements yields A = —1.27293(46). The resulting value for V ud is 



V ud = 0.97757(65) 



compared with the value V u d = 0.97192(65) obtained using the previous world average for 
r n and the new average for A. 

Several efforts are underway in order to obtain a value of V u d using a combination of the 
neutron lifetime, r n , and correlation coefficient measurements, and a review of these studies 
can be found in Ref. [5j. In light of the new r n result, additional precise determinations of 
the neutron lifetime are either underway or are being planned at ILL, NIST, and LANSCE. 

Precise measurements of neutron decay correlation coefficients may also probe the SUSY 
loop-induced non-(V — A) <g> (V — A) interactions. For example, the (3 energy-dependent 
component of the neutrino asymmetry parameter B depends linearly on both the scalar 
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and tensor interactions of Eq. (j!29p : 



B 



SUSY box 




(142) 



As with the Fierz interference term, the SUSY contributions to B can approach the 10~ 3 
level for nearly maximal L-R mixing. Experimentally, future measurements using cold or 
ultracold neutrons may be a ble t o determine the energy-dependent component of B with 
a sensitivity of a few x 10~ 4 [l05l ] . At present, the results of nuclear /3-decay correlation 
measurements do not appear to be sensitive to the linear interference of scalar and tensor 
interactions with the SM amplitude; improvements in superallowed sensitivity to bp and in 
neutron decay correlation measurements appear to hold the brightest prospects for probing 
the SUSY loop- induced non-(V — A) ® (V — A) interactions. 



E. Pion /3-decay 

In addition to the use of nuclear and neutron decay, measurements of the rate for pion 
/?-decay (irp) also yield a determination of G v . In this case, the non-universal hadronic 
contributions to the radiative correction Afp for Tip differs from that for neutron and nuclear 
decays, and in the past it has been argued that the corresponding theoretica l unc ertainties 
for TTp are smaller. However, a recent analysis using ;\PT reported in Ref. FloI quotes a 



theoretical uncertainty in V u d of ±0.0005 - a value that is larger than the new theoretical 
uncertainty associated with neutron and nuclear decays. 
The rate for pion /3-decay is given by 

r (^) = g^a > ( 143 ) 

where f±(t) are the two pion form factors and /(A+) is a phase space integral that results 
from inclusion of real photon emission and that is a function of the slope A+ of /+(£) at 
the photon point. The non- universal, long-distance 0(a) corrections can be included by 
replacing /+(£) by 

Fl(t, u) = [l + £l>, ml ml, A 7 )] , (144) 

where A 7 is an infrared regulator whose effect on the total rate is cancelled by the corre- 
sponding A 7 -dependence of /(A+). 

Contributions to T(u,m^,ml, A 7 ) that are non-analytic in the various masses and mo- 
menta can be computed unambiguously at one-loop order and carry no dependence on a 
priori unknown parameters. However, there exist analytic contributions that arise at the 
same chiral order that are parameterized by three constants in the effective Lagrangian: 
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Kl 2 (fjL)j Xi, and Xq(/i). Here, K{ 2 and Xg depend on the renormalization scale \i since 
one-loop graphs generate divergences having the same structure as the corresponding terms 
in the Lagrangian. A theoretical prediction for K[ 2 (m p ) with ~ 10% uncertainty has been 
given in Ref. [l07 ], while bounds on X\ and Xg(m p ) were obtained by the authors of 



Ref. [106] using dimensional analysis. The uncertainty associated with these constants 
dominates the theoretical error in the extraction of V u d from r(irp). 

Experimentally, the PIBETA collaboration has recently obtained t he m ost precise deter- 



mination of the ratio of branching ratios for the 7rg( 7 ) and 7r e2 ( 7 ) decays [108l ] . Multiplyin g by 
the current world average for the 7r e2 ( 7 ) branching ratio leads to the irp branching ratio [63] 

S^ (7 ) = [1.036 ± 0.004(stat) ± 0.004(sys) ± 0.003(vr e2 ( 7 ))] x 10" 8 (145) 

for a total experimental error of ±0.006 x 10 -8 . Note that the fractional uncertainty is 
roughly a factor of ten times larger than the theoretical error and fifteen times larger than 
the combined experimental and nuclear structure theory uncertainty in the superallowed 
J-t value. Thus, considerable experimental and theoretical progress is necessary before a 
value for V u d can be obtained from 7Tg with precision competitive with the superallowed 
value. 



F. Kaon decays and V u , 



In order to use the value of V u d (H36I) as a probe of SUSY, one must compare it with 
the SM expectation. In this case, the SM implies that the CKM matrix is unitary, so that 



|VU 2 + |VU 2 + 



V, 



ub 



SM 



(146) 



The value of V u b = 0.0032 ± 0.0009 is obtained from the decays of 5-mesons. Given the 
level of uncertainty in V u d, both the magnitude of V u h and its error are too small to affect a 
test of Eq. (11461) . On the other hand, the value of V us as well its uncertainty are critically 
important. The theoretically cleanest determination of V us is obtained from the branching 
ratios for decays, K — > 7r£v. The partial rate for this decay mode is given by [63] 



dV{K\ 



G 2 ,m 5 



K 



(.?,) 



1287T 3 



S J wC{t)\V w \ i \f+ (0)1 



1 + 



Aft 



mi 



1 + 2Af : f7 ( 2 ) 



4- 1A KI 1 
+ zia em\ 



(147) 



where f+(t) is the K-to-n transition form factor with t = (j>k — Pn) 2 , A^/m^. is the slope 
of the form factor at t = 0, and C(t) is a kinematic function that depends on the kaon 
form factor, Sew contains the short-distance electroweak radiative corrections, and Agj^ 
and &em indicate corrections generated by the breaking of flavor SU(2) and long-distance 
electromagnetic corrections, respectively [86[]. To carry out a test of CKM unitarity at the 
0.1% level, one must include the A^^ and A^^ corrections and determine the value of 
/+■ (0) with one percent uncertainty or better. 

As with the determination of V u d using the /9-decay half lives, arriving at a value for V us 
from the partial rates requires both experimenta l and theoreti c al input. Experimen- 
tally, new determinations of the branching ratios [l 09 . lid . 111 . 112L Ill3 | - combined 
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with experimental values for and C(t) [l 14 . 115 . 116]- have yield a shift in the world 
average for the product V us x f+(0). The 2005 Particle Data Group value is 



V us [/* (0)/0.96l] = 0.2257(9) 



(148) 



Here the value of f+{0) has been normalized to the c ombi ned chiral perturbation theory 
(xPT)-quark model prediction of Leutwyler and Roos [ll7 | that had been used for many 
years. That prediction includes contributions from non-analytic, one-loop terms through 
chiral order p 4 , analytic 0(p 4 ) terms that can be obtained from fits to other pseudoscalar 
meson observables, and a quark model estimate of the 0(p 6 ) contribution. 



85|. In this 



Recently, the 0(p e ) loop contributions have been computed in Refs. 
context, the dominant, remaining uncertainty is associated with the 0(p 6 ) analytic contri- 
butions that depend on the square of the 0(p 4 ) constant L r 5 ([i) and two of the 94 unknown 
constants appearing in the 0{jP) chiral Lagrangian, C[ 2 (/i) and C^/i): 



K (6), analytic _ Q ( m l ~ 



f- 



F A 



F 2 



c r M-c r M 



+ 



(149) 



where /i is the renormalization scale, usually taken to be ~ A x . The constant L r 5 (/j,) is 
presently well known from experiment, whereas the determination of C[ 2 3 4(/i) require ad- 
ditional input. In principle, new measurements of the pion and kaon scalar form factors 
fo' K (t) could allow a determination of these unknown constants, removing this remaining 
uncertainty. In particular, a 5% measurement of the slope of f$ (t) and a 20% determina- 
tion of its curvature - coupled with a 1% theoretical determination of the ratio of decay 
constants F K /F n from theory - would be sufficient to determine C[ 2 34 (/x) at a level needed 



for the first row CKM unitarity test |63l. 185 



Alternately, values can be determined theoretica lly. R ecent work using large Nc QCD 
has yield a value for C{ 2 34 leading to the prediction [ll^ 



/+ (0)large7V c 



0.984 ±0.012 



(150) 



A number of lattice QCD computations of /+ (0) have been carried out that yield, in 
effect, the sum of the nonanalytic and analytic terms in the analysis. The results 
tend to favor a smaller value for f+(Q) that is consistent with the Leutwyler and Roos 
estimate 117 : 



lattice 



0.960 ± 0.005 stat ± 0.007 ; 
0.962(6)(9) 
0.952(6) 
0.955(12) 



svs quenched, Wilson[119] 

unquenched, staggered [120] 
unquenched, Wilson [121] 
unquenched, doman wall [122] 



:i5r 



where recent quenched and unquenched results are shown for different lattice fermion 
actions (Wilson, staggered, domain wall). Note that the quenched results of Ref . 119| do 



not include any systematic uncertainty associated with the quenched approximation. 
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An alternative determination of V„, s can be made by comparing the rates for the leptonic 
decays of the charged kaon and pion |l23| . From Eq. (I106p and the analogous expression 
for the charged kaon decay rate one has 



r[7r+ -> n+v(ri)\ 



ml 



2\2 



V 2 d Fl m\ 



m 



2\2 



a 



A 



7T 



Kit 
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where Akw gives the difference in the radiative corrections entering the RHS of Eq. ( 11061) . 
The uncertainty from the hadron structure-dependent contributions to this difference has 
been estimated to be ±0.75, corresponding to an uncertainty of ~ 0.1% in the ratio 
{Vus/Vud) 2 ■ Preliminary lattice QCD results for the pseudos calar decay constants obtained 
by the MILC collaboration give F K /F n = 1.198±0.003j£!£!(l23. The resulting theoretical 
uncertainty in V us obtained by this technique is comparable with that entering the analysis 
of Kiz branching ratios. Using a new result for the branching ratio obtained by KLOE 
one obtains V us = 0.2245±Ho3i- 

Using these results and those for V ud from the superallowed decays, one obtains for the 
first row of the CKM matrix 



v ud \ 2 + \v us \ 2 + \v 



ub 



0.9968 ± 0.0014, large jV c [118] 

0.9998 ± 0.0015, unquenched lattice, domain wall [122] 



(153) 



G. CKM Unitarity Tests: Implications for SUSY 

The implications of the CKM unitarity test for SUSY can be significant. In the 
presence of RPV interactions, a CKM unitarity deficit could be remedied by having 
[Afp - AfJ R pv < 0, thereby implying that A m (e|) > A' uk (d%) > [see Eq. ffT2Tj) ]. On 

the other hand, consistency with CKM unitarity would allow both A 12 fc(e^) and A' llfc (c?^) 
to be nonzero, but would imply a strong correlation on their magnitudes. RPV correc- 
tions of order ~ 0.1% are not unreasonable from the standpoint of expected magnitudes 
of superpartner masses or couplings. Indeed, having 

A^(4)~0.1 (154) 

implies that 

rrizk 



- 40A i?fc (155) 
100 GeV 3 V ; 

or Mgfc ~ 1 TeV for A^fc ~ V 'Anxa. As discussed in Ref. [94J, RPV corrections of this 



magnitude are not inconsistent with analogous bounds on RPV corrections obtained from 
other precise measurements, such as the studies of rare decays or flavor-changing neutral 
current processes. 

The implications for superpartner loop corrections for CKM unitarity have been studied 



in Ref. [98|. The contributions from W propagator Weu e vertex and external leg corrections 
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cancel from the difference [Afg — Af M ], leaving only the dependence on the W /iz/^ and Wdu 
vertex and external leg corrections as well as the box graphs entering the (3- and /x-decay 
amplitudes. As with the corrections to R e /n, these box graphs are subdominant in the 
presence of gaugino-higgsino mixing, as the vertex and and external leg corrections receive 
logarithmic enhancements. As a result, the dependence of the correction [Af^ — Af M ] on 
the parameters in £ so ft simplifies considerably, and the authors of Ref. [98] were able to 
perform an analytic study of the corresponding parameter dependence. 

To illustrate the dependence of [AfJ — Af M ] on the SUSY parameters, it is useful to 
consider several representative cases: 

(i) For situations in which the scalar superpartners of left- and right-handed fermions 
mix via the /i-term and tri-scalar SUSY-breaking interactions into mass eigenstates 
/i )2 , one has 

\Ar v - Af 1 ~ a{& ~ ^ In ^ + . . . h 56 ) 

^ 13 SUSY loop 32n 2 d 2 § 2 ' ^ } 

where the m ? are the corresponding physical masses and where this expression holds 

in the limit \mj 2 — mjj » m f 1 - Note that the overall sign of [Af J - Af „] SUSY loop 
depends on the relative degree of mass splitting among the sleptons and squarks. At 
the time when the analysis of Ref. [98] was carried out, the CKM unitarity deviation 
favored a negative sign, leading to the requirement 

/ \ 1/2 

m ih> fm^\ (157) 



For nearly degenerate squarks, one would need m^ 2 £ 3m^ 1 . For values of m^ close 
to the current direct search bounds , a m ass splitting of this magnitude is ruled out 



by the recent muon (g — 2) M results [125 



(ii) The constraints from (g — 2) /t may be evaded by taking M\ R w for the scalar 
fermions 25 . In this case, only the superparters of the left-handed fermions contribute 
to the CC weak interaction, and the requirements on the scalar fermion spectrum 
are rather different. To illustrate, we consider the limit of large sfermion masses, 
yielding the asymptotic expression 



[AfJ-AfJ SUSY _ loop ~^cos2/3 



1 Ml . m\ Ml rn~ 



ln^T-^Mn 



3 m (M|) ml (M|) 



+ ■■• , (158) 



where (MJ) 1 / 2 is the mass scale associated with the chareinos and neutralinos. For 
tan p > 1 as currently favored, one requires m~ > 3m? in order to obtain a negative 



25 



From Eq. (|157p we observe that the degree of mixing among squarks must always be less than that for 
smuons. Hence, setting M\ R = for the smuons implies a similar condition for the squarks. 
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sign for [Ar^ — Af M ]susY-ioo P (up to small logarithmic corrections). Note that such 
a sfermion spectrum would conflict with models that assume a universal sfermion 
mass at high scales, since in this case SU(3)c contributions to the renormalization 
group evolution of the first generation squark masses increases their magnitude at 
the weak scale relative to that of the first and second generation slepton masses. 

(iii) A final possibility arises when there exists significant mixing among either the u or 
o?-type squarks, but in a way that is not identical for both. In this case, gluino loop 
effects dominate the SUSY corrections to the Wud vertex, and one can accommodate 
a negative sign for [AfY — Af^Jgusy-ioop without requiring significant mixing among 
the smuons. In order to evade the (g — 2)^ constraints, the dominant contribution 
to M\ R for the d squarks must arise from a large value for the triscalar coupling 
Ad- In this case, avoiding color or charge-breaking minima in the scalar potential 
implies that all by the CP-even Higgs h° must be quite heavy. One can suppress these 
supersymmetric SU(3)c loop corrections for M 3 > 500 GeV. In this case, however, 
one returns effectively to the situation characterized by item (i), leaving item (ii) as 
the only viable option. 

At present, the CKM first row unitarity situation is unclear, so one cannot draw strong 
conclusions about either PRV or R-parity conserving scenarios. However, the prospective 
implications for the superpartner spectrum or presence of RPV interactions underlines the 
importance of resolving the various theoretical and experimental uncertainties germane to 
CKM unitarity outlined above. 



V. NEUTRAL CURRENT EXPERIMENTS 



A. Introduction 



Historically, the study of parity-violating (PV) neutral current interactions has played 
an important role in elucidating the structure of the electroweak interaction. In the 1970's, 
PV deep inelastic scattering (DIS) measurements performed at SLA C confirmed the SM 
prediction for the structure of weak neutral current interactions 126| . These results were 



consistent with a value for the weak mixing angle given by sin dw ~ 1/4, implying a 
tiny U(electron)x A(quark) neutral current interaction. Subsequent PV measurements - 
performed at both very low scales using atoms as well as at the Z-pole in e + e~ annihilation 
- have been remarkably consistent with the results of the SLAC DIS measurement. 

The value of s 2 at scale Q = ^/\q 2 \ = Mz has been determined precisely via Z-pole 
precision measurements at LEP and SLD. The global fit to all the precision observables 



gives 
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s 2 (M z ) = 0.23118 ± 0.00017 (159) 

The SM predicted scale-dependence of s 2 , on the other hand, has never been established 
experimentally to a high precision. The solid curve in Fig. [H] (see Section IIIip shows the 
running of weak mixing angle sin 2 9yy [in the modified minimal subtraction (MS) scheme] 
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as a function of the scale Q. The dip around M w is due to the decoupling of W boson 
when Q < My/. When we approach lower and lower energy, quarks with mass m q > Q 
decouple and the slope of the s 2 running changes. For Q < 1 GeV, all the heavy quarks 
decouple and s 2 is roughly a constant. The difference between s 2 at Z-pole and s 2 at 
low energy is: s 2 (0) - s 2 (M z ) = 0.00749 ± 0.00015 ± 0.00007, where the first error is the 
experimental error 26 s 2 (Mz) and the second is the theoretical error associated with the 
running to Q — 0. 

Recently, determinations of s 2 at various low energy scales have been performed, al- 
though the experimental error is still relatively large. The cesium atomic parity-violation 
(APV) experiment measured the parit y- forbidden transition between atomic states by ex- 
ploiting Stark interference effects Yli\ . The cesium weak charge, which depends on s 2 at 



Q ~ 0, appears to be consistent with the SM prediction. At h ighe r energies, the NuTeV 



collaboration measured v- (z/-) nucleus deep inelastic scattering [1281 ] . and studied the ratio 
between the cross section of the neutral current and that of the charged current. The result 
can be interpreted as a determination of s 2 at Q ~ 3 GeV. The observed deviation of cross 
section ratios from the SM predictions implies a +3<r deviation in s 2 at that scale. 

More recently, the SLAC E15 8 collaboration measured the electron weak charge via the 
PV M0ller (ee) scattering (l29) . It determined s 2 at Q 2 « 0.026 GeV 2 and obtained a 
result that is consistent with the SM prediction at the 1.1 a level. The Qweak experiment, 
which plans to measure the proton weak charge via elastic PV ep scattering using polar- 
ized electron beam at the Jefferson Laboratory (JLab), is currently under construction. 
The expected 4% measurement of proton weak char ges a t Q 2 ~ 0.03GeV 2 corresponds to 
0.3% determination of sin 6 W : 5sm 2 e w = 0.0007 poL better than any of the current 
measurements. 

Several future PV experiments are being considered. One proposal involves a more 
precise version of PV M0ller (ee) scattering at JLab using the planned 12 GeV upgrade 
of the accelerator. With a 2.5% precision in electron weak charge m easu rement at Q 2 ~ 
0.008 GeV 2 , s 2 could be determined with an error of 5s 2 = 0.00025 |l3l| . comparable to 



the Z-pole measurements. Another proposal is to study the PV electron-Deuterium DIS 
at JLab with current 6 GeV electron beam, or at future 12 GeV upgrade. Although the 
expected precision is worse than the Qweak experiments, such measurements have unique 
opportunity to probe the V(e) xA(q) interactions, QCD higher twist effects, possible charge 
symmetry violation in the parton distribution functions, and the ratio d(x)/u(x) at x — > 1. 
An DIS electron-Deuterium experiment has been approved for running with the present 



6 GeV beam [1321 ]. a nd s everal options for experiments using the future 12 GeV beam are 



under consideration 133 



In Table [V] we list the precision of various current measurements and possible future 
experiments, along with the sensitivity of the measurements of s 2 . Precise determinations 
of the value of s 2 at different scales, obtained from different types of experiments, will 
provide a consistency check of the SM at loop level. Any significant deviation from the SM 



26 Here, we have taken the value of s 2 (Mz) from a ht to precision data rather than the value computed 
from Eq. ([81]). 
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Measurements 


<5sin 2 9w/sin 2 


9w Ssin 2 9w 


Z-pole 


0.07% 


0.00017 


0.6% APV Q w (Cs) 


0.7% 


0.0016 


NuTeV iz-DIS 


0.7% 


0.0016 


13.1% SLAC E158 Q w (e) 


0.5% 


0.0013 


*2.5% JLab M0ller Q w (e) 


0.1% 


0.00025 


4% JLab Qweak Qw(p) 


0.3% 


0.00072 


*0.8% JLab eD DIS-parity 


0.45% 


0.0011 



TABLE V: Precision of various experiments which are sensitive to the value of sin 2 9\y at low 
energies. Entries with * are ideas for future experiments. 



prediction would constitute striking evidence for new physics. These high precision, low 
energy measurements will be sensitive to the new physics up to TeV scale. For example, we 
can write down the e ffective four fermion operators that contribute to the parity-violation 



eq scattering as [134 



C% = CZ + CZ = -^TfrW E ^wffQ + J2 K?fq, (160) 

q i 

where the first term is the SM contribution, and the second term gives new physics effects, g 
is the typical new physics coupling, A is the new physics scale, and h v is an 0(1) coefficient 
that parameterizes the new physics contributions for different quarks. A 4% measurement 
of proton weak charge Q v w corresponds to a probe of the new physics scale of 



• w 

A 



4.6 TeV. (16L 



9 JV2GJ6QP 



in 



Table [VTl summarizes the s ensit ivity to various new physics scale from current and future 
parity- violation experiments |l3l| . Also shown are the direct collider search limits for new 
physics from current colliders (LEP2, CDF and HERA) and indirect search limit from 
current electroweak precision fit. The scales that the low energy precision measurements 
can probe is close to - and in some cases even exceeds - the ones accessible in high energy, 
direct searches. Even after LHC begins running, the low energy measurements will be 
able to probe new physics scale comparable to the LHC reach. Once any new physics is 
discovered at LHC, low energy precision measurements can be used to probe details of the 
new physics, such as the couplings and charges of new particles. 

In the following sections, we discuss in detail three different types of neutral current 
experiments. Sec JVBI is devoted to parity violating electron scattering (PVES), which 
includes ee M0ller scattering, ep elastic scattering, and eD deep inelastic scattering. Atomic 
parity violation is discussed in Sec. IV D\ and neutrino-nucleus DIS is discussed in Sec. IV El 
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Z' 


models 


leptoquark 


compositeness 




m(Z x 


) m{Z LR ) 


m LQ (up) m LQ (down) 


e — q e — e 


Current direct search limits 


0.69 


0.63 


0.3 


0.3 




Current electroweak fit 


0.78 


0.86 


1.5 


1.5 


11-26 8-10 


0.6% Q w (Cs) 


1.2 


1.3 


4.0 


3.8 


28 


13.1% Q w (e) 


0.66 


0.34 






13 


* 2.5% Q w (e) 


1.5 


0.77 






29 


4% Qw(p) 


0.95 


0.45 


3.1 


4.3 


28 



TABLE VI: Sensitivity of new physics scale from various current and future low energy precision 
measurements. Entry with * is an idea for future experiment. Also shown are direct search limits 
from current colliders (LEP, CDF and Hera) and indirect search limit from current electroweak 
precision fit. The various new physics scale presented here are mass of Z' with extra U(l) [m(Zx)], 
or in left-right models [m(Zi,ii)], mass of leptoquark in up quark sector [uilq (up)], or down quark 
sector [mLQ (down)], composite scale for e — q composite, or e — e co mpos ite. Entries with 



either do not exist or do not apply. This Table is updated from Ref. 13 1| . 



B. Parity Violating Electron Scattering: M0ller and Qweak 



The SLAC E158 experiment measured PV ee M0ller scattering at Q 2 ~ 0.026GeV 2 (l29| . 



whil e the Qweak experiment at JLab will measure PV ep scattering at Q 2 ~ 0.03GeV 2 



1301 ] . In both cases, polarized electron beams are used to measure the PV asymmetry 



N R + N L ' y ' 

where Nr (Nl) is the number of detected events for incident electrons with positive (nega- 
tive) helicity. Although the dominant contribution to the scattering is via parity conserving 
photon exchange, the interference between the photon exchange and the parity violating Z 
exchange generates the PV asymmetry and filters out the much larger electromagnetic scat- 
tering effects. At leading order in Q 2 , the contributions to Apy are governed by A(e) xV(f) 
operator, with the coefficient being Q^, the "weak charge" of the target fermion, / (/ = e 
for ee scattering, and / = p for ep scattering, where Q v w = 2Q^ + Q%/). 

G 

£eff = -T^Qwn^efjJ . (163) 

At tree-level in the SM the weak charges of both the electron and the proton are suppressed: 
Qw = —Qw = 1 — 4sin 2 Qw w 0.1. One-loop SM electroweak radiati ve correct ions further 
reduce this small number, leading to the predictions Qw = —0.0449 135 . 136| and Q v w — 



0.0716 [i~36l . The factor of ^10 suppression of these couplings in the SM renders them more 
transparent to the possible effects of new physics. Consequently, experimental precision of 
order a few percent, rather than a few tenths of a percent, is needed to probe SUSY loop 
corrections. The advantages of measuring sin 2 6 W in these two PVES experiments are that 
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hydrogen target is used in both experiments, which is a relatively clean environment since 
the QCD contamination of heavy nuclei can be avoid. It is also theoretic ally clean since 
the hadronic uncertainties are relatively small and under control |l35l . 136l . 



The ee M0ller scattering experiment measured a parity violating asymmetry 129j |: 
Alr = (—131 ±14 (stat.)±10 (syst.)) x 10 -9 , leading to the determination of the weak mix- 
ing angle sin 2 9 W = 0.2397 ± 0.0010 (stat.) ± 0.0008 (syst.), evaluated at Q 2 = 0.026GeV 2 . 
Comparing to the SM prediction of sin 2 9 W = 0.2381 ± 0.0006 at this energy scale, the 
El 58 results agreed with the SM value at 1.1 a . When expressed in terms of the electron 
weak change, the difference between the measured value and the theoretical predicted one 
is 5Q e w = (Q e w ) exp - (Q e w ) su = 0.0064 ± 0.0051. 

In the Qweak experiment, the SM prediction for the left-right asymmetry is 290 x 10~ 9 . 
The expected experimental precision for ep scattering at Qweak experiment is 4%. When 
translated to a precision in the sin 2 9w measurement, it corresponds to Ss 2 of 0.0007, which 
is two times smaller than the precision in Cs APV and in NuTeV measurement. 

In order to interpret these high precision electron scattering experiments in terms of 
possible new physics, it is crucial to have the hadronic uncertainties under control. There 
are two kinds of QCD uncertainties that one must consider: (a) QCD corrections to the 
weak charge itself, and (b) QCD effects that impact the extraction of the weak charge from 
the experimentally measured PV asymmetry. The QCD corrections to weak charges have 
been discussed in Sec. IHII In the case of Qweak, the latter set of uncertainties is brought 
under control by a judicious choice of kinematics and by extrapolating the results of other 



experimental measurements. At forward angle 9, the PV asymmetry can be written as [137 



A PV = ^jS^[Q p w + F(9,Q 2 )} } (164) 
4v 



ira 



where F{9, Q 2 ) is the unknown form factor, that is proportional to Q 2 at low energy. The 
form factor F(9, Q 2 ) depends on a linear combination of isovector and isoscalar electromag- 
netic (EM) form factors as well as those associated explicitly with strange quarks. While 
the results of parity conserving electron scattering experiments provide the needed isovec- 
tor and isocalar EM contributions to F(9,Q 2 ), the strange quark contributions can only 
be determined by additional PV electron scattering measurements. These contributions to 
F(9, Q 2 ) - which vanish with Q 2 as Q 2 — » - have bee n stu died wit h an extensive pro- 
gram of measurements by the SAMPLE [l38|, HAPPEX [H^, PV A4 [l4~oj . and GO [l41 



Collaborations. The results yield tight constraints on the strange quark contributions to 
F(9,Q 2 ), which one can then extrapolate to the much smaller Q 2 relevant to the Qweak 
experiments. On the other hand, one cannot take Q 2 to be too small since Apy itself is 
proportional to Q 2 and the statistical error increases for smaller Q 2 . Experimentally, these 
two effects are optimized and Q 2 is chosen to be about 0.03 GeV 2 . The resulting hadronic 
uncertainty from the form factor is about 2%, half of the total experimental uncertainties. 

The precise measurements of the weak charges could probe both supersymmetric loop 
effects as well as tree-level, RPV contributions. For the .R-parity conserving MSSM, contri- 
butions to the weak charge and Q v w appear at loop level. With higher-order corrections 
included, the weak charge of a fermion / can be written in terms of the parameters Ptvc(O), 
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k, and \ vA introduced in Section HTT1 

Q f w = Pnc(0) [2l( - 4Q f k(0,p)s 2 {p)\ + \ f v + fa* - 4Q /S 2 ) \{ + box , (165) 

where 4 an d Qf are, respectively, the weak isospin and the electric charge of the fermion 
/, and s 2 = siia 2 9w(Mz). The quantities pivc(O) and k(0, p) are universal in that they do 
not depend on the fermion / under consideration. Detailed expressions for p and k can be 
found in Sec. IIHI The corrections X VA , on the other hand, depend on the fermion species. 
They include the vertex and external leg corrections to the weak charge. The relev ant 
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expressions for supersymmetric contributions to the X VA can be found in Ref. 142 

Note that at tree-level, one has p^c — 1 = k and \ VA = 0, while both SM and new 
physics contributions enter at loop level. Consequently, in what follows we will refer to 
5p = Pnc — 1 an d 5k = k — 1. 

The MSSM loop cont ributions to the electron and proton weak charges have been an- 
alyzed in detail in [1421 ] . which is reproduced in Fig. [9j Here, we plot the MSSM loop 
contributions to the shift in the weak charge of the proton, 5Qyy = 25Qw + 5Qy V , versus 
the corresponding shift in the electron's weak charge, 5Qyy, normalized to the respective 
SM values. The dots show the results of a random scan over a range of MSSM parameters. 
The loop corrections in the i?-parity conserving MSSM can be as large as ~ 4% (Qyy) an d 
~ 8% (Qw) ~~ roughly the size of the experimental errors for the two PVES measurements. 
Given the current results on E158, the SUSY loop contributions is consistent with the mea- 
surement at about 2 a level. In general, the MSSM effects are larger for large tan /3, light 
SUSY particles, and large splitting between sfermions, although the latter is an isospin 
breaking effect and therefore constrained by oblique T parameter. 

The shifts SQ^ are dominated by Sk SVSY since the corrections to Qyy due to shifts 
in the ppy parameter are suppressed by 1 — 4s 2 . In addition, the non-universal correc- 
tions involving vertex corrections and wavefunction renormalization experience significant 
cancellations. Since £k susy is universal, and thus identical for both and Q^y, the 
5k dominance produces a linear correlation between the two weak charges. The correc- 
tion to 5k is nearly always negative, corresponding to a reduction in the value of 
sin 2 0w £ (Q 2 ) — k(Q 2 , p)s 2 (p) for the parity-violating electron scattering experiments 
[see Eq. (1T6%]) ]. 

As evident from Fig. [91 the relative sign of the corrections to both Q v w and Qyy - 
normalized to the corresponding SM values - is nearly always the same and nearly always 
positive. Since Q v w > (Qw < 0) in the SM, SUSY loop corrections give 5Q V W > 
(5Qyy < 0). This correlation is significant, since the effects of other new physics scenarios 
can display different signatures. The combined measurements of the electron and p roton 



weak charge can be a probe to distinguish other new physics, as illustrated in Fig. [TD] [143 . 



In these studies, the electron anapole contribution was included in the parameter k and we denoted 
k pv ; the quantity pnc{®) was denoted ppy; and the sum of the terms containing the \ f V A in Eq. (IBB1 
were denoted by A/, which did not include the electron anapole moment contribution. 
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FIG. 9: R elativ e shifts in electron and proton weak charges due to SUSY effects (updated plot 



from Ref. 142j |). Dots indicate MSSM loop corrections for ~ 3000 randomly-generated SUSY- 



breaking parameters. Interior of truncated elliptical regions (a) and (b) give possible shifts due 
to R-parity non-conserving SUSY interactions (95% confidence), using value of 5|U u d| 2 /|U u d| 2 (a) 
and (b) in Table. IIIIl respectively. 



The arrows indicate correlated effects. For example, while both the superpartner loops and 
leptoquark exchange give a positive contribution to the proton we ak charge, only MSSM 
give rise to a sizable effect on the electron weak charge [134 . 1361. for the general class 



of Z' theories based on Eq gauge group, with neutral gauge bosons having mass ^ 1000 
GeV, the effect s on Q v w and Qyy also correlate, but 5Qw '/Qw can have either sign in this 



1341 In the case when Eq Z' models and MSSM have similar effects on the 



case 

electron and proton weak charge, measurement of the cesium weak charge using atomic 
parity violation can further tell these two apart, as explained in Sec. IV Dl below. 

If we relax the assumption of .R-parity conservation, tree level corrections to the weak 
charges are generated RPV interactions. In this manner one obtains the following effective 
four-fermion Lagrangian, using the notation of A^ k as defined in Eq. f|48l) : 

£rpv = -A'iki(q k L)d R j' M d R e L ^e L + A' llk (d k R )u L ^u L e L -f^e L 

-A m (e^) [p^Ll^y-LeLl^eL + h.c] , (166) 

where we have taken \q 2 \ w?~ and have retained only the terms relevant for the PVES 
scattering. The last term contributes to the muon decay, which affects the extraction of 
the fermi constant from the muon decay lifetime. Note the absence from Eq. (11661) of 
the parity-violating contact four-electron interaction. This is because the superpotential in 
Eq. (fT2|) can only produce parity-conserving contact interactions between identical leptons. 
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FIG. 10: Comparison of anticipated errors for and Q w with deviation from the S M expected 
from various extensions and allowed range (at 95% CL) by fits to existing data 1431 ] . 



The relative shifts in the weak charges are [144] : 



5Q 



e 
W 



Q 



5Q P 



1 -4s 



2 ' ^ 



Ai 2fc (4) = -29.8A 12fc (4 



W 



w 



-2A,A m (et) + 2A' m (4) - A' lfcl (qJ) - A m (e* 



1 -4S 2 , 

18.7A m (4) + 55.9A' 1U .(4) - 27.9A' lfcl (g*) 



' ~fe 



s 2 (l - s 2 ) 



1 



l-2s 2 l-Ar SM 



0.35 . 



(167) 



Since the A^ fc are non-negative, Eq. (I167j) indicates that the relative shift in is negative 



semidefinite. On the other hand, the relative shift in Q v w can have either sign depending 
on the relative magnitudes of A 12 fc, A' llfc , and A' lfel . 

T he qu antities A^, etc. in Eq. (I167P are constrained from the existing precision 
data |l44| . A summary of the existing constraints is given in T able II III in Sec. [Til which in- 



cludes superallowed nuclear /5-decay that constrains \V u d\ [1451 ] . at omic PV measurements 



of the cesium weak charge |l46j |. the ratio R e /^ of tti 2 decays |l47 |. and a com paris on 



of the Fermi constant G„ with the appropriate combination of a, Mz, and sin 9w [148 



The 95% CL region allowed by this fit in the &Q v w jQ v w vs. SQw/Qw pl &ne is shown 
by the closed curves (a) and (b) in Fig. [93 corresponding to the RPV fit with the value of 
^iKtdlVlKdl 2 of case (a) and (b) in Table. IIII[ respectively. Note that the truncation of the 
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initially elliptical curves is due to the sign requirements Aij k (f), AJ- fc (/) > [see Eq. (1351)]. 
The correction to Q w and Q v w from RPV SUSY could be two to three times larger than 
the SUSY loop effects. In addition, the prospective effects of Pr non- conservation are 
quite distinct from SUSY loops. The value of SQ W /Q W * s never positive in contrast to the 
situation for SUSY loop effects, whereas SQw/Qw can have either sign. Thus, a comparison 
of results for the two parity-violating electron scattering experiments could help determine 
whether t his e xtension of the MSSM is to be favored over other new physics scenarios (see 
also Ref. [136]). If SUSY is the new physics beyond the SM, it is in particular important 
to know whether .R-parity is conserved or not. Indeed, if -R-parity is conserved, then the 
neutral LSP (for example, the lightest x°) would be a suitable candidate for dark matter. If 
any deviation of the electron and proton weak charge is observed, the correlation between 
these two would help us to identify whether or not there is .R-parity conservation and, thus, 
shed light on the feasibility of SUSY dark matter. 

Ideas for measuring sin 2 Ow at low energy with higher precision have been explored 
recently. Ther e is a proposal of a similar M0ller (ee) scattering measurement at JLab 12 
GeV upgrade fl3l| . The estimated precision for electron weak charge is 2.5%. It could 



be used to determine the value of sin 2 Ow at Q 2 ~ 0.008 GeV 2 , with a 0.1% precision: 
£sin 2 6w — 0.00025, which is com parable to the precision of sin 2 0w determined from Z- 



pole precision measurements [131] . Such high precision enables us to constrain the SUSY 



parameter space, whether or not a deviation of Q w is observed. 



C. Electron-Deuterium Parity Violating Deep Inelastic Scattering 

In light of the recent developments in parity-violating ee and ep scattering discussed in 
Sec. IV B\ a new generation of PV PIS measurements wit h de uterium targets has been con- 
sidered at JLab 6 GeV beam [l32 ] and 12 GeV upgrade jl33l |. The DIS-parity experiments 



seek to study the deep inelastic scattering of a longitudinally polarized electron beam on 
unpolarized deuterium target. Neglecting target mass and higher-twist corrections as well 
as contributions from sea quarks, the PV asymmetry for eD DIS has the simple form: 

„e£>, dis 3G^Q 2 2Ci u — C\d + Y(2C 2u — Cm) aQ -, 

A ' v = ^Ph a 5 ' (168) 

where 

Y = 1 T 2 (1 ~ 2 €m i ff v and R(^Q 2 ) = — «o.2, (169) 

1 + (1 - y) 2 - y 2 R/{l + R) a R 

and y G [0, 1] is the fractional energy transfer to the target in the lab frame. The quantities 
Ci q parameterize the low-energy, PV electron-quark interaction 

£ pv = % E V^rfl^riiA + C 2q e^eq W q] . (170) 
q 

Using the SM values for C iq at tree level, one obtains 



A e Jy DIS « 10~ 4 Q 2 



l(l + Y)-(f + 6Y)s^0 w 



(171) 
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FIG. 11: T he li mits on C2 U and C m lis ted by the particle data g: 
experiment 138 ], and by DIS-parity [133]. The plot is taken from Ref. 



2J, by the SAMPLE 

im. 



The DIS asymmetry is much larger than Apy in the M0ller scattering and Qweak exper- 
iments: for Cj 2 =3.7 GeV 2 , A e ^J DIS = 0.0003. An expected 0.8% measurement in A e Jy' DI 
corresponds to 0.45% precision in s 2 : 5s 2 = 0.0011. While the sensitivity to s 2 in eD DIS is 
not as good as in the SLAC E158 and Qweak experiments, it has the unique opportunity to 
constrain the combination of 2C2 U ~ CW Assuming the successful completion of the Qweak 
experiment, an absolute uncertainty of 8C lu r d \ = 0.005 will be obtained. With this prospec- 
tive limit, DIS-parity experiment places an absolute uncertainty of b{2 Ciu — Cia) = 0.026. 



When taken together with the results from the SAMPLE experime nt I l38l| . much tighter 



bounds are placed on C^u and than were previously available [27J, as illustrated in 

Fig, m nig. 

For the R-parity conserving MSSM, loop corrections to Ci q appear. The Ci q are conve- 
niently computed using the expressions 

C lq = 2p Arc J 3 e (/| - 2Q q ks 2 ) - h\ (172) 

C 2q = 2p NC Ii(II - 2Q e ks 2 ) - l -\ q 2 , (173) 

where the A* 2 contain the appropriate combinations of the \ v A that are process dependent. 
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Detailed expressions for expressions for the \\ can be found in Ref. (TBI . 

As before, tree level contributions to Ci q arise for RPV SUSY. In terms of the Aijk(f) 
and AJ . fc (/), one has the following shifts in the Ci q 



RPV 



AC- 



Ac*r 



where X x is defined in Eq. ( I167p . 



\Ciu - -X x ]A 12k (e R ) - A' uk (d k R ), 

{C ld +^X x ]A 12k (e k R ) + A[ kl (q k ), 

\C 2u — 2A x ]A 12 fc(e^.) 
[C 2d + 2X x ]A 12k (e k R ) 



A' llfc (4), 



A' 



1 feiCoi) ? 



(174) 

(175) 
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FIG. 12: 95 % CL allowed region for RPV contribution to A e p y' D1S (y = l,Q 2 = 3.7 GeV 2 ) 
vs. electron weak charge (a) and proton weak charge (b). The dots indicate the SUSY loop 
corrections. The figures are reprinted from Ref. |l5Cl| ] with permission from Elsevier. 



we illustrate the sensitivity of A e Jy' DIS to the effects of MSSM loop contri- 
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We plot the relative shifts in A e p y vs. 



those 



In Fig. 

butions and tree-level RPV effects 
in Q e w and Q v w . The interior of the truncated ellipse gives the 95% CL. region from RPV 
effects allowed by other precision electroweak data. A deviation of about 1% could be ex- 
pected from MSSM loop effects while the maximum correction from R PV e ffects would be 



— 1.5%, corresponding to about 2a for the precision proposed in Ref. 1331 ] . The presence 
of RPV effects would induce negative relative shifts in both A e p ^j DIS and Q^, whereas 
the relative sign of the loop corrections is positive in both cases. A sizable positive shift 
in Q v w (up to 3cr for the proposed Qweak measurement) due to RPV contributions could 
correspond to a tiny effect on Ap 1 ^ DIS whereas a substantial negative shift in the proton 
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weak charge could also occur in tandem with a substantial negative correction to A e pV ' 
On the other hand, even a result for Q p w consistent with the SM would not rule out a 
sizable effect on A e J^ DIS . The addition of an eD DIS measurement would provide a useful 
complement to the PV ee and elastic ep measurements, assuming it can be performed with 
~ 0.5% precision or better. 



D. Atomic Parity Violation 



The effects of atomic parity violation (APV) can be measured either by observing the 
rotation of the polarization plane of linearly polarized light, or by measuring the rate for 
a Stark induced transition. The most precise measurement of an APV effect has been 
perfor med b y the Boulder group using the Stark interference method on a beam of Cesium 



atoms [127]. The PV transition between two atomic states can be expressed as 



M(n'P 1/2 — > nSi/2 



G/j, 



C SP (Z)Q W (Z,N) + - 



(178) 



where the dots are th e eff e cts f rom finite nuclear size, nucleon substructure, and nuclear 
spin-dependent term |l36 . 151 ] . The uncertainty associated with these effects is small, 



about 0.15%. Qw(Z, N) is the weak charge for atom, which is a combination of the weak 
charge of the up and down quark: 

Q W {Z, N) = (2Z + N)Q U W + {Z + 2N)Q d w « Z{\ - 4sin 2 B w ) - N « —N (179) 



The coefficient Csp{Z) parametrizes the contribution related to the atomic structure |l52l . 
A precise measurement on the cesium tran sit dipole amplitude can be used to determine 



Csp with relatively small uncertainty [153 



The cesium weak charge extracted from the APV measurement is [127 



Q^ s (exp.) = -72.69 ± 0.48, 



(180) 



with the combined experimental a nd t heoretical uncertainty to be about 0.6%. This is 
consistent with the SM prediction jl36l 



Q^ S (SM) 



-73.16 ±0.13. 



:isr 



The correction to the atom weak charge can be written as a sum of the weak charge of 
the up and down type quarks: 



6Q W (Z, N) = (2Z + N)5Q U W + (2N + Z)6Q d w . 



:i82) 



Using the results of the MSSM corrections to the up and down quark weak charge as 
described in Sec. IV B] we can obtain the contribution of SUSY loop corrections to the weak 
charge of heavy nuclei probed with APV. Since the sign of SQ^ / Q^ w due to superpartner 
loops is nearly always the same, and since > and < in the SM, a strong 
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cancellation between SQ^ and 5Qyy occurs in heavy nuclei. This cancellation implies that 
the magnitude of superpartner loop contributions to 5Qw(Z, N)/Q W (Z, N) is generally less 
than about 0.2% for cesium and is equally likely to have either sign. Since the presently 
quoted uncertainty for the cesium nuclear weak charge is about 0.6% |l54j - cesium APV 
does not substantially constrain the SUSY parameter space. Equally as important, the 
present agreement of Qyfi with the SM prediction does not preclude significant shifts in Q e ^ 
arising from SUSY. The situation is rather different, for example, in the Eq Z' scenario, 
where sizable shifts in would also imply observable deviations of from the SM 
prediction. 

There are several ongoing atomic parity violation experiments fol lowin g the cesium 
APV results. A more precise cesium APV measurement is underway |l55 |. The Seattle 
group is measuring the APV effects using Ba + ions, which is able to achieve the same level 
of precision as for cesium [l56 |. The measurement of APV effects along an isotope chain 



would elimi nate the large theory uncertainties from atomic structure. Efforts are underway 
at Berkeley [l57 ] to measure APV along Yb isotope chain. 

A measurement of the Helium weak charge via + — > + transition is under investi- 



gation at JLab, which could be used to cross check the Cesium APV experiments |158 . 
Furthermore, weak charge data from JLab on both Hydrogen and Helium would have the 
advantage of correlated errors which would cancel in the ratio Q v w jQ^-, potentially yield- 
ing significantly tighter constraints on new physics. A preliminary study showed that a 
1% measurement of the Helium weak charge c an constr ain RPV SUSY at nearly the same 
level as a 4% proton weak charge experiment |l58l . 159l . 



E. Neutrino-Nuclei Deep Inelastic Scattering 



Neutrino scattering experiments have played a key role in elucidating the structure of 
the SM. Recently, the NuTeV collaboration has performed a precise determination of the 
ratio R v (Rp ) of neutral current and charged current deep-inelastic (z/ M )-nucleus cross 
sections [l28 |. which can be expressed as: 



R 



v(v) 



a{v{v)N -> vX) 
a{y(v)N -»■ H+)X) 



;i83) 



where r = <Jp^ /a^ and (gf R ) are effective hadronic couplings (defined below). Compar- 
ing the SM predictions [27J for (g^) 2 with the values obtained by the NuTeV Collaboration 
yields deviations: 5R y (v) 



R cxp _ R SM 
U u(u) U u(P) 



5R U = -0.0033 ± 0.0007, 5R P = -0.0019 ± 0.0016. 



;i84) 



Within the SM, these results may be interpreted as a test of the scale-dependence of 
the s 2 since the (gf^) 2 depend on the weak mixing angle. The results from the NuTeV 
measurement imply a +3<7 deviation at Q ~ 3 GeV. This interpretation of the NuTeV 
results has been the s ubject of considerable debate. Unaccounte d fo r effe cts, such as 
NLO QCD corrections [l60L 161 ], Electroweak radiative corrections 161 . 162| . strange sea 
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asymmetry [1631 ] . isospin violation |164j, nuclear shadowing [1651 ]. o ther nuclear effects such 
as neutron excess in the ta rget and nuclear cross section effects 1661 ]. electron neutrino 



content in the NuTeV beam 167l ] have been proposed as poss i 
Alternatively, one may consider physics beyond the SM [161 



)le remed ies for the anomaly. 



168, 1691. In this review, we 



16* 



focus on the MSSM effects on the neutrino-nuclei scattering 

For momentum transfers satisfying |g 2 | <^ Mf , the neutrino-quark interactions can 
be represented with sufficient accuracy by an effective four fermion Lagrangian: 



r NC 



G^pNC 



C 



cc 

vq 



^f%7 A (l - 75>„ E ^[4(1 " 75) + 4(1 + J 5 ))q , (185) 
^ 7 A (1 - 75 )^7a(1 - 75 )d + h.c. , (186) 



V2 



where 



-QqK p 



;i87) 
:i88) 



The parameters p cc and X q L R are similar to t hose quantities that are defined in Sec. IIII1 
The relevant expressions can be found in Ref. 1681 ] 28 . 

The NC to CC cross section ratios R u and R p can be expressed in terms of the above 
parameters via the effective couplings 
way: 

,,\ 2 / N2 



v 

f R ) 2 appearing in Eq. (1 1 8 3 [) in a straightforward 



(9 



L.R) 



Ml 



M 2 W 



( pNC \ 
\PccJ 



E 



~L,R) 



(189) 



The SM values for these quantities are [27| (gf) 2 = 0.3042 and (gf ) 2 = 0.0301 while the 
NuTeV results imply (gf) 2 = 0.3005 ± 0.0014 and (gf) 2 = 0.0310 ± 0.0011. 

The MSSM l oop contributions to R u and Rp are highly correlated, in the range of about 
— 1.5 x 10 -3 |l68| . The sign of the SUSY loop corrections is nearly always positive, in 
contrast to the sign of the NuTeV anomaly. There is one corner of the parameter space 
which admits a negative loop contribution. This scenario involves gluino loops, whose 
effect can become large and negative when the first generation up-type squark and down- 
type squarks are nearly degenerate and left-right mixing is close to maximal. Although 
the gluino contribution could be as much as fewxlO -3 in magnitude, equal and large left- 
right mixing for both up- and down-type squarks is inconsistent with th e oth er precision 
electroweak inputs, such as the M w and charged current universality [l70 |. and color 
neutrality of the vacuum (unless the H°, A , and H ± become super heavy). In addition, 
the negative contribution from gluino to R u and R 9 could not account for the apparent 
deviation of sin 2 9yy from the SM prediction implied by the NuTeV analysis when a Paschos- 
Wolfenstein type relation is used. 



28 



In that reference a subscript u vN" on pnc.cc was included. 
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When .R-parity is not conserved, one obtains the effective Lagrangian for neutrino-quark 
scattering via tree-level contribution: 

£rpv = -A' 2kl {d k L )d R ^d R 9^^ L + A' 21k (d k R )d L ^d L 9^^ L 

-A' 21fc (4) [UL-fdLPLj^ + h.c] . (190) 

The corresponding shifts in Rup\ are 

SRvp) = K[-~e u L + |ei][l + rH>]A m (4) - 2[J% + e d L ]A' 21k (d k R ) + 2r™ > e d R A' 2kl (d k L ) 
« -0.25[1 + r^)A 12h (e R ) - 2[J% - 0.43]A' m (4) + l.Br^A^^). (191) 

As we discuss in Section [TT1 Ai 2 fc(e^) and A 2lk (d R ) are constrained by other precision 
electroweak data, while A' 2kl (d k L ) is relatively unconstrained. In Eq. ( 11911) . the coefficients 
of A 21k (d R ) and A' 2kl (d k L ) are positive, while the coefficient of A 12 fc(e^) is negative. Since 
the Aij k are non- negative, we would require sizable value of A 12fe (e^) and rather small 
values of A 21k (d R ) and A' 2kl (d k L ) to account for the negative shifts in R v and Rp implied by 
the NuTeV result. The present constraints on Ai2k(e R ) from other precision electroweak 
observables, as listed in Table HTl| however, are fairly stringent. The possible effects on R v 
and Rp from RPV interactions are by and large positive. While s mall n egative corrections 



are also possible, they are numerically too small to be interesting [168 



In short, the MSSM - with or without R-parity conservation - is likely not responsible 
for the NuTeV anomaly. The culprit, apparently, is to be found elsewhere. 

Finally, we note that another proposal to measure the weak mixing angle at Q 2 = 
4x 10~ 6 GeV 2 with a reactor-based experiment via v e e~ elastic scattering has been proposed 
m The estimated error on s 2 is about 1%, comparable to APV and NuTeV results, 

but with substantially different systematic contributions. Such measurement could explore 
the electroweak corrections, though the possible implications for SUSY have not been 
analyzed. 



VI. FLAVOR, CP, NEUTRINOS, AND COSMOLOGY 

The issues of flavor and CP symmetries are generally challenging for SUSY phenomenol- 
ogy. In the case of the MSSM with R-parity conservation, the structure of the soft SUSY- 
breaking Langrangian allows for a variety of flavor changing neutral current (FCNC) pro- 
cesses that must be suppressed in order to be consistent with experiment. By itself, the 
general structure of the soft Lagrangian does not provide for this suppression, so models 
of SUSY-breaking mediation must be constructed that provide it in a natural way. Simi- 
larly, jC so ft contains a host of new CP-violating phases beyond the phase of the Standard 
Model CKM phase that accounts for CP-violation (CPV) in the neutral kaon and B-meson 
systems. If these phases are 0(1) and if the soft masses are on the order of a TeV, the 
associated CPV interactions can give rise to permanent electric dipole moments (EDMs) 
of the electron, neutron, and neutral atoms that are up to two orders of magnitude larger 
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than experimental EDM limits. Short of any fortuitous cancellations between various CPV 
effects, there is no a priori reason to expect large suppressions of these phases as needed 
for consistency with experiment. The corresponding "SUSY CP problem" again provides 
a challenge to model builders. 

Both the SUSY flavor and CP problems have been reviewed extensi vely elsewhere , and 
we r efer the reader to excellent recent discussions (see, e.g., Refs. (iR lira flTi [rTi [l75l 



176| ). Here, we focus on aspects of these issues most relevant to the current experimental 



efforts in the low-energy sector as well as on recent theoretical developments pertaining 
to their broader implications for particle physics and cosmology. After reviewing general 
features of SUSY flavor physics and CPV, we concentrate on three areas of interest: (a) 
lepton flavor violation and the corresponding implications for neutrino physics; (b) EDM 
searches and their theoretical interpretation; and (c) implications for SUSY baryogenesis 
and dark matter. 



A. General Considerations 

Flavor 

Within the Standard Model, the GIM mechanism provides an elegant explanation for 
the suppression of FCNCs among quarks. In the limit of degenerate quarks, a sum over all 
intermediate quark states in loop contributions to FCNC processes yields a vanishing result 
due to the unitarity of the CKM matrix. The natural scale for FCNC effects - such as 
K°-K° mixing and b — > 37 - is thus governed by differences in the squares of quark Yukawa 
couplings (and products of off diagonal elements of the CKM matrix). The presence of 
these factors provides a natural way to understand the observed suppression of FCNCs. 

In general, superpartner loop contributions to FCNCs can upset the GIM suppression 
mechanism. For example, the difference of scalar quark masses need not be small compared 
to the weak scale, so that the corresponding loop contributions can be enhanced relative 
to those arising from quark loops. Similarly, flavor mixing among squarks need not be 
suppressed since there exists no a priori reason to expect flavor non-diagonal terms in the 
squark mass matrix to be small compared to the diagonal terms. Thus, studies of FCNC 
semileptonic or hadronic weak interactions can provide important constraints on the flavor 
structure of £ so ft- 

It has become conventional to characterize these SUSY flavor-violating effects by as- 
suming that they are small enough to be described by single insertions of the relevant 
flavor-violating soft m ass p arameter. In this "mass insertion" approximation, one may 

m 



consider the parameter [17 



^ab) 13 = - ABhj _ 1/2 (192) 

{M\ A ) it {M* BB )„ 

where A, B denote L or R and i and j are flavor indices. Note that flavor can be violated 
separately among the left- and right-handed fermion superpartners as well as in mass terms 
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that mix them after electroweak symmetry breaking, viz 



(Mini 



Vd 

V2 



fi taa/3 + {a d ) i 



(193) 



for down-type squarks. While the term containing the Yukawa matrix can be diagonalized 
by performing the same rotation on L- and i?-squarks that diagonalize the quark mass 
matrix, the term containing the triscalar coupling (ddjij will generally remain flavor non- 
diagonal after this rotation. 

Although the mass insertion approximation may not always accurately reflect the scale 
of flavor- violation in a given scenario, it provides a useful framework for comparing the 
constraints on SUSY flavor violation obtained from different experi ments . A summary of 
present limits on the (Sab)^ for various processes can be found in Ref. |l72j |. As emphasized 
by the authors of that work, there does not exist a sufficient set of experimental observables 
to completely determine the flavor-violating parameters that enter the (Sab)^ even within 
the MSSM. Consequently, one must use the experimental limits as input for model building. 
To this end, several broad approaches have been pursued. Among the most popular are: 



(i) Universality, which assumes that the soft terms are flavor diagonal and universal. 
For example, one may take 10 

M) = m 2 1 (194) 

for f = Q,U, D, L, E and 

a, = A f Y, (195) 



for f = U,D, E. Assuming that Eqs. (1194111951) hold at the SUSY-breaking scale, RG 
evolution to the electroweak scale will induce corrections to the relations fll94H195l) at 
the electroweak scale. To the extent that the RG evolution is dominated by Yukawa 
interactions, one would expect the soft parameters at th e electroweak scale to have 
an expansion in the Yukawa matrices(see, e.g., Ref. 1771 ] and references therein): 



Mi 



Ml 



M 



D 



m 2 



m 2 



m 2 



Sxl + h x Y\Y u + b 2 Y\Y d + b 3 ( Y\Y d Y{ Y u + Y{Y u Y d Y c 
a 2 l + b 4 Y u Y[ 
a 3 l + b 5 YdY\ 



(196) 



with similar expressions for the triscalar couplings Au t D to third order in the Yukawa 
matrices. Equations (11961) illustrate an alternative approach known as "minimal 
flavor violation" (MFV) in which all of the flavor violation in the soft sector is 
dictated solely by the structure of the Yukawa interactions. 

Alignment, a scenario in which the soft interactions can be diagonalized by the same 
rotations that diagonalize the SM Yukawa interactions. 
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Both the universality and alignment approaches build a "super GIM mechanism" into 
the soft SUSY-breaking Lagrangian and protect one against the appearance of large FCNC 
effects. Specific models for SUSY-breaking mediatio n m ay or may not lead to either univer- 



sality or alignment (for a discussion, see e.g., Ref. [1781 ]). and neither of these approaches 



may ultimately be correct. Nevertheless, they provide a useful starting point for the phe- 
nomenology of low-energy precision tests of SUSY in the flavor sector. 

Lepton Flavor and Number 

In lepton sector, total lepton number (LN) is an exact, accidental symmetry of the SM, 
while lepton flavor violation (LFV) involving charged leptons is highly suppressed by the 
scale of neutrino mass. Neither feature generally carries over to SUSY. In addition to the 
long-standing scrutiny of the flavor structure of £ so ft, there has been considerable recent 
interest in the possibility of total lepton number violation (LNV) in SUSY models. In the 
MSSM, LNV arises when the requirement of Pr conservation is relaxed, allowing for the 
existence of the A, A', and // terms in the superpotential of Eq. (fl2j) . As discussed earlier, 
such interactions may generate tree-level corrections to SM CC and NC interactions, some 
of which provide rather stringent constraints on the associated coupling to mass ratios. 
The interaction s in Wa r,= i may also give rise to a non-zero rate for neutrinoless double 
/3-decay (Oudd) [l 79 . 180| | - a AL = 2 process - as well as to LN conserving but LFV 



proc esses such as /i — > ej, /i — > e conversion, and /i — > 3e (for a recent discussion, see, e.g 
Ref. [l81| and references therein). For TeV scale superpartner masses, the interactions in 
W / al=i may have important consequences for the interpretation of these LNV and LFV 
processes as we discuss below. In addition, the existence of AL = 1 SUSY intera ction s 



imply the existence of a radiatively-induced Majorana mass term for the neutrino 182 
while extensions of the MSSM that explicitly allow for RH neutrino superfields have taken 
on renewed interest (the l iterature on th e topic is vast; for representive discussions, see 



e.g., Refs. (TH, H, H, H, E3). In SUSY models that explicitly contain tree- 



level Majorana masses, one may also find a possibility of relatively low-scale leptogenesis. 
Below, we review some of these recent developments involving LNV and the neutrino sector 
in SUSY. 



CP 

The plethora of new CPV phases that arise in £ so f t leads to similar complications for 
phenomenology as in the case of flavor. There simply do not exist a sufficient number of 
experimentally accessible CPV observables to independently constrain all of the phases, 
and one is generally forced to adopt simplifying, model assumptions. The situation is 
simplest in the Higgs and gauge sectors, where all CPV phases may be rotated into the 
following relat ive phases between the /i parameter and the gaugino mass parameter (see, 
e.g., Ref. [l89] ): 

Arg {fiMib*) Arg (i^M*) (197) 
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where i,j run over the three gauge groups of the MSSM (leading to a total of three 
independent phases in this sector). The analysis of CPV in this sector is often further 
simplified by assuming a common gaugino mass parameter at high scales, thereby reducing 
the number of independent phases to one. When discussing this special case, we refer to 
the common relative phase of the //-parameter and gaugino mass parameters as M . 

In the most general situation, the parameters in the scalar sector of £ so ft allow for an 
additional 37 independent CPV phases (for a discussion of parameter counting, see, e.g., 



Ref. [16j). As with the Higgsino-gaugino sector, the number of independent phases can be 
reduced by adopting a version of the flavor universality or MFV scenarios. In the latter 
case, for example, most of the CPV phases can be absorbed by sfermion field redefinitions, 
leaving only the parameters in Au as complex. Alternatively, assuming flavor diagonality 
for the scalar ma ss m atrices Mq UD but a general set of triscalar couplings, one has the 



additional phases 189l 



Arg (AfM*) Arg (A f A},) . (198) 

A third, commonly employed scenario is that of a common triscalar coupling and a single 
CPV phase for the gaugino masses at high scales, leaving only one additional phase <pA- in 
what follows, we will consider this minimal scenario wherein only two phases - M and 4>a 
- need be considered. 

The presence of non-zero RPV terms in the superpotential can introduce a large number 
of additional CPV phases. We will note review the implications of these additi onal sources 



of CPV in this review and refer the reader to recent literature (see, e.g., Ref. 190]). 

Connections with Cosmology 

While the search for CPV beyond that of the SM is interesting in its own right, there 
exists additional motivation from cosmological considerations. In particular, an explana- 
tion of the small (but anthropically relevant) baryonic component of the energy density of 
the universe points to the need for C PV beyond that of the SM. Indeed, as first observed 



by Sakharov nearly four decades ago [191| . arriving at a particle physics-based accounting 
for the baryon asymmetry of the universe (BAU) requires three ingredients in the particle 
physics of the early universe: (a) violation of baryon number (B); (b) violation of both C 
and CP symmetry; and (c) a departure from thermal equilibrium, assuming that CPT is an 
exact symmetry. In principle, SM contains all three ingredients. Baryon number violation 
arises through anomalous "sphaleron" processes that cause transitions between different 
electroweak vacua having different Chern-Simons number and, therefore, different total 
B + L. At temperatures above the electroweak scale, these transitions are mediated by 
the excitation of gauge field configurations called sphalerons. At lower temperatures, the 
probability of sphaleron excitations is Boltzmann suppressed, and transitions between dif- 
ferent vacua can only occur via exponentially suppressed tunneling. The SM also contains 
electroweak CP violation as well as C violating interactions (the gauge-boson couplings to 
axial vector currents). Finally, a departure from thermal equilibrium occurs as the universe 
cools through the electroweak temperature and the gauge symmetry of the SM is sponta- 
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neously broken. The strength of the CPV effects are strongly suppressed by the l ight quark 



Yukawa couplings and Jarlskog invariant associated with the CKM matrix [192|, 1 1931 . 1194 



while the LEP II lower bound on mass of the SM Higgs boson precludes a strongly first 
order elec troweak phase transition as needed to prevent washout of the BAU (see, e.g., 
Ref. [l95j l 



A variety of particle physics scenarios have been proposed that attempt to circumvent 
these SM shortcomings in explaining the BAU via baryogenesis at different cosmic epochs. 
At present, have no conclusive evidence favoring either early time/high scale scenarios 
such as leptogenesis, or relatively late time, electroweak scale baryogenesis. From the 
standpoint of phenomenology, consideration of the latter is particularly attractive, since a 
combination of CPV searches, precision electroweak measurements, and collider studies can 
highly constrain and possibly even rule out electroweak baryogenesis (EWB). In addition, 
it is interesting to consider the possibility that new physics at the electroweak scale may 
provide both an explanation of the BAU and viable candidate for cold dark matter (DM). 
In this case, DM considerations can provide additional constraints on EWB. 

The study of SUSY DM remains an active, on-going field that has been reviewed ex- 
tensively elsewhere (see, e.g., Refs. fl96l. fl97j]\ Below, we focus on the viability of SUSY 
EWB, taking into account recent field theoretical developments and the phenomenology 
of EDM searches, precision electroweak measurements, collider studies, and DM consid- 
erations. In this respect, we pro vide an up-date to the exte nsive reviews of baryogenesis 
provided by Trodden and Riotto |l98l | and Dine and Thomas [l 99 . 



B. Lepton Flavor and Number Violation 



The best-known probe of LFV is to sea rch for the SM-forbidden decay 
current best limit on the branching ratio is 



B 



7) 



< 1.2 x 10 



-li 



90%C.L. 



e^vv) 



ej. The 
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obtained by the MEGA collaboration. A similarly interesting bound on the r ate fo r /i — > e 



converstion in gold nuclei has been obtained by the SINDRUM collaboration [201 



T[fi- + A(N, Z) -> e~ + A(N, Z)\ 
r[/i- + A(N, Z) -> vA(Z -l,N + 1)] 



< 8 x 10 



-13 



90%C.L. 



(200) 



In addition, stringent limits have be en o btained for other LFV rates: 1.0 x 10 



-12 



for 



B 



202 



A new experi- 



4.3 x 10~ 12 for B^. c \20% and 4.6 x 10~ n for 5™ e |204 
ment is being performed by th e MEG collaboration at PSI that hopes to reach a sensitivity 
of ~ 5 x 10~ 14 for -B M ^e7[205], while until recently there had been an effort by the MECO 
collaboration to probe Bf^ e at the level of 5 x 10~ 17 using the AGS at Brookhaven. Al- 
though that experiment has now been derailed by the U.S. funding agencies, efforts are 
underway to pursue an experiment at an alternate site, possibly using a future muon 
strorage ring at JPARC in Japan. 
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The prospective implications of these experiments for supersymmetric lepton flavor 
structure has be en a n alyz e d by several authors (for a discussion within non-SUSY models, 



see, e.g., Refs. 206 . 207l . 208 ] and references therein). A comprehensive analysis in a 
minimally-extended MSSM that includes RH neutrino supermultiple ts an d neutrino mass 
generation via the see-saw mechanism has been carried out in Ref. jl88j |. While the RH 
neutrino sector decouples from low-energy LFV observables due to the large RH neutrino 
mass {Mr ~ 10 12 GeV), the LFV effects of this sector can be communicated to the charged 
slepton mass matrices and triscalar couplings in £ so f t through RG running from high scales. 
The authors obtain general expressions for the rates for fi — > and /i — > e conversion in 
terms of the slepton and light sneutrino rotation matrices, Z 1 ^ and Z T J , introduced earlier 
without relying on the mass insertion approximation. One has 
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(201) 



where A\ ' r are the dipole amplitudes appearing in the amplitude for /i — > 



(*). 
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eG, e a *m P - q) [(q 2 la - iq a ) (if Pr + A L X P z ) 
+im„a aP qP (A r P r + A%P L )} e(p) 



(202) 



where, in contrast to Ref. [l88 ] we have normalized the amplitudes to G^. 

Contributions to these amplitudes are generated by the diagrams of Fig. [131 The 
chargino loop contribution to A R , for example, gives 



A 




mi. 



\U kl \ 2 Zl>Z 2 r F^) 



(203) 



-Uki Vk2 



Fi(x 



jk) 



where Xjk = vr? x Jmi~., and the F\^{x) are loop functions defined in Ref. |188j ]. Analogous 

expressi ons f or the other chargino and neutralino loop contributions to the A^ 1 are given 
in Ref. |l88j |. Assuming no cancellations among the various contributions, one may use 
Eqs. ( I201H203]) and experimental limits on B^ ei to derive bounds on the LFV parameters 
Z X JZ 2 J* etc.. Note that these LFV couplings also arise in the non (V - A) x (V — A) SUSY 
box graph contributions to the muon decay parameter, g^ R as in Eq. fl95|) . The bounds 
on these quantities are sufficiently stringent that only the L-R mixing terms can make 
appreciable contributions to g RR . In the case of fi — > e conversion in nuclei, the "penguin" 
amplitudes proportional to A x ' contribute to the four fermion e/jqq conversion operators, 
as do LFV Z°-exchange amplitudes and box graphs. As we discuss below, the A±' R may 
give the dominant contributions in the presence of Pr non- conservation. 

An older anal ysis using the mass insertion approximation has been performed by the 
209. ecial case in which the amplitude is dominated by 



authors of Refs. 
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FIG. 13: Contributions to the LFV amplitudes entering /i — * cy 
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photino loops one has 

24a sin 2 9 W 
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(204) 



where is an average slepton mass, x = M?/mf, and the M^x) are defined in Ref. 12091. 
The limits on the parameters (5 e AB )2i vary with x. At x = 1.0 one has 

(205) 



\(6l L ) 21 \ <1.9xl0- 3 



|(^) 21 |<4.3xl0- 7 



for mi = 100 GeV. These authors did not analyze the /1 — > e conversion process, so no 
limits on the (^5)21 are available for this process. 

In some scenarios, the mass insertion approximation may not provide a realistic gui de to 



the magnitude of LFV effects. A well-known illustration occurs in SUSY GUT models [2 11 



212| . wherein contributions from the large, third generation up- type Yukawa coupling to 
the RG running of the sfermion and triscalar couplings from the Planck scale to the GUT 
scale lead to splittings among these parameters. After implementing symmetry-breaking 
at the GUT scale and evolving these parameters to the weak scale one has for an SU(5) 
GUT 
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where Iq and I' G are contributions arising from the Yukawa induced running from Mp 
to Mq. It is convenient to rotate the sfermion fields by the same transformations that 
diagonalize the Yukawa matrices, 



V e e R 



il = UlL 



so that the term in Mj R containing I' G becomes 



^I'aVlM^ 



+ h.c. 



(207) 



(20^ 



76 



where M e is the charged lepton mass matrix. For an SO(10) GUT, the LH slepton mass 
matrix picks up an additional contribution from running above Mq, while the induced 
triscalar contribution becomes 

- ^4r G V*eMeVie' L + h.c. (209) 

In general, the numerical impact on of these LFV couplings is not well reproduced 

using the mass approximation. Importantly, the expected magnitude of this LFV ob- 
servable is expected to be well within the reach of the MEG experiment for superpartner 
masses on the order of a few hundred GeV. For sufficiently large m e - fi , however, scales 
as 1/mf , so a null result for the MEG experiment would imply TeV-scale superpartner 
masses in this GUT scenario. 

The foregoing analyses assume conservation of Pr, so that the potentially significant 
LFV effects at the weak scale are not accompanied by corresponding LNV. If Pr is not 
conserved, however, then the terms in the superpotential Wal=i can lead to observable 
low-scale LNV as well as LFV. An analyses of this possibility of low- scale LFV that is 



accompanied by low-scale LNV has recently been carried out in Ref. |207l | using an ef fective 
operator approach, following on the earlier analysis of Raidal and S ant amaria [2 1 3| . The 
leading LFV and LNV operators have dimension six and nine, respectively, and appear in 
the effective Lagrangians valid below the electroweak scale 29 . 

£ lfv = £ %Of > + • • • (210) 

i 

^NV = E|^ (9) + --- ( 211 ) 

i 

where the + • • • indicate terms containing higher dimension operators and where 

Of L = e iL a fa ,ipe jL F> a ' + h.c. (212) 

O lL = iiL^L^kL^mL (213) 

= t lL Tit qiX q q (214) 

with Tg q is a shorthand for all possible 7-matrix insertions and where the corresponding 
operators with RH fields have not been explicitly included 30 . In the case of £lnv one has 
operators of the general form that will contribute to Qvf3f3 

O i9) = qT iq qT 2 q eT 3 e c + h.c. (215) 



29 Consequently, the operators do not generally respect the SU(2)lxU(1)y symmetry of the SM 

30 In the semileptonic operator of^ we have not included chirality labels on all the fields since si 
tensor interactions flip chirality while vector and axial vector interactions preserve it. 
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where a complete list of operators of this form can be found in Ref. 214|. All searches 
for * //3/?- decay involve + — > + transitions, and for these cases the dominate operator is 



1791. 1214 



O 



(9)- 



[qRT + q L q R r + q L + q L T + q R q L r + q R ] e(l + 7 5 )e c 



with coefficient in the case of Pr non- conservation given by 

87ra s |A' m | 2 
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(216) 
(217) 

indicate contri- 



where rrtq is an average squark mass, nig is the gluino mass, and the + ■ 
butions proportional to the semiweak coupling ai2 ■ 

The chiral structure of Ox} ++ implies that it can contribute to an effective inree operator 



that contains no derivatives [2 14J , th ereby leading to an enhanced, long-range pion-exchange 
contribution to the 0i//3/3-decay rate |l79j |. From experimental upper limi ts on this rate and 



taking into account this long-range contribution, the authors of Ref. 179| | obtained the 
constraint 



A' m < 2 x 10~ 
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(218) 



lOOGeW V 100 GeV/ 
The coefficients of the LFV operators can similarly be expressed in terms of the A and 
A' couplings: 
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(219) 



where the various combinations of the A and A' entering c a are given in Ref. |18l| . Limits 
on various combinations of the Pr- violating couplings are given in Tables II and III of that 
work. From the present limits on B^ ni one has for example 
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while from B^ e one obtains 
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FIG. 14: Contributions from RPV interactions to (a) Oz^/3/3-decay and (b) neutrino mass. Figure 

(a) gives representative contribution from semileptonic trilinear RPV interactions parameterized 
by the coupling A' m ; the "+ • • • indicate contributions involving neutralino exchange. Figure 

(b) shows analogous neutrino mass contribution from semileptonic trilinear RPV; the cross and 
shaded circle denote fermion mass and sfermion L-R mixing insertions, respectively. In both cases 
generation indices are suppressed. 



Analogous limits on other combination of the couplings can be found in Ref. [181 . 

The results in Eqs. (I218H222]) lead to several observations. First, for superpartner masses 
of C(100) GeV, the bounds on A' m obtained from OvfiP are stronger than those derived 
from the LFV observables, though the gap between them shrinks as the superpartner 
masses are increased. Second, the LFV /i — > e conversion limits are more stringent than 
those derived from /i — > ej, even though B^ e contains an extra factor of e 2 suppression 
compared to B^ n . The reason is that the penguin operators generated by the terms 
proportional to A 1 ' in Eq. (12021) contain large logarithmic enhancements while the dipole 
operators do not, and these large logs overcome the nominal 4ira suppression of the conver- 
sion process. In addition, tree-level exchange of superpartners can occur in the presence of 
Pr non-conservation, leading to the presence of four-fermion operators not directly related 
to the LFV electromagnetic amplitudes o f Eq . (12021) . This situation differs from that of 
the Pft-conserving GUT scenario of Refs. |m[ l2lj . wherein the magnetic amplitudes are 
dominant and the naive expectations for the relative magnitudes of conversion and \x — > &y 
branching ratios obtains. 

This observations can have important consequences for the interpretation of 0vf3f3 in 
terms of light Majorana neutrino exchange. Indeed, For < O(TeV) masses in £ S oft, the 
contribution from cO^ ++ /A 5 to the rate for 0v/3j3 can be comparable to the contribution 
from the exchange of a light Majorana neutrin o. D enoting the amplitudes for the former 
and latter as A H and A L , respectively, one has |207l 



A 
A, 



ii 



(224) 



where rapp is the effective mass associated with the exchange of the light Majorana neutrino 
and k ~ 50 MeV is its typical virtuality. For A ~ 1 TeV and mpp ~ 0.1 — 0.5 eV, the ratio 
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in Eq. (12241) can be order unity. In general, one would like to know whether or not there 
exist important heavy particle contributions when attempting to ex tract robust bounds on 
mpp from the 0vf3f3 rate limits. To this end, it was observed in Ref. 2071 ] that the presence 
of important heavy particle (superpartner), Pr- violating exchange contributions is also 
accompanied by logarithmically-enhanced contributions to the conversion rate. Should 
the results of future studies find that the relative sizes of the LFV branching ratios are in 
accord with the naive expectations, one would conclude that there are no large Pr- violating 
contributions to the 0i//3/3-decay process. A similar conclusion holds for other, non-SUSY 
scenarios, pointing to the usefulness of LFV studies as a diagnostic for the presence of low 
scale LNV. 



C. R Parity Violation and Neutrino Mass 



One may incorporate massive neutrinos in SUSY models by supersymmetrizing any of 
the S M ex tensions that incorporate non-vanishing m v . As discussed above, the authors of 
Ref. 1881 ] carried out such an analysis in an extension of the MSSM that contains right- 
handed neutrino superfields iVj 7 that are singlets under the MSSM gauge groups. The 
small masses of the light neutrinos arise via the see saw mechanism with a right-handed 
neutrino mass parameter Mr ~ 10 12 GeV. Except for the effects of the resulting non-zero 
light neutrino Majorana mass, the effects of the fields in decouple from low-energy 
observables. Variations on t his c l ass of SUSY mode l s wit h massive neutrinos can be found 
in the literature (see Refs. pHl [H3, OH, EH, and references therein) . 

An alternate mechanism for generating neutrino mass involves the LN V in teractions in 
W&l=i (for a recent analysis and survey of the literature, see e.g., Ref. [HEl). Tree-level 
Majorana masses can arise from the bilinear terms that lead to mixing between neutrinos 
and neutral Higgsinosjsee Eq. ({12]) ]. Specifically, one has the contribution to the light 
neutrino mass matrix [2 151 ] 



m 



(225) 



where m is a characteristic soft mass parameter. In the absence of fine-tuning between 
the contribution in Eq. (I225j) and other tree-level neutrino mass terms, one obtains the 
following rough upper bound on the scale of the P^-violating bilinear coupling 
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3 x 10~ 6 

cos (3 
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100GeV\ 
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(226) 



The neutrino mass implications for the magnitude of the dimensionful Pr violating SUSY 
couplings are thus quite severe. 

The triscalar couplings in Wal=i may also induce contributions to m v through one-loop 
radiative corrections. Such contributions have been studied extensively, and we give only a 
flavor of these analyses here (for an extensive list of the literature, see Ref. [4] of Ref. |215| ). 
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Roughly speaking, one finds 

,(AA) 



l — 2 ^k\ 3kl m Ll (5 L LR ) kk g fc (227) 
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where 



{Ml)L + (M 2 R ) kk 



is a number typically of 0(1) and where we have assumed that (M|)£ fc ~ (M|,){ fc >> 
{Ml R ){ k . From Eqs. (I227p we may derive neutrino mass naturalness bounds on products 
of the triscalar couplings for a given value of the flavor diagonal LR mixing parameters, 
{8 LR ) kk . Assuming that (M| H )^ oc rrif, so that {5^ LR ) kk ~ rrif/m, we obtain the most 
restrictive bounds for third generation fermions (£ — 3): 



^-4X10-^)^) (229) 



^sS4xlO-»^j^ i ^j (230) 

Note that these limits are comparable to those obtained from LFV and LNV observables 
for rh ~ 100 GeV, though the dependence on the soft masses differ in the various cases. 



D. EDM Searches: Implications for SUSY 



Recent advances in experimental techniques have put the field of E DM s earches o n th e 
brink of a revolutio n. T he present limits on the EDMs of the electron jUJ, neutron [2 171. 



and mercury atom 218| - shown in Table IVHI - are already remarkably stringent. The 
pursuit of EDMs tha t began w ith the pioneering studies of the neutron by Purcell and 
Ramsey in the 1950's |219lj220j is poised to enter a new era with prospective improvements 
in sensitivity of up to four orders of magnitude. New e fforts are underwa y that ai m to 
push the sensitivity of the EDM searches for the electron |22l|, H, H], neutron 1 2 25 



2261, and neutral atoms [2271 . 12281 . |229| ]. as well as for the muon and deuteron |230j ]. These 
prospective experimental advances, as well as relat ed t h eore t ical issues and developments, 
have been discussed in seve ral r ecent reviews (Bl, Il89 . 231 . 232j | as well as the text by 



Lamoreaux and Khriplovich[233|, and we refer the reader to those publications for an 
extensive survey of the literature. 

While the potential improvements shown in Table IVHI would not provide access to 
EDMs associated with the CPV phase of the CKM matrix 234, 235 . 236| . they could allow 
one to observe an EDM associated with SUSY CPV. As noted above, present EDM limits 
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TABLE VII: Present and prospective EDM limits. Expectations based on SM (CKM) CP violation 
are also shown. 



System 


Present Limit (e-cm) 


Group 


Future Sensitivity 


Standard Model (CKM) 


e~ 


1.6 x 1(T 27 (90% CL) 


Berkeley 




< io- 38 


e 




Vein fPKO^I 

i aie \r dkj ) 


~ 1U 




e 




Indiana/ Yale 


^ m-30 




e 




Amherst 


~ m-30 




e~ 




Sussex (YbF) 


~ 10~ 29 






y.o x iu ^yu/o 






i n— 36 








in-24 




n 


2.9 x 10~ 26 (90% CL) 


ILL 


1.5 x 10~ 26 


1.4 x 10~ 33 - 1.6 x 10~ 31 


n 




PSI 


7 x 10~ 28 




n 




SNS 


2 x 10~ 28 




n 




ILL 


2 x 10~ 28 




199 Hg 


2.1 x 10~ 27 (95% CL) 


Seattle 


5 x 10~ 28 


< io- 33 


225 Ra 




Argonne 


1Q -28 




129 Xe 




Princeton 


io- 31 


< io- 34 


D 




BNL 


~ io- 27 




223 Rn 




TRIUMF 


~ io- 28 





already preclude 0(1) CPV phases and TeV scale superpartner masses in the MSSM, and 
the future measurements will make these constraints even more stringent. 

The theoretical interpretation of the present and prospective EDM searches in terms 
of the parameters of £ so f t requires a careful delineation of a variety of effects. The most 
straightforward analysis occurs at the level of operators involving the SM fermion and 
gauge boson fields. In the strong sector of the SM, the lowest- dimension, gauge invariant 
operator that can generate an EDM is the QCD #-term: 
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where is the SU(3)c field strength tensor and = (l/2)e a/3/iy G Q , /3 . The pre sent 



limits from the EDM of the neutron, d n , lead to the most stringent bounds on 9 [217 



|0| < (1.2 ± 0.6) x 1(T 10 (90% C.L.) 



(232) 



where the (1.2 ± 0.6) prefactor is obtained using the QCD sum rule result of Ref. |l89l ] and 
includes the theoretical error quoted in that work. 

CP-violation in the electroweak sector arises from the complex phase in the CKM ma- 
trix that enters the renormalizable interactions of quarks with gauge bosons. The 
magnitude of its effects is governed by the Jarlskog invariant (l7| . 

J = cos 9\ cos 92 cos sin sin 9i sin 9-x sin 5 = (2.88 ±0.33) x 10" 5 , (233) 

with the 9i and 8 being the three angles and complex phase in the CK M matrix. The 



EDMs ge nerated by SM electroweak CP violation arise at multi-loop level[ 234j . 12351 . 1237 



239l | and are proportional to J, suppressing their effects to the levels indicated in 
VIII (for a discussion of the corresponding SM effects in atoms and nuclei, see, e.g., 



ilo, EH, H, H, H, HH, EH, E3, B}). 



2381 
Table 
Refs. 

The effects of supersymmetric CPV arise through loop corrections to operators involving 
SM fields. For example, one-loop SUSY corrections to quark propagators can generate 
complex phases in the quark mass matrix, and through redefinitions of the quark fields, 
these phases can be absorbed into 9. In the absence of a Peccei-Quinn (PQ) symmetry 
that allows one to absorb these contributions into the axion field and maintain a vanishing 
9 prior to symmetry-breaking, the experimental limits on 9 given in Eq. f 1 2 3 2 j) lead to 



tight constraints on CPV in the SU(3) C sector of the MSSM[189j. Generally speaking, 



however, phenomenological analyses of supersymmetric CPV implicitly assume such a PQ 
mechanism, leading one to consider SUSY contributions to higher dimension operators. 
The lowest dimension non-renormalizable, gauge-invariant CPV operators arise at di- 
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mension six: 
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where the + ■ ■ • indicate gauge invariant operators involving lepton fields. Here, we have 
chosen to normalize the operators in terms of a new physics scale A taken to be greater 
than the scale of electroweak symmetry breaking. After electroweak symmetry-breaking, 
the terms containing the SU(2)l and U(l)y field strength tensors give rise to the electric 
dipole moments of the elementary fermions: 
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(235) 

(236) 
(237) 
(238) 



where cw = cos^vk and sw = sin^. 

The terms coupling G^ v to quarks are the chromoelectric dipole moment operators; the 
term containing three powers of G^ u is the CPV Weinberg three gluon operator; and the 
four fermion operator involves products of SU(2)x doub let and singlet fields (with SU(2)^ 
indices i, j) of flavors a, . . . , d. As discussed in Ref. [2491], the terms containing products of 
HlH u and GG etc., become topological when the Higgs fields are replaced by their vevs, 
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(a) (b) 

FIG. 15: Representative one loop super symmetric contributions to elementary fermion (a) electric 
dipole moment and (b) chromoelectric dipole moment (quarks only). 



and these effects cannot be analyzed in perturbation theory. They do, however, contribute 
to the operator C?^ and amount to a shift in the value of 6 that is constrained by d n . The 

parts of the HlH u GG operators containing physical scalar degrees of freedom contribute to 
a renormalization of £S? V and, thus, to EDMs. Naive dimensional analysis suggests that 
values of the operator coefficients w u ~ a s /4ir would lead to EDMs that are consistent with 
present limits for the mass scale A > 1 TeV. A comprehensive study of these operators, 
however, has not been performed at present. 

One-loop SU SY contributio ns to the operator coefficients in Eq. H234I) have been com- 



puted in Refs. [236 . 250 . 25 1| (for older one-loop computations, see the literature cited 
in these studies), while two- loop con t ributions to th e elementary fermion EDMs dj have 



recently been analyzed in Refs. HH, H, H, 2551. Illustrative contributions are show 
in Figs. [15] and [T71 Approximate results for the one-loop quark and lepton EDMs and 



quark chromo-EDMs have been given in Ref. (Tgl for a simplified scenario in which only 
two CP-violating phases contribute: </> M , a common relative phase between the electroweak 
gaugino masses and the \i parameter, and 4>Ai a common phase associated with the triscalar 
couplings: 



d 2 9 2 i ■ , fs& 9 2 A . , 

= sin 04 + 1 sm <p u tan a 

en e 12 y V 24 24/ M 

d 1 2a 2 

-2- = -f (sin^tan/^ 1 + sin0 A ) + 0(g 2 2 , fi) 

d G ^f7 3 

= -f (sin^[tan/3] ±1 + sin0 A ) + 0{g\,gl 

Kq lo 



where e q is the quark charge, 



K = m f 
16ir 2 m 



(239) 
(240) 
(241) 

(242) 
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and where for simplicity we take A = m, a common soft mass, and where the upper (lower) 
sign corresponds to negatively (positively) charged quarks. These results, together with 
two- loop contributions to the Weinberg three gluon operator coefficient w, can be used 
to compute the EDMs of charged leptons, the neutron, and neutral atoms. In the case 
of hadrons and atoms, one must contend with a variety of theoretical issues involving 
non-perturbative QCD, nuclear structure, and atomic s t ructure t heory. Extensive, recent 
reviews of these issues can be found in Refs. jl, Il89l . 231 . 232 . 233| . so we do not reproduce 



those discussions here. Instead, we illustrate the sensitivity of these EDMs to the CPV 
phases. Considering first the electron, we scale the SUSY mass scale m to 100 GeV and 
use g\ = 4na tan / cos 2 9w and g\ = Ana/ sin 2 9w to obtain 

-f « 5 x 10~ 25 ( -z— J [tan (3 sin M - 0.05 sin (j) A \ e - cm . (243) 



A \ m 



In this case, we see that for tan/3 £ 1, dj is overwhelmingly sensitive to the relative phase 
of the /i parameter and electroweak gaugino soft masses and that the present experimental 
bounds in Table IVHI imply 

/ ~ \ 2 

sin ^ 3xl0 "(j^v) coU3 • < 244 > 

In short, only for m > a few TeV can one accomodate an (9(1) phase and remain consistent 
with present limits. 

In Fig. [16l we illustrate the complementary sensitivity of various EDMs to the SUSY 
phases using illustrative results for the electron and neutron. Here, we have assumed a 
common sfermion mass mj — 1 TeV for the first two generations and show results for 
two different scenarios for the values of and M 2 . The widths of the bands correspond 
solely to the experimental error and contains no theoretical uncertainty associated with 
the non-perturbative methods used to obtain d n . The shaded band indicates the regions 
required to produce the baryon asymmetry during the supersymmetric electroweak phase 
trans ition . The band associated with the 199 Hg EDM limit is similar to that for d n (see, 
e.g., |l89j ]). The bands in Fig. [16] illustrate the general feature that various EDMs display 
complementary dependences on SUSY CP-violating phases and that results from a variety 
of EDM searches are needed to obtain meaningful constraints in a given scenario. Moreover, 
the scale of the allowed phases is quite small: ~ few x 10 -2 . This small scale does not 
appear to be a priori n atur al, and leads to the SUSY CP problem indicated earlier. 

As discussed in Ref. [1891 ] two-loop "Barr-Zee" contributions to the EDMs of elementary 
fermions may become important for TeV scale sfermions in the limit of large tan (3. For 
the electron, these contributions are given approximately by 

« ^ In tan/? sin(0 M + <f> A ) , (245) 

where Y t is the top Yukawa coupling and itia is the mass of the CP-odd Higgs. The de- 
pendence on these parameters arises from the coupling of the CP-odd Higgs to scalar top 
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FIG. 16: One loop constraints on CPV phases (horizontal axis) and (j>A (vertical axis) from 
present 95% C.L. limits on the EDMs of the electron (solid lines), neutron (dashed lines) and the 
baryon asymmetry (colored bands). Constraints consistent with the WMAP value for the baryon 
asymmetry are given by the blue band and those obtained from BBN are given by the blue 
+ green bands. Left panel corresponds to choosing (\/i\,M2) = (250,200) GeV (non-resonant 
baryogenesis) and the right panel corresponds to (200,200) GeV (resonant baryogenesis). In 
obtaining both figures, a common sfermion mass of 1 TeV was used. These figures are courtesy 
of C. Lee. 

quarks in the two-loop amplitude. Comparing with the dominant one-loop contribution 
we observe that - apart from the logarithmic factor - the two-loop Barr-Zee contribution 
is suppressed relative to the one- loop EDM by a little more than the 1/167T 2 loop factor. 
Since the dependence on tan/3 is similar in both cases, the two-loop Barr-Zee contribu- 



tion will be comparable with the one-loop contribution only if 



<< 



*)a\ and a large 



splitting between the masses of the scalar fermions and CP-odd Higgs generate logarithmic 
enhancements of the two-loop amplitude. 

The tight one-loop limits on the SUSY CP violating phases may be relaxed for suffi- 
ciently heavy sfermion masses, as in "split supersymmetry" scenarios. In this case, the 
dominant contributions arise from the two-loop graphs of Fig. [TTJ. Omitting the Barr-Zee 
contribution and letting 



d}(2 loop) = dj( 7 /i) + d}(Zh) + d}(WW) 



(246) 



corresponding to the three different graphs of Fig. [T71 one has, for example, [2521 . 12531 . 1255 



A 2 



^ . m(A ?) m '^ 



2 e 2 



M w m 2 h0 



2 7 '- //('•///) 



(247) 



87 



where rf H = (m +/ m /i°) 2 and Djj involves combinations of the chargino diagonalization 
matrices and Higgs-Higgsino-Gaugino couplings. In the simplified scenario discussed above 
in which the electroweak gaugino mass parameters have a common phase, Im(D^) oc sin tp^. 
Analogous expressions can be obtained for the other contributions in Fig. [T71 However, 
the l itera ture does not agree on the results for these graphs. Numerically, the authors of 
Ref. 252| obtain in the heavy chargino limit 



dJ(Zh) « OMdJ^h) (248) 

(T e {WW) « -0.3^(7/1) (249) 

dl(Zh) » cT^h) (250) 

dl(WW) » -0.7^(7^) (251) 

where the results for <i n include QCD evolution from the electroweak scale to the hadronic 
scale. Note that in the case of the electron, the Zh contribution is suppressed by the small 
vector coupling of the Z-boson to the electron. Moreover, the authors find substantial 
cance llati ons between the WW and 7/1 contributions to d n . In contrast, the authors of 



Ref. [2531 ] find an opposite sign for the WW contribution and no cancellation. Since 
the relative importance of this disagreement is less severe for d e we concentrate on the 
electron EDM below. The sensitivity of the two-loop EDM to the CPV phases is given 
approximately by 

«.o.m-2rK^ /100GeV\ 



A 



97 ( m^± \ ( 100 GeV\ / PN . , . , 

2 x 10- 27 -M — (fcnl® F,n{rt u ) e - cm . (252) 

V m h° J \ mh° J 



For F 1 n{j'f H ) ~ 0(1) and ImD^ ~ sin^^ the two-loop EDM is roughly 300 times less 
sensitive to sin M than the one-loop contribution. Thus, the present d e limits may accom- 
modate 0(1) phases for sfermion masses of O(10 TeV). For sin^ = 0.5, for example, the 
one-loop contributions become suppre ssed relative to the two-loop effects for mj > 4 — 5 
TeV and for 200 GeV < //, M 2 < 1 TeVH . We discuss the corresponding implications for 
SUSY baryogenesis and dark matter below. 

While the foregoing discussion relies largely on flavor-diagonal CPV, elementary fermion 
EDMs may also provide information on the flavor structure of £ so ft- For example, one-loop 
gluino contributions to the d-quark EDM can be enhanced by a factor of m^jm^ in the 
presence of flavor non-diagonal CPV |l89l 



J7 n ,d ( m b\ a s tan/3 

dl = eQ d ^ {- 2 ) (253) 

where 

^31 =Arg[fe) 13 fe) 33 (254) 

and where the {5^ B )^ have been normalized to an average sfermion mass-squared rather 
than to the denominator appearing in Eq. (11921) . The quantity 8f 31 also enters the imaginary 
part of the one-loop quark mass renormalization, and so will be constrained to the 10~ 9 
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FIG. 17: Two-loop contributions to elementary fermion EDMs. The figures are reprinted from 
Ref. H with permission from Elsevier. 



level in the absence of a PQ symmetry. Imposing the latter and taking the present EDM 
limits one finds bounds on the 5[ 31 ranging from 10~ 3 to 10 -6 for m = 1 TeV. 

Apart from the Weinberg operator proportional to TtGGG, the other dimension six 
operators appearing in Eq. (12341) have received less scrutiny than the the dipole operators. 
It has been observed in Ref. [2571 ] that the EDMs of the neutron and neutral atoms may 
receive important contributions from the four fermion operators in the limit of large tan /3, 
as these contributions grow with tan 3 /3. 



E. SUSY Baryogenesis and Dark Matter 

Explaining the origin of the matter density of the universe is an ongoing task that lies 
at the interface of particle physics, nuclear physics, and cosmology. As is well known, 
the smallest - but most anthropically relevant - component of the energy of the universe 
consists of baryonic matter, while the next largest component is the non-baryonic cold 
dark matter (DM). In principle, SUSY could provide a particle physics basis for explaining 
how both of these matter components came to be and take on their observed abundance. 
The literature pertaining to both supersymmetric dark matter and baryoKenesis is vast, 
and comprehensive reviews have appeared over the past decade [l~9ol fl97l uM fl99j . Here, 



we review recent theoretical developments and the corresponding implications for SUSY 
CP-violation. 

The baryon asymmetry of the universe (BAU) can be characterized by the ratio of the 
baryon density to the entropy density of photons: 



us = f (7.3 ± 2.5) x 10- 11 , BBN ^ 
s 1 (9.2 ± 1.1) x 1Q- 11 . WMAP [258 1 ' ' 



where the first value (BBN) is obtained from observed light element abundances and stan- 
dard Big Bang Nucleosynthesis and the second value is obtained from the cosmic microwave 
background as probed by the WMAP collaboration. 



89 



Despite the p resen ce of all three Sakharov ingredients discussed above, it was shown 
by Shaposhnikov [l92[ ] that they are not sufficiently effective within the SM to account for 
the observed BAU. In order to prevent "washout" of any baryon number created at the 
electroweak temperature, the electroweak phase transition (EPWT) has to be strongly first 
order. The strength of the EWPT depends on the parameters of the Higgs potential, Vh, 
that also govern the mass of the Higgs boson. In the SM, the mass of the Higgs must be 
below 45 GeV to allow for a strong first order EWPT, so the present LEP 2 direct search 
limit m h > 114.4 GeV precludes this possibility. Roughly speaking, this bound arises from 
a competition between thermal contributions to Vjj and those generated by parameters in 
the Lagrangian at T = 0: 

V H (T, 0) = Kfr(0) + V th (T, 0) (256) 
V cS (<p) = i m V0+^(0t 0)2 + ... (257) 

where V e g((f)) is the zero temperature one- loop effective potential and V t h(T, 0) gives ther- 
mal contributions to the potential. One may then use the high temperature e xpan sion to 



determine the shape of Vh{T, 0) as a function of T, leading to (see, e.g., Ref. [2591 ]) 



V H (T, 0) = D(T' 2 - T o 2 )0 2 - ET(f) 3 + ^0 4 (258) 

where To is the temperature below which the quadratic term becomes negative and where 
D, E, A(T) and To are functions of the scalar, vector boson, and fermion masses, v, and 
T. The EWPT is characterized by two additional temperatures in addition to T : 

9 8ADT 2 , . 

= SXD^E 2 (259) 
T 2 _ ADT 2 



c " XD-E 2 ' 

For T < T, a potential barrier forms between the phases of broken and unbroken elec- 
troweak symmetry, associated with minima of the potential at 0° = and 0° = v(T)/y/2 = 
m . The formation of a barrier is accompanied by the onset of tunneling between the two 
phases and the formation of bubbles of broken electroweak symmetry. For T < Tc, one 
has Vff(0 m , T) < Vh(4> = 0,T) while for T < T < Tc, the potential at = is a local 
maximum. 

In order to avoid washout of the baryon asymmetry during the EWPT, the energy 
associated with the sphaleron configurations, T sp h, must be sufficiently large as to suppress 
the associated transition rate between different vacua. The sphaleron energy is also tied 
to the value of v(T) since the gauge and Higgs fields are coupled through the associated 
equations of motion. From numerical studies, one obtains 

T c g T c 
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where B(X/g 2 ) is a constant of 0(1) that depends gently on X/g 2 and where the last 
inequality is imposed in order to ensure a value of Yb no smaller than given by Eq. (12551) . 
The latter requirement, thus, leads to 

^2 1 . (261) 
J-c 

The value of v(Tc) is, in turn, determined by minimization condition V'((j),Tc) = that 
yields 



v(T ) 2E A _ / v, 



T c A(T C ) \m 



AE Mr + • • • » ( 262 ) 



where the the "H — • " indicate small corrections arising from the logarithmic T-dependence 
of A. From Eqs. (I261II262 I ) one has 

AE Q) > 1 (263) 

as the condition on the parameters in Vh{4>,T) that must be satisfied in order to obtain a 
strong first order EWPT that prevents washout of the baryon asymmetry. 
In the SM, the cubic coupling is given by 

E = 2n?w + « 0.01 (264) 

implying that must be lighter than about 45 GeV in order to satisfy the condition (12631) . 
As this bound is well below the present LEP II lower bound, the electroweak baryogenesis 
in the SM is clearly ruled out experimentally. 

In the MSSM, the bound on the lightest CP-even Higgs mass can be relaxed either 
through additional scalar contributions to Vh{4>,T) that increase the value of E or choices 
of the other parameters that allow for a lighter Higgs boson (see the discussion in Section 
IVIIip . The former possibility is realized in the presence of a light, right-handed stop that 
couples strongly to the Higgs fields, contributes to E at the one-loop level, and enhances 
E by roughly an order of magnitude : 

yf sin 3 (\-A\lMl ) 3/2 
Emssm « E SM + \^ ~ 9 EsM ( 265 ) 

At = A t — n cot (3 . 

In principle, light left-handed stops could also generate large enhancements, but precision 
electroweak data rule out the existence of such a light II- The feasibility of the light tu- 
induced enhancement of E also depends on the avoidance of color- and charge-breaking 
minima, leading t o co nditions on the soft parameter M^ 3 . From the analysis of this re- 
quirement in Ref. 12601, one finds that RH stops lighter than about 130 GeV are disfavored. 



Assuming these conditions are satisfied, the enhanced value of Emssm increases the upper 
bound on the lightest Higgs from Eq. (12631) to ~ 120 GeV. 
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Going beyond the MSSM, one may st reng t hen the EWP T by introducing singlet Higgs 
supermultiplets S via the superpotential |261 . 262 . 263 . 264 



^ singlet = (ji + aS) H x Hz + 8S + ^S 3 



and the corresponding soft Lagrangian that contains triscalar couplings 
Aingiet = - {aA a H 1 eH 2 S + h.c.) - (^kA k S 3 + h.c. J + • • 



(266) 



(267) 



where e is the antisymmetric SU(2) tensor, and A a and A K are the soft SUSY breaking 
parameters. Writing the vevs of the neutral fields as 



Re(H?) = —= cos 7 cos 8 
V2 

Re (if 2) = -7= cos 7 sin /5 
V2 

the corresponding singlet contribution to the cubic term is given by 



(268) 



2v / 2sin 7 / A 2 n 2 A 2 

-C/smgiet = ^ I a^cos 7 sm 28 + -«A«sin 7 



(269) 



Super symmetric models with singlet Higgs are generically referred to as "next-to-minimal" 
and are motivated largely by a desire to generate the \i parameter as the vev of the singlet 
field 31 . Mod els of this type have received considerable attention recently, as in the studies 
of Refs. [2~6ol 



The feasibility of choosing the soft parameters to obtain a sufficiently larger value of 
-Ekngiet while respecting t he ex perimental lower bounds on the lightest CP-even Higgs were 
initially studied in Ref. |262j |. It was shown in that there exist considerable regions of 
the singlet parameter space leading to a sufficiently strong first order EWPT. The more 
stringent LEP II lower bounds on the mass of the Higgs reduce this available parameter 
space, but there remain considerable regions that admit v(T c )/T c > 1. 

In addition to these EWPT considerations, the size of the CP-violating asymmetries 
generated by particle physics interactions at the phase boundary must also be studied. In 
the SM, these asymmetries are highly suppressed by the Jarlskog invariant f)233p and by 



{vl - vlM - vlM - vlM - vDivl - vDivl - vl) « 4 x 10 



-17 



(270) 



31 In writing Eq. (|266p . we follow the conventions of Ref. [263(, where a possible quadratic term in eliminated 
by a constant shift of the field S — > S + c. In this case, an explicit jj, parameter appears. 
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since CP-violating effects shou l d van ish in the limit that any two quarks become degenerate. 



Farrar and Shaposhnikov [193l . |194| subsequently argued, however, that the relevant CP- 
violating asymmetry depends solely on the difference between the probabilities for reflection 
and transmission of s and (i-quark currents and the phase transition boundary, so that Yb 
is proportional to y s — yd rather than the combination in Eq. (1270)) . Assuming that a 
non-SM mechanism generates a strong first order EWPT, the resulting expectation for the 
BAU in the SM is much closer to the observed value than in Shaposhnikov's original work. 

As with the EWPT, the presence of supersymmetric interactions at the phase bound- 
ary may also lead to larger CP-violating asymmetries, as there exist a plethora of CPV 
interactions that are not Jarlsko g suppress ed. The computation of these asymmetries was 



initiated by the authors of Refs. [2661 . 12671 ] using conventiona l transpo rt methods and was 



followed up by a number of subsequent studies in the MSSM 2681 . 12691 ] . The results gener- 
ally indicate d th at SUSY CPV phases of £ 0(1) would be needed to obtain the observed 
values of Yr (268| . In the late 1990's, however, Riotto pointed out - using more sophisti- 
cated non-equilibrium field theory methods - that memory effects in the EWPT plasma 
could resonantly en hance the CP- violating sources needed for successful EWB with smaller 
CP-violating phases (270t | for appropriately tuned values of the MSSM paramete rs. Subse - 
quent detailed analyses of these resonant enhancements were performed in Refs. |27ll . 272l . 



These resonant enhancements allow for successful EWB in the MSSM with significantly 
smaller CPV phases than implied by earlier work, thereby allowing for consistency with 
the corresponding EDM bounds on these phases. 

The CPV sources appear in the transport equations for Higgs and quark supermultiplet 
densities that govern the projection of chiral charge at the phase boundary. In the case of 
the Higgs supermultiplet current density, H^, for example, one has 

= _ r " - r Y ( « - L + » ) _ f r ( * _ « + * ) + r b A + sf . (27!) 

k H \k Q k T k H J \k B k Q k H J k h H 



Q and (B,T) are the number densities of particles in the third generation left- and right- 
handed quark supermultiplets, respectively; the kjj,h,Q,T,B are statistical weig hts; Sf is 

a CP-violating source; and Th, IV, Ty, and r\ are transport coefficients. The terms 
proportional to Th and f\ cause any non-zero Higgs supermultiplet asymmetry to relax 
to zero, as favored by minimization of the free energy. The terms containing Ty and Ty 
favor the transfer of the Higgs asymmetry into the baryon sector and are, thus, essential 
for the generation of a non-vanishing Yg from CPV in the Higgs sector. The net baryon 
asymmetry depends on a detailed competition between the effects of the CPV sources and 
the CP-conserving relaxation and Higgs-to-baryon transfer rates. 

Analogous equations obtain for the quark supermultiplet densities. Carena et al ob- 
served that the enhancements of the sources occur when LH and RH scalar top 

H,T,Q 

quarks or Higgsinos and gauginos are near ly degenerate, leading to resonant scattering 
from the spacetime varying Higgs vevs [271 |. As noted above, however, the requirements 



of a strong first order EWPT and of precision electroweak data preclude the occurrence of 
such degeneracies in the scalar top sector, implying that resonant supersymmetric EWB 
may only occur via gauginos and Higgsinos. 
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274, 275 , who found 



276, 277 



These developments were followed up by studies in Refs. |273 
somewhat smaller resonance effects from the sources, and by the work of Refs 
in which the relaxation and Higgs-to-baryon transfer coefficients terms in the quantum 
transport equations were computed using the same non-equilibrium methods. In particu- 
lar, the authors of the latter work observed that the various terms in Eq. (I27ip (and the 
analogous quark supermultiplet equations) could be derived by expanding the Greens func- 
tions and self-energies entering the non-equilibtrium Schwinger-Dyson equations in ratios 
of physical scales present during the EWPT. It was also pointed that resonance effects 
could also enhance Higgsino and chiral charge relaxation coefficients (r#, etc.) , thereby 
mitigating the effect of enhanced CP violating sources. In addition, the lowest order contri- 
butions to the Higgs-Quark transfer coefficient, Ty, had been neglected in earlier analyses. 
The net impact of these refinements implies that the viability of supersymmetric EWB in 
the MSSM is a quantitative question, depending in detail on the parameters of the the- 
ory and requirement input from EDMs, present and future collider studies, and precision 
electroweak data. 

Although formal theoretical issues in EWB remain to be addressed, it is instructive to 
consider the phenomenological i mpli c ations of the r ecent studies for EDM searches. Recent 
analyses have appeared in Refs. |25Rl27i . l279ll280| . Illustrative results are given in Fig. [TH] 
, where we show the regions of the [i-M\ parameter space that correspond to resonant 
MSSM EWB (light blue bands), assuming GUT relation between the gaugino masses. The 
funnel-like structure corresponds to the resonance conditions: /i ~ M\ or /i ~ M2. The 
red region is excluded by LEP 2. The semicircular bands indicate the exclusion region 
derived from the present (dark blue) and prospective (blac k) el e ctro n EDM measurements 
for sin^ ~ 0.5 using the two-loop computations of Refs. |252 . 253 ]. An analogous set of 
plots for various values of tan/3 and rriA are given in Ref. (256j~ for AMSB type gaugino 
mass relations. 

The study of Ref. 256| indicates that the viability of resonant EWB in the MSSM 
implies \d e \> 10~ 28 e-cm and gaugino/Higgsino mass parameters < 1 TeV. Although some 
portion of the remaining parameter space can be explored with LHC studies, full collider 
exploration of the needed MSSM parameters will await the ILC. In addition, future dark 
matter detection experiments may provide a complementary probe, as the character of 
the LSP is governed by the same parameters \i and that determine the viability of 
resonant EWB. Assuming that the relevant portions of the MSSM parameter space leads 
to the observed DM relic abundance, one could expect resonant EWB to be accompanied 
by enhanced production of, and detection rates for, high energy neutrinos produced by 
neutralino annihilation in the sun. The absence of evidence for high energy solar neutrinos 
in Super Kamiokande data implies that a portion of the parameter space in Fig. [18] near 
the base of the EWB funnels would be excluded under such scenarios. Future neutrino 
telescopes and ton-sized direct detection experiments will considerably extend the reach of 
these DM probes of EWB. In most cases, however, obtaining the observed relic abundance 
requires modifications from standard cosmology, such as the presence of an additional 
energy-density that increases the Hubble parameter and allows for earlier decoupling of 
neutralinos. 
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FIG. 18: Constraints on MSSM parameters from resonant electroweak baryogenesis, two-loop 
electric dipole moment of the electron, and LEP II . This figure is courtesy of S. Profumo. 



VII. Z-POLE ELECTROWEAK PRECISION MEASUREMENTS 



The electron-positron colliders SLC and the LEP have produced millions of Z bosons 
at \fs ~ Mz- The measurements of Z lineshape, forward-backward asymmetries and 
polarized asymmetries at Z-pole lead to a precise determination of the Z-boson pole mass, 
the total Z-width and the Z couplings to fermions pairs. The Z-pole precision observables 
are also combined with the results from other experiments, like CDF, DO, and low energy 
atomic parity violation (APV) and scattering measurements to confront theories such as 
the Standard Model or other new physics extensions. The global fit to the electroweak 
precision observables are in excellent agreement with the Standard Model. Any new physics 
beyond the SM that could contribute to these precision observables is, therefore, tightly 
constrained by the existing data. In the i?-parity conserving MSSM, new corrections to 
the precision observables arise from superparticle loops. The SUSY contributions could be 
significant when the superparticles are light. Global anal yses of precision observables i n th e 
framework of MSSM have appeared in the literature 
which will be briefly reviewed in this section. 



281, 282, 283, 284, 285, 286, 287 



A. Precision Observables 



The Z pole observables can be organized into several groups |28al288j: 

• 9 lineshape observables 

A fit to the Z lineshape and the leptonic forward-backward asymmetries determines 
the Z boson mass Mz, the total Z- width Tz, the hadronic cross section a° ad , and, 
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for each lepton flavor, I — e, /i, r, the ratio Ri = r had /r H 32 , and the pole asymmetry 
A°/ B . Defining 

A = l-4g /S in 2 ^ ff 
f 1-4Q /S in 2 fl{ ff + 8Q 2 sin 4 # e y 

where sin 2 9^ s is the effective weak mixing angle for fermion / at the Z-pole, we have 

A% = \A e A f . (273) 

• 3 further LEP asymmetries 

These include the r polarization V{r) = A T , its forward-backward asymmetry 
Vfb(t) = A e , and the hadronic charge asymmetry, (Qfb), which is used as a deter- 
mination of sin 2 



• 6 heavy flavor observables 

These include the ratios R q = T qq /T ha , d , the forward-backward pole asymmetries A 
and the left-right forward-backward asymmetries A FB (q) = A g , each for q — b,c. 



0, g 
FBi 



• 3 further SLD asymmetries 

The precise measurements at SLD using the polarized electron beam determine A e , 
A^ and A T 33 . 

Besides the precision measurements performed at the Z-resonance, the W boson mass 
(M\y) and width (r^) have been measured to a relatively high precision at both the 
Tevatron and the LEP 2. These quantities are usually included in the list of precision 
observables. 

Low energy precision observables are sometimes included in global fits. These observ- 
ables include 

• Two weak charge measurements from atomic parity violation: Qj^ and Q^- 

• Deep inelastic scattering experiments that yield k, which is a linear combination of 
effective 4-Fermi operator coefficients. 

• u^e scattering experiments that yield the leptonic 4-Fermi operator coefficients gy 
and g v £. 

• NuTeV results on the neutrino-nucleus deep inelastic scattering. 

• The branching ratio of B — > X s ^. 

• Muon anomalous magnetic moment (g — 2) M . 



32 The results are combined into one Ri if lepton flavor universality is assumed. 

33 The results are combined into one Ai if lepton flavor universality is assumed. 
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B. SM Global Fit 



The Standard Model contributions to the precision observables depend on the free pa- 
rameters in the model. In the gauge and Higgs sectors one has five parameters: the 
three gauge couplings, the Higgs vev, and the physical Higgs boson mass (ra^). From 
among these, one ordinarily chooses the fine structure constant (a) and Fermi constant 
{Gfj, = l/\^2v 2 ) as independent inputs because they are known precisely from low-energy 
measurements. The remaining quantities in this sector that are relevant for Z-pole observ- 
ables are the mass of the Z-boson (Mz), TUh, and the SU(3)c coupling (a s ). In addition, Z- 
pole observables depend on the value of the running QED coupling a(Mz)- Since a(Mz)/a 
receives hadronic contributions associated with the five light quarks that are not calculable 
with the same precision as other SM contributions, one treats these contributions, denoted 

(5) 

as Ac^ a ; d , as an separate parameter to be obtained from the fits. The top quark does not 
contribute to the running couplings at \i = M z and below because m t > M z . However, 
one-loop radiative corrections depend strongly on m t so it is also treated as a fit parameter. 

The Z-boson mass M z has been determined at LEP to a high precision comparable 
to that of G^. Therefore, M z is sometime taken as fixed input instead of as a fitting 
parameter. The top quark mass mt has been measured directly at CDF and DO. Its 
value from the direct measurement is included in the global fit as a constraint. The strong 
coupling constant a s , can be dete rmin ed from from non-Z lineshape |289j ]; and Aa^ d , is 
obtained from hadronic r decay |2~90j . Both determinations are included in the global 
fitting as extra constraints. The lower bound on the Higgs mass from LEP Higgs searches 
is sometimes included as well. 

A global fit to the precision observables with extra constraints on m t , a s , Aaj^j and rrih 
determines the fitte d values of the input parameters. The studies from LEP electroweak 
working group [28a ] included the latest results on the mass of the top quark from Teva- 
tron: m t = 171.4 ± 2.1 GeV, the width of the W boson from the Tevatron and LEP-2: 
T w = 2.147 ±0.060 GeV, and the mass of the W from LEP-2: M w = 80.392 ±0.029 GeV. 
The program ZFITTER is used to calculate the SM predictions for those precision observ- 
ables, inclu ding full one-loop radiative corrections 34 and higher order QCD and electroweak 
corrections [2911 ] . The global fit to all the high Q 2 precision observables showed excellent 
agreement between the measurements and the SM fitted values. The discrepancies are 
usually less than 2 o for almost all the observables except for A b FB , where the deviation 
is about 3 a. The fitted Higgs mass is in the range of 85±28 GeV at 68% C.L. An upper 
limit of rrih < 166 GeV is obtained at 95% C.L. This limit increases to 199 GeV when the 
LEP-2 direct Higgs search limit of 114 GeV is included in the fit. 



ZFITTER uses on-shell renormalization. 
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C. MSSM Contributions to the Precision Measurements 



The success of the SM global fit to the electroweak precision measurements imposes 
strong constraints on any new physics extension beyond the Standard Model. Super- 
symmetric models can always avoid such constraints since supersymmetric corrections de- 
couple in the rh —>■ oo limit. Thus, supersymmetric models look just like the Standard 
Model if the supersymmetric mass scale is large enough. As long as the Standard Model 
with a light Higgs boson provides a good fit to the data, supersymmetric models can 
as well. On the other hand, large contributions from SUSY are possible when a super- 
symmetric spectrum includes light superparticles. There have been numerous attempts 
to identify constraints on the MSSM param eters from precision observables, both in the 
general framework of MSSM l28ll B l286j . or in a sp ecific SUSY breaking scenario, e.g., 



mSUGRA@, H, H, H, EJ, or GMSB[283]. The results from different group 



differ slightly, depending on the choices of the set of precision observables that are in- 
cluded in the fit, the order of loop corrections, and the experimental values that are used. 
We will review the general feature of SUSY contributions to the precision observables in 
this section, and leave the discussion of the global fit in a particular SUSY framework to 
section IV11DI 

When either the supersymmetric scalars (squarks and sleptons) or the supersymmetric 
fermions (charginos and neutralinos) are sufficiently heavy, the radiative corrections to 
precision observables are dominated by the universal gauge boson propa gator corrections, or 



the oblique corrections S, T and U [see Eqs. ( l86]) j. The authors of Ref. [2811 ] systematically 
studied the MSSM contributions to the oblique parameters from four different sectors: 
squarks, sleptons, neutralinos/charginos, and the Higgs sector. They found that relatively 
light squarks/sleptons generally make the fit to the electroweak data worse than the SM 
fit. The squark sector contributes essentially to the positive T direction. The slepton 
sector contributes negatively to S, but T remains constant or slightly positive for large 
tan/?. Both tend to be disfavored by data. The contributions from light charginos and 
neutralinos make both S and T negative, which slightly improves the fit. The best fit 
is obtained when the lightest chargino mass is near its experimental lower bound. The 
contributions from the MSSM Higgs sector are found to be small when the light CP-even 
MSSM Higgs mass is taken to be the SM Higgs boson mass. 

When both the supersymmetric scalars and fermions are light, additional non-oblique 
(vertex, external leg, and box graph) corrections to all the Zff vertices become important. 
In addition to considering the non-oblique corrections to specific processes, one must also 
include the the MSSM contributions to the muon decay amplitude since we express the 
coupling g 2 in terms of (see Section UTTj) . Morevoer, for large tan/3, the bottom and tau 
Yukawa couplings are large. In this case, the MSSM Higgs boson loops can appreciably 
affect the Zbb, Ztt and Zv T v T vertices, especially for small tua- As a result, the overall fit 
to the precision observables for large tan (3 is worse than the SM. On the other hand, the 
order a s SUSY-QCD contributions to Zqq vertex via gluino- squark- squark loop are found 
to be negligibly small when the mass for gluino and squarks are bigger than about 200 
GeV. The electroweak contributions to Zqq vertices due to light squarks and light neu- 
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tralinos/charginos are insignificant when their masses are above the current direct search 
limit. When the masses of left-handed slepton and neutralinos/charginos are light, the Zll 
vertices as well as the muon-decay amplitude are affected significantly. The fit is improved 
slightly when the left-handed slepton mass is around 200-500 GeV H- In contrast, 
the right-handed slepton and squark masses are not constrained significantly, and hence 
smaller masses are still allowed. 

In summary, for SUSY spectrum with light left-handed sleptons and 
chargino/neutralinos, the global fit in the MSSM has a lower \ 2 value than in the 
SM. Since the MSSM fit has fewer degrees of freedom than the SM fit 35 , the overall fit 
probability in the MSSM is very similar to that in the SM. 



D. Global Analysis in mSUGRA and GMSB 

In a particular SUSY breaking scenario such as mSUGRA or GMSB, the complete SUSY 
spectrum can be determined from the RGE running of only a few parameters from the high 
energy scale down to the weak scale. Due to the relatively small numbers of parameters 
in these scenarios, a global fit to the precision observables can be used to exclude certain 
region o f the parameter space if the fit is significantly worse than the SM. 

Ref. [283) studied the global fit of precision observables (including certain low energy 
measurements) in the framework of mSUGRA and GMSB. It is shown that significant 
portions of the parameter spaces of mSUGRA and GMSB are excluded. Requiring Xmssm - 
Xsm < 3.84, a lower limit on the mass of the light CP-even Higgs: rrih > 78 GeV can be 
obtained. Also, the first and second generation squark masses are constrained to be above 
280(325) GeV in the mSUGRA(GMSB) model. 

The global fit in mSUGRA performed in Ref. [282] [including (g — 2)^ and b — > 57] 



showed that at 95% C.L., the value of tan/3 is constrained to be above 6.5, while the value 
of the gaugino masses at the GUT scale has to be above ~ 220 GeV, which corresponds 
to a lower limit on the lightest neutralino(chargino) of 95 (175) GeV. 

Several analyses focusing on th e co l d da r k m atter r elic d ensity favorite region in 
mSUGRA have been performed @, H, H, Ref. H includes all the high Q 



2 



precision observables, muon (g — 2) M and b — > S7. The analysis shows a favored region 
for /i > and small m and m 1 / 2 - Ref. ((| focuses on the SUSY contributions to the W 
boson mass M\y, the effective leptonic weak mixing angle sin 2 6 e e, the anomalous magnetic 
moment of the muon (g — 2)^, and b — > 57. Higher order loop corrections are included and 
both theoretical and experimental errors are treated. A fit to these precision quantities in 
mSUGRA shows a clear preference for relative small value of mi/2, with a best-fit value of 
about 300 GeV for tan/3 = 10. An upper bound o f abo ut 600 GeV on m 1( / 2 is obtained at 
90% C.L.. We note that only the fits of Refs. 0, [287] take into account the most recent 



(g — 2)^ results (for the final report of the Brookhaven E821 experiment, see Ref. |3| 



35 The increase in the number of parameters in the MSSM compared to the SM reduces the number of 
fitting degrees of freedom. 



99 



VIII. THE EXPERIMENTAL LIMIT ON THE MSSM NEUTRAL HIGGSES 



As emphasized in Section IVIt the Higgs potential plays a crucial role in determining 
the viability of electroweak baryogenesis. Both the shape of the potential and the mass 
of the lightest, Standard Model-like Higgs boson are important in this respect. Here, we 
review what is known about the neutral Higgs bosons in the MSSM and highlight the 
various assumptions associated with the corresponding Higgs mass limits. Importantly, 
some scenarios allow for weaker experimental lower bounds on mifi than the bound for the 
mass of the SM Higgs that makes it too heavy to accommodate a strong first order EWPT 
in the SM. 

There have been extensive searches for Higgs bosons at the LEP. No signals have been 
found so far and a lower limit of 114.4 GeV has been set for the mass of the SM model 
Higgs boson at 95% C.L. H . In MSSM, two Higgs doublets need to be introduced and 
there are five physical Higgses in the spectrum: three neutral ones and two charged ones. 
When CP is conserved in the Higgs sector, the three neutral Higgses are CP eigenstates: 
two CP-even ones h° and H° and one CP-odd one A . However, there is no reason to 
exclude the CP violation in the Higgs sector. In particular, CP violation in MSSM could 
provide o ne of the i ngredients to explain the observed matter-antimatter asymmetry in the 



Universe [1981 . |199| . as discussed above. In the CP violating scenario, the three neutral 
Higgs mass eigenstates Hi, H 2 and H 3 are a mixture of the CP-even and the CP-odd 
states. The Higgs production and decay might differ significantly from the CP conserved 
scenario Analysis of t he LEP Higgs searches have been performed in both scenarios 
293 . 294 . 295 . 296l . 297 ]. In this section, we will briefly review the Higgs searches and the 



exclusion bounds on the neutral Higgs masses and other parameters. 

In the CP conserving scenario, the main production processes of h°, H° and A at 
the LEP are Higgsstraahlung e + e~ — > h°Z and e + e~ — > H°Z (if kinematically possible) 
and pair production e + e _ — > h°A° and e + e~ — ► H°A° (if kinematically possible). These 
two processes are complementary since a^z is governed by the hPZZ coupling, which 
is proportional to sin(/3 — a), while <Jh°A° is governed by the Zh°A° coupling, which is 
proportional to cos(/3 — a). For the Heavy Higgs H°, sin and cos are e xcha nged. 

The light neutral Higgs h°, whose mass is typically below 140 GeV |298| . decays domi- 
nantly into fermion pairs since its mass is below the threshold of WW and ZZ. Although, 
for particular choices of parameters, the fermionic decay may be strongly suppressed, which 
will be discussed below in the large fi, gluophobic and small a e ff benchmark models |299l . 
For the CP-odd state A , it also dominantly decays into fermion pairs, since its coupling 
to the gauge boson vanishes at tree level. For not too small tan j3, the Higgses decay into 
bb or t + t~ , while for tan/? < 1, decays to cc might be important. 

In the CP violating scenario, the mass eigenstates Hi, (i = 1, 2, 3) are a mixture of CP 
eigenstates. Each of them can be produced by Higgsstrahlung e + e~ — > H{Z via the CP- 
even field components, and also in pairs e + e~ — > HiHj(i ^ j). The relative rates depend 
on the CP-even/odd mixing. Such mixing does not occur in the tree-level potential, but 
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does appear at one-loop order. The degree of CP-mixing is proportional to the quantity 

m\ Im(/M) 



v 2 m 2 



(274) 



that arises from stop loops. Large CP violation in the Higgs sector is expected for small 
rh and large Im(fiA). The CP violation effects are also very sensitive to the precise value 
of the top quark mass, which is known with a few GeV experimental error. 

The Higgs searches in the CP-violating scenario are, in general, more challenging than 
the CP-conserving case. The reason is that the production of the lightest Higgs Hi is 
reduced due to its suppressed coupling to Z. While the production of the heavier Higgses 
are suppressed or forbidden due to kinematics. The decay of the Higgses in the CP- violating 
scenario is very similar to the CP-conserving case that discussed above. 

Higgs searches have been performed at L EP up to the highe st LEP energy of 209 GeV, 
carried out by the four LEP collaborations |294 . 1295 . 296, 2971. The searches include the 



Higgsstrahlung process and pair production process, which are sensitive over the accessible 
MSSM parameters due to their complementarity. For the Higgsstrahlung process, the 
principle signal topologies are Higgs decays to fermion pairs bb, t + t~ or flavor independent 
qq, while the Z decays into pair of jets (qq), vv (associated with missing energy), or lepton 
pairs e + e~, t + t~. The reconstruction of the Z mass offers a discrimination of 

the signal over the background. Searches including Higgs cascade decay e + e~ — > H 2 Z — > 
(H\Hi)Z have also been performed, which might play an important role when this decay 
mode is open. For Higgs pair production process, b pairs and r pairs final states have been 
studied when Higgs masses are above the t + t~ threshold. When the bb decay mode of the 
Higgs is suppressed in certain parameter spaces, flavor-independent searches are used as a 
supplementation or replacement. 



The combined LEP data show no significant signal for Higgs boson production |293 



Therefore, the search results are used to set an upper limit on the Higgs production cross 
secti on, and they are interpreted in a set of representative MSSM "benchmark" models 
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For the CP-conserving scenario, the m^-max benchmark [299| occurs when the stop 
mixing parameter is set to a large value, X t = A— ft cot (3 = 2Msusy- This model is designed 
to maximize the theoretical upper bound on m° h for a given tan (3, thereby providing the 
largest parameter space and therefore the mos t conservative exclusion limits among all the 
CP-conserving scenarios studied. Fig. fT9l |29~3j shows the excluded region in (m", m^) (left 
plot) and in (m°,tan/?) (right plot). For tan/5 < 5, the 95% C.L. exclusion bound on the 
Higgs mass is about 114 GeV, provided by the Higgsstrahlung process, while for higher 
values of tan/?, pair production process dominates and the bounds is about 93 GeV for 
both m\ and m° A . A certain region of tan/5 between 0.5 and 3 is also excluded, which, 
however, depends on the precise value of the top quark mass. The excluded region gets 
smaller for larger m t , and no lim it ca n be set for tan/3 for m t > 183 GeV. 



For the no-mixing benchmark 299], the stop mixing parameter X t is set to zero, thereby 
minimizing the stop contribution to the Higgs mass. The theoretical bounds of the pa- 
rameter space are more restrictive than in the m^-max case, although the experimental 
bounds are similar. It is worth mentioning that a small domain at m° h w 80 GeV, m° A < 3 
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FIG. 19: The MSSM exclusions, at 95% CL. (light-green) and 99.7% C.L. (dark-green), for 
the mh -max benchmark scenario, with mt=174.3 GeV. The theoretically inaccessible regions are 
shown in yellow. The dashed lines indicate the boundaries of the regions expected to be excluded 
on the basis of the Monte Carlo simulations with no signal. In the right plot, the upper edge of 
the parameter space is indicated for various top quark masses; fr om le ft to right: m t = 169.3, 
174.3, 179.3 and 183.0 GeV. The figures are reproduced from Ref. 2931 ] with kind permission of 
Springer Science and Business Media. 



GeV and tan/? < 0.7 is still allowed. This domain is not covered by the current searches 
since the branching ratio of h° — > bb is suppressed while A — > t + t~ is not kinematically 
allowed. 

In the large- fj, scenario the experimental detection is a priori challenging due to 

the suppressed decay of h° — > bb, t + t~. The dominant decay modes are h° — > cc, gg 
and W + W~ . The flavor- and decay- mode-independent searches are used instead, which 
exclude almost all of the a ccess ible MSSM parameter space. 

The gluphobic scenario [2991 ] is constructed so that Higgs-gluon-gluon coupling is sup- 
pressed, leading to a reduced Higgs production by gluon fusion at the LHC. The small a e ff 
scenario refers to the case when h° — > bb and r + r~ are suppressed, since the corresponding 
couplings are proportional to a e ff- Note that a e ff is the effective mixing angle of the neu- 
tral CP-even Higgs sector (defined in Eq. (HOI in Sec. [IT]) including radiative corrections. 
Both scenarios were invented to test situations that might be problematic at LHC. In both 
case large parameter spaces are excluded by the L EP s earches. 

The parameters of the CP- violating benchmark have been chosen to maximize the 
difference with respect t o the CP-conserving scenario: rh = 500 GeV, // = 2000 GeV and 
arg(/xA) = 90°. Fig. [20] [293 ] shows the excluded region in (m^, mjj 2 ) (left plot) and in 
{tuhx-, tan (5) (right plot). For large m# 2 , the Hi is almost completely CP-even and the 95% 
C.L. exclusion bound on the H\ mass is about 113 GeV. For lighter mass of m# 2 < 130 
GeV, Hi has a large CP-odd mixture, leading to unexcluded domain. For tan j3 between 
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FIG. 20: Exclusions, at 95% CL. (light-green) and 99.7% C.L. (dark-green), for the CP-violating 
scenario, with mt=174.3 GeV. The theoretically inaccessible regions are shown in yellow. The 
dashed lines indicate the boundaries of the regions expected to be excluded, at the 95% C.L., 
on t he ba sis of the Monte Carlo simulations with no signal. The figures are reproduced from 
Ref. H with kind permission of Springer Science and Business Media. 

about 3.5 and 10, the exclusion is particularly weak. Nonetheless, the region of < 114 
GeV and tan/5 < 3.0 are excluded by the data. Furthermore, at 95% C.L. tan/5 < 2.6 is 
excluded for all values of Higgs masses. 

The exclusion region for the Higgs mass, however, depends strongly on the top quark 
mass. For the limits discussed above, m t = 174.3 GeV was used. The exclusion power is 
reduced for larger top quark mass, especially in the region of tan/3 between 4 and 10. The 
bound on tan j3 quoted above, however, is barely sensitive to the precise value of tan (3. The 
exclusion region also depends on the CP- violating phase, arg(/xA), the \x parameter, and 
fh. The exclusion region is somewhat larger for values deviate away from the benchmark 
parameters. 

IX. CONCLUSIONS AND OUTLOOK 

Precision measurements of electroweak observables played an important role in devel- 
oping and testing the Standard Model, and they will undoubtedly be a crucial tool in 
determining the larger framework in which the SM lies. If that framework includes low- 
energy supersymmetry, then one would expect a rich array of effects to be discernible in 
precision measurements carried out at both high and low-energies. In this review, we have 
concentrated on the low-energy domain, where the "precision frontier" will lie at least until 
the era of a future e + e~ linear collider. We hope to have demonstrated that through pre- 
cise measurements of both SM observables as well as those forbidden or highly-suppressed 
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in the SM, studies of these low-energy observables will offer important information about 
SUSY that can complement what we may learn from the Large Hadron Collider. 

We also hope to have illustrated the opportunities and challenges in this field. Experi- 
mentally, recent advances have made the prospects for carrying out ~ 0.1% measurements 
of SM electroweak observables - as needed to probe SUSY - quite realistic, and a number 
of efforts are underway with such precision as a goal. Recent theoretical developments 
have also made it possible to interpret measurements at this level in terms of SUSY, as 
a number of strong interaction uncertainties have been circumvented or reduced. In both 
cases, going beyond the ~ 0.1% precision level for SM observables represents the next hori- 
zon, one that both experimentalists and theorists are beginning to approach. At the same 
time, the prospects for performing measurements of rare and forbidden observables, such 
as electric dipole moments and lepton flavor violating effects, have improved dramatically, 
with several orders of magnitude increases in sensitivity now within reach. As we hope 
to have shown, the "physics reach" of such experiments can match and even exceed what 
will be achievable at both the LHC and a linear collider, assuring that they will remain 
important avenues of study well into the future collider era. 

Finally, we emphasize that though we have concentrated here on the Minimal Super- 
symmetric Standard Model, it is by no means assured that the MSSM (together with the 
conventional assumptions about SUSY-breaking mediation) is the right low-energy man- 
ifestation of SUSY. There remains a wealth of possible variations and a correspondingly 
rich field of low-energy, precision electroweak phenomenology to be explored. This oppor- 
tunity, and the attendant experimental and theoretical challenges, will surely keep particle, 
nuclear, and atomic physicists busy for many years. 
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